On the calculation of solubilities via direct coexistence simulations: Investigation of NaCl aqueous solutions and Lennard-Jones binary mixtures

Espinosa, J. R.; Young, J. M.; Jiang, H.; Gupta, D.; Vega, C.; Sanz, E.; Debenedetti, P. G.; Panagiotopoulos, A. Z.

doi:10.1063/1.4964725

Title: On the calculation of solubilities via direct coexistence simulations: Investigation of NaCl aqueous solutions and Lennard-Jones binary mixtures

Journal Article · Wed Oct 19 00:00:00 EDT 2016 · Journal of Chemical Physics

DOI:https://doi.org/10.1063/1.4964725· OSTI ID:1465107

Espinosa, J. R. ^[1];

^[2];

^[2]; Vega, C. ^[1]; Sanz, E. ^[1];

^[2];

^[2]

Univ. Complutense de Madrid, Madrid (Spain). Dept. de Quimica Fisica I, Facultad de Ciencias Quimicas
Princeton Univ., NJ (United States). Dept. of Chemical and Biological Engineering

Direct coexistence molecular dynamics simulations of NaCl solutions and Lennard-Jones binary mixtures were performed to explore the origin of reported discrepancies between solubilities obtained by direct interfacial simulations and values obtained from the chemical potentials of the crystal and solution phases. We find that the key cause of these discrepancies is the use of crystal slabs of insufficient width to eliminate finite-size effects. We observe that for NaCl crystal slabs thicker than 4 nm (in the direction perpendicular to the interface), the same solubility values are obtained from the direct coexistence and chemical potential routes, namely, 3.7 ± 0.2 molal at T = 298.15 K and p = 1 bar for the JC-SPC/E model. Such finite-size effects are absent in the Lennard-Jones system and are likely caused by surface dipoles present in the salt crystals. We confirmed that μs-long molecular dynamics runs are required to obtain reliable solubility values from direct coexistence calculations, provided that the initial solution conditions are near the equilibrium solubility values; even longer runs are needed for equilibration of significantly different concentrations. We do not observe any effects of the exposed crystal face on the solubility values or equilibration times. For both the NaCl and Lennard-Jones systems, the use of a spherical crystallite embedded in the solution leads to significantly higher apparent solubility values relative to the flat-interface direct coexistence calculations and the chemical potential values. Finally, our results have broad implications for the determination of solubilities of molecular models of ionic systems

View Accepted Manuscript (DOE)

View Accepted Manuscript (Publisher)

Cite

Export

Save

Research Organization:: Princeton Univ., NJ (United States)

Sponsoring Organization:: USDOE Office of Science (SC)

Grant/Contract Number:: SC0002128

OSTI ID:: 1465107

Alternate ID(s):: OSTI ID: 1329331

Journal Information:: Journal of Chemical Physics, Vol. 145, Issue 15; ISSN 0021-9606

Publisher:: American Institute of Physics (AIP)Copyright Statement

Country of Publication:: United States

Language:: English

Citation Metrics:

Cited by: 74 works

Citation information provided by
Web of Science

References (50)

Consensus on the solubility of NaCl in water from computer simulations using the chemical potential route Benavides, A. L.; Aragones, J. L.; Vega, C. The Journal of Chemical Physics, Vol. 144, Issue 12 https://doi.org/10.1063/1.4943780	journal	March 2016
A molecular dynamics method for simulations in the canonical ensemble Nosé, Shūichi Molecular Physics, Vol. 52, Issue 2 https://doi.org/10.1080/00268978400101201	journal	June 1984
Nucleation of NaCl from Aqueous Solution: Critical Sizes, Ion-Attachment Kinetics, and Rates Zimmermann, Nils E. R.; Vorselaars, Bart; Quigley, David Journal of the American Chemical Society, Vol. 137, Issue 41 https://doi.org/10.1021/jacs.5b08098	journal	October 2015
Determining the three-phase coexistence line in methane hydrates using computer simulations Conde, M. M.; Vega, C. The Journal of Chemical Physics, Vol. 133, Issue 6 https://doi.org/10.1063/1.3466751	journal	August 2010
Mean spherical model for asymmetric electrolytes: I. Method of solution Blum, L. Molecular Physics, Vol. 30, Issue 5 https://doi.org/10.1080/00268977500103051	journal	November 1975
LINCS: A linear constraint solver for molecular simulations Hess, Berk; Bekker, Henk; Berendsen, Herman J. C. Journal of Computational Chemistry, Vol. 18, Issue 12 https://doi.org/10.1002/(SICI)1096-987X(199709)18:12<1463::AID-JCC4>3.0.CO;2-H	journal	September 1997
Computer simulations of NaCl association in polarizable water Smith, David E.; Dang, Liem X. The Journal of Chemical Physics, Vol. 100, Issue 5 https://doi.org/10.1063/1.466363	journal	March 1994
A method for computing the solubility limit of solids: Application to sodium chloride in water and alcohols Paluch, Andrew S.; Jayaraman, Saivenkataraman; Shah, Jindal K. The Journal of Chemical Physics, Vol. 133, Issue 12 https://doi.org/10.1063/1.3478539	journal	September 2010
Calculation of the melting point of alkali halides by means of computer simulations Aragones, J. L.; Sanz, E.; Valeriani, C. The Journal of Chemical Physics, Vol. 137, Issue 10 https://doi.org/10.1063/1.4745205	journal	September 2012
A Kirkwood–Buff derived force field for sodium chloride in water Weerasinghe, Samantha; Smith, Paul E. The Journal of Chemical Physics, Vol. 119, Issue 21 https://doi.org/10.1063/1.1622372	journal	December 2003
Chemical Potentials, Activity Coefficients, and Solubility in Aqueous NaCl Solutions: Prediction by Polarizable Force Fields Moučka, Filip; Nezbeda, Ivo; Smith, William R. Journal of Chemical Theory and Computation, Vol. 11, Issue 4 https://doi.org/10.1021/acs.jctc.5b00018	journal	March 2015
Direct Coexistence Methods to Determine the Solubility of Salts in Water from Numerical Simulations. Test Case NaCl Manzanilla-Granados, Héctor M.; Saint-Martín, Humberto; Fuentes-Azcatl, Raúl The Journal of Physical Chemistry B, Vol. 119, Issue 26 https://doi.org/10.1021/acs.jpcb.5b00740	journal	June 2015
Rate of homogeneous crystal nucleation in molten NaCl Valeriani, C.; Sanz, E.; Frenkel, D. The Journal of Chemical Physics, Vol. 122, Issue 19 https://doi.org/10.1063/1.1896348	journal	May 2005
Computing the free energy of molecular solids by the Einstein molecule approach: Ices XIII and XIV, hard-dumbbells and a patchy model of proteins Noya, E. G.; Conde, M. M.; Vega, C. The Journal of Chemical Physics, Vol. 129, Issue 10 https://doi.org/10.1063/1.2971188	journal	September 2008
GROMACS 3.0: a package for molecular simulation and trajectory analysis Lindahl, Erik; Hess, Berk; van der Spoel, David Journal of Molecular Modeling, Vol. 7, Issue 8 https://doi.org/10.1007/s008940100045	journal	August 2001
Charged hard spheres in dipolar hard sphere solvents. A model for electrolyte solutions Levesque, D.; Weis, J. J.; Patey, G. N. The Journal of Chemical Physics, Vol. 72, Issue 3 https://doi.org/10.1063/1.439333	journal	February 1980
The crystal-fluid interfacial free energy and nucleation rate of NaCl from different simulation methods Espinosa, Jorge R.; Vega, Carlos; Valeriani, Chantal The Journal of Chemical Physics, Vol. 142, Issue 19 https://doi.org/10.1063/1.4921185	journal	May 2015
Water activity as the determinant for homogeneous ice nucleation in aqueous solutions Koop, Thomas; Luo, Beiping; Tsias, Athanasios Nature, Vol. 406, Issue 6796 https://doi.org/10.1038/35020537	journal	August 2000
Kinetic aspects of the thermostatted growth of ice from supercooled water in simulations Weiss, Volker C.; Rullich, Markus; Köhler, Christof The Journal of Chemical Physics, Vol. 135, Issue 3 https://doi.org/10.1063/1.3609768	journal	July 2011
Prediction of the phase equilibria of methane hydrates using the direct phase coexistence methodology Michalis, Vasileios K.; Costandy, Joseph; Tsimpanogiannis, Ioannis N. The Journal of Chemical Physics, Vol. 142, Issue 4 https://doi.org/10.1063/1.4905572	journal	January 2015
Molecular dynamics study of salt–solution interface: Solubility and surface charge of salt in water Kobayashi, Kazuya; Liang, Yunfeng; Sakka, Tetsuo The Journal of Chemical Physics, Vol. 140, Issue 14 https://doi.org/10.1063/1.4870417	journal	April 2014
Mean ionic activity coefficients in aqueous NaCl solutions from molecular dynamics simulations Mester, Zoltan; Panagiotopoulos, Athanassios Z. The Journal of Chemical Physics, Vol. 142, Issue 4 https://doi.org/10.1063/1.4906320	journal	January 2015
Polymorphic transitions in single crystals: A new molecular dynamics method Parrinello, M.; Rahman, A. Journal of Applied Physics, Vol. 52, Issue 12 https://doi.org/10.1063/1.328693	journal	December 1981
Recent progress in molecular simulation of aqueous electrolytes: force fields, chemical potentials and solubility Nezbeda, Ivo; Moučka, Filip; Smith, William R. Molecular Physics, Vol. 114, Issue 11 https://doi.org/10.1080/00268976.2016.1165296	journal	March 2016
Determination of the melting point of hard spheres from direct coexistence simulation methods Noya, Eva G.; Vega, Carlos; de Miguel, Enrique The Journal of Chemical Physics, Vol. 128, Issue 15 https://doi.org/10.1063/1.2901172	journal	April 2008
Simulation of crystal growth from Lennard-Jones solutions Huitema, H. E. A.; van Hengstum, B.; van der Eerden, J. P. The Journal of Chemical Physics, Vol. 111, Issue 22 https://doi.org/10.1063/1.480374	journal	December 1999
Pseudo hard-sphere potential for use in continuous molecular-dynamics simulation of spherical and chain molecules Jover, J.; Haslam, A. J.; Galindo, A. The Journal of Chemical Physics, Vol. 137, Issue 14 https://doi.org/10.1063/1.4754275	journal	October 2012
Melting point and phase diagram of methanol as obtained from computer simulations of the OPLS model Gonzalez Salgado, D.; Vega, C. The Journal of Chemical Physics, Vol. 132, Issue 9 https://doi.org/10.1063/1.3328667	journal	March 2010
Simulation of nucleation in almost hard-sphere colloids: The discrepancy between experiment and simulation persists Filion, L.; Ni, R.; Frenkel, D. The Journal of Chemical Physics, Vol. 134, Issue 13 https://doi.org/10.1063/1.3572059	journal	April 2011
A new polarizable force field for alkali and halide ions Kiss, Péter T.; Baranyai, András The Journal of Chemical Physics, Vol. 141, Issue 11 https://doi.org/10.1063/1.4895129	journal	September 2014
Temperature dependence of crystal growth of hexagonal ice (Ih) Rozmanov, Dmitri; Kusalik, Peter G. Physical Chemistry Chemical Physics, Vol. 13, Issue 34 https://doi.org/10.1039/c1cp21210a	journal	January 2011
Simulation of homogeneous crystal nucleation close to coexistence ten Wolde, Pieter-Rein; Ruiz-Montero, Maria J.; Frenkel, Daan Faraday Discussions, Vol. 104 https://doi.org/10.1039/fd9960400093	journal	January 1996
Molecular Dynamics Simulations of the Dynamic and Energetic Properties of Alkali and Halide Ions Using Water-Model-Specific Ion Parameters Joung, In Suk; Cheatham, Thomas E. The Journal of Physical Chemistry B, Vol. 113, Issue 40 https://doi.org/10.1021/jp902584c	journal	October 2009
Molecular dynamic studies of the solubility of sodium chloride: fast calculations using seed crystalline cluster probe Wiebe, Heather; Louwersheimer, Johan; Weinberg, Noham Molecular Physics, Vol. 113, Issue 21 https://doi.org/10.1080/00268976.2015.1011247	journal	February 2015
Structural transformation in supercooled water controls the crystallization rate of ice Moore, Emily B.; Molinero, Valeria Nature, Vol. 479, Issue 7374 https://doi.org/10.1038/nature10586	journal	November 2011
On fluid-solid direct coexistence simulations: The pseudo-hard sphere model Espinosa, Jorge R.; Sanz, Eduardo; Valeriani, Chantal The Journal of Chemical Physics, Vol. 139, Issue 14 https://doi.org/10.1063/1.4823499	journal	October 2013
The melting lines of model systems calculated from coexistence simulations Morris, James R.; Song, Xueyu The Journal of Chemical Physics, Vol. 116, Issue 21 https://doi.org/10.1063/1.1474581	journal	June 2002
Solubility of NaCl in water by molecular simulation revisited Aragones, J. L.; Sanz, E.; Vega, C. The Journal of Chemical Physics, Vol. 136, Issue 24 https://doi.org/10.1063/1.4728163	journal	June 2012
Molecular force fields for aqueous electrolytes: SPC/E-compatible charged LJ sphere models and their limitations Moučka, Filip; Nezbeda, Ivo; Smith, William R. The Journal of Chemical Physics, Vol. 138, Issue 15 https://doi.org/10.1063/1.4801322	journal	April 2013
Thermodynamic and Transport Properties of H ₂ O + NaCl from Polarizable Force Fields Jiang, Hao; Mester, Zoltan; Moultos, Othonas A. Journal of Chemical Theory and Computation, Vol. 11, Issue 8 https://doi.org/10.1021/acs.jctc.5b00421	journal	July 2015
A smooth particle mesh Ewald method Essmann, Ulrich; Perera, Lalith; Berkowitz, Max L. The Journal of Chemical Physics, Vol. 103, Issue 19 https://doi.org/10.1063/1.470117	journal	November 1995
Exact Solution of an Integral Equation for the Structure of a Primitive Model of Electrolytes Waisman, Eduardo; Lebowitz, J. L. The Journal of Chemical Physics, Vol. 52, Issue 8 https://doi.org/10.1063/1.1673642	journal	April 1970
Solubility of KF and NaCl in water by molecular simulation Sanz, E.; Vega, C. The Journal of Chemical Physics, Vol. 126, Issue 1 https://doi.org/10.1063/1.2397683	journal	January 2007
The missing term in effective pair potentials Berendsen, H. J. C.; Grigera, J. R.; Straatsma, T. P. The Journal of Physical Chemistry, Vol. 91, Issue 24 https://doi.org/10.1021/j100308a038	journal	November 1987
Molecular Force Field Development for Aqueous Electrolytes: 1. Incorporating Appropriate Experimental Data and the Inadequacy of Simple Electrolyte Force Fields Based on Lennard-Jones and Point Charge Interactions with Lorentz–Berthelot Rules Moučka, Filip; Nezbeda, Ivo; Smith, William R. Journal of Chemical Theory and Computation, Vol. 9, Issue 11 https://doi.org/10.1021/ct4006008	journal	October 2013
Temperature-dependent solubilities and mean ionic activity coefficients of alkali halides in water from molecular dynamics simulations Mester, Zoltan; Panagiotopoulos, Athanassios Z. The Journal of Chemical Physics, Vol. 143, Issue 4 https://doi.org/10.1063/1.4926840	journal	July 2015
A molecular dynamics method for simulations in the canonical ensemble NosÉ, SH[Ubar]ICHI Molecular Physics, Vol. 100, Issue 1 https://doi.org/10.1080/00268970110089108	journal	January 2002
Structural transformation in supercooled water controls the crystallization rate of ice Moore, Emily B.; Molinero, Valeria arXiv https://doi.org/10.48550/arxiv.1107.1622	text	January 2011
Calculation of the melting point of alkali halides by means of computer simulations Aragones, J. L.; Sanz, E.; Valeriani, C. arXiv https://doi.org/10.48550/arxiv.1209.3692	text	January 2012
Recent Progress in Molecular Simulation of Aqueous Electrolytes: Force Fields, Chemical Potentials and Solubility Nezbeda, Ivo; Moučka, Filip; Smith, William R. arXiv https://doi.org/10.48550/arxiv.1602.04065	text	January 2016

Cited By (25)

Solubility prediction from first principles: a density of states approach Boothroyd, Simon; Kerridge, Andy; Broo, Anders Physical Chemistry Chemical Physics, Vol. 20, Issue 32 https://doi.org/10.1039/c8cp01786g	journal	January 2018
Nucleation of urea from aqueous solution: Structure, critical size, and rate Mandal, Taraknath; Larson, Ronald G. The Journal of Chemical Physics, Vol. 146, Issue 13 https://doi.org/10.1063/1.4979141	journal	April 2017
Computational methodology for solubility prediction: Application to the sparingly soluble solutes Li, Lunna; Totton, Tim; Frenkel, Daan The Journal of Chemical Physics, Vol. 146, Issue 21 https://doi.org/10.1063/1.4983754	journal	June 2017
Charge compensation at the interface between the polar NaCl(111) surface and a NaCl aqueous solution Sayer, Thomas; Zhang, Chao; Sprik, Michiel The Journal of Chemical Physics, Vol. 147, Issue 10 https://doi.org/10.1063/1.4987019	journal	September 2017
A potential model for sodium chloride solutions based on the TIP4P/2005 water model Benavides, A. L.; Portillo, M. A.; Chamorro, V. C. The Journal of Chemical Physics, Vol. 147, Issue 10 https://doi.org/10.1063/1.5001190	journal	September 2017
Forward flux sampling calculation of homogeneous nucleation rates from aqueous NaCl solutions Jiang, Hao; Haji-Akbari, Amir; Debenedetti, Pablo G. The Journal of Chemical Physics, Vol. 148, Issue 4 https://doi.org/10.1063/1.5016554	journal	January 2018
NaCl nucleation from brine in seeded simulations: Sources of uncertainty in rate estimates Zimmermann, Nils. E. R.; Vorselaars, Bart; Espinosa, Jorge R. The Journal of Chemical Physics, Vol. 148, Issue 22 https://doi.org/10.1063/1.5024009	journal	June 2018
Chemical potentials of alkaline earth metal halide aqueous electrolytes and solubility of their hydrates by molecular simulation: Application to CaCl ₂ , antarcticite, and sinjarite Moučka, Filip; Kolafa, Jiří; Lísal, Martin The Journal of Chemical Physics, Vol. 148, Issue 22 https://doi.org/10.1063/1.5024212	journal	June 2018
Addressing hysteresis and slow equilibration issues in cavity-based calculation of chemical potentials Wand, C. R.; Totton, T. S.; Frenkel, D. The Journal of Chemical Physics, Vol. 149, Issue 1 https://doi.org/10.1063/1.5036963	journal	July 2018
Computational methodology for solubility prediction: Application to sparingly soluble organic/inorganic materials Li, Lunna; Totton, Tim; Frenkel, Daan The Journal of Chemical Physics, Vol. 149, Issue 5 https://doi.org/10.1063/1.5040366	journal	August 2018
Communication: Nucleation rates of supersaturated aqueous NaCl using a polarizable force field Jiang, Hao; Debenedetti, Pablo G.; Panagiotopoulos, Athanassios Z. The Journal of Chemical Physics, Vol. 149, Issue 14 https://doi.org/10.1063/1.5053652	journal	October 2018
Calculation of the water-octanol partition coefficient of cholesterol for SPC, TIP3P, and TIP4P water Espinosa, Jorge R.; Wand, Charlie R.; Vega, Carlos The Journal of Chemical Physics, Vol. 149, Issue 22 https://doi.org/10.1063/1.5054056	journal	December 2018
Nucleation in aqueous NaCl solutions shifts from 1-step to 2-step mechanism on crossing the spinodal Jiang, Hao; Debenedetti, Pablo G.; Panagiotopoulos, Athanassios Z. The Journal of Chemical Physics, Vol. 150, Issue 12 https://doi.org/10.1063/1.5084248	journal	March 2019
Solubility of paracetamol in ethanol by molecular dynamics using the extended Einstein crystal method and experiments Bellucci, Michael A.; Gobbo, Gianpaolo; Wijethunga, Tharanga K. The Journal of Chemical Physics, Vol. 150, Issue 9 https://doi.org/10.1063/1.5086706	journal	March 2019
Carbon solubility in liquid silicon: A computational analysis across empirical potentials Luo, Jinping; Alateeqi, Abdullah; Liu, Lijun The Journal of Chemical Physics, Vol. 150, Issue 14 https://doi.org/10.1063/1.5088495	journal	April 2019
Thermodynamic analysis of the stability of planar interfaces between coexisting phases and its application to supercooled water Singh, Rakesh S.; Palmer, Jeremy C.; Panagiotopoulos, Athanassios Z. The Journal of Chemical Physics, Vol. 150, Issue 22 https://doi.org/10.1063/1.5097591	journal	June 2019
Solubility prediction for a soluble organic molecule via chemical potentials from density of states Boothroyd, Simon; Anwar, Jamshed The Journal of Chemical Physics, Vol. 151, Issue 18 https://doi.org/10.1063/1.5117281	journal	November 2019
A force field of Li ⁺ , Na ⁺ , K ⁺ , Mg ²⁺ , Ca ²⁺ , Cl ⁻ , and SO42− in aqueous solution based on the TIP4P/2005 water model and scaled charges for the ions Zeron, I. M.; Abascal, J. L. F.; Vega, C. The Journal of Chemical Physics, Vol. 151, Issue 13 https://doi.org/10.1063/1.5121392	journal	October 2019
Ion dissolution mechanism and kinetics at kink sites on NaCl surfaces Joswiak, Mark N.; Doherty, Michael F.; Peters, Baron Proceedings of the National Academy of Sciences, Vol. 115, Issue 4 https://doi.org/10.1073/pnas.1713452115	journal	January 2018
Molecular dynamics simulation of the surface tension of aqueous sodium chloride: from dilute to highly supersaturated solutions and molten salt Wang, Xiaoxiang; Chen, Chuchu; Binder, Kurt Atmospheric Chemistry and Physics, Vol. 18, Issue 23 https://doi.org/10.5194/acp-18-17077-2018	journal	January 2018
Charge compensation at the interface between the polar NaCl(111) surface and a NaCl aqueous solution Sayer, Thomas; Zhang, C.; Sprik, Michiel Apollo - University of Cambridge Repository https://doi.org/10.17863/cam.14313	text	January 2017
Computational methodology for solubility prediction: Application to sparingly soluble organic/inorganic materials. Li, Lunna; Totton, Tim; Frenkel, Daan Apollo - University of Cambridge Repository https://doi.org/10.17863/cam.31809	text	January 2018
Addressing hysteresis and slow equilibration issues in cavity-based calculation of chemical potentials. Wand, Cr; Totton, Ts; Frenkel, Daan Apollo - University of Cambridge Repository https://doi.org/10.17863/cam.30174	text	January 2018
Computational methodology for solubility prediction: Application to the sparingly soluble solutes Li, L.; Totton, T.; Frenkel, Daan Apollo - University of Cambridge Repository https://doi.org/10.17863/cam.12702	text	January 2017
Molecular Dynamics Simulation of the Surface Tension of Aqueous Sodium Chloride: from Dilute to Highly Supersaturated Solutions and Molten Salt Wang, Xiaoxiang; Chen, Chuchu; Binder, Kurt Atmospheric Chemistry and Physics Discussions https://doi.org/10.5194/acp-2017-1013	journal	November 2017

Similar Records

Molecular simulation of disjoining-pressure isotherms for free liquid , Lennard-Jones thin films

Journal Article · Mon Oct 01 00:00:00 EDT 2001 · Journal of Physical Chemistry B: Materials, Surfaces, Interfaces, amp Biophysical · OSTI ID:1465107

Bhatt, Divesh; Newman, John; Radke, C J

Molecular simulation of phase coexistence: Finite-size effects and determination of critical parameters for two- and three-dimensional Lennard-Jones fluids

Journal Article · Tue Nov 01 00:00:00 EST 1994 · International Journal of Thermophysics · OSTI ID:1465107

Panagiotopoulos, A Z

Effect of dispersive long-range corrections to the pressure tensor: The vapour-liquid interfacial properties of the Lennard-Jones system revisited

Journal Article · Fri Nov 14 00:00:00 EST 2014 · Journal of Chemical Physics · OSTI ID:1465107

Martínez-Ruiz, F. J.; Mendiboure, B.; Moreno-Ventas Bravo, A. I.

Related Subjects

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY

Title: On the calculation of solubilities via direct coexistence simulations: Investigation of NaCl aqueous solutions and Lennard-Jones binary mixtures

Citation Formats

References (50)

Cited By (25)

Similar Records

Related Subjects