(377j) Bridging Two-Liquid Theory with Molecular Simulations for Electrolytes: An Investigation of Aqueous NaCl Solution
The eNRTL model requires two adjustable binary interaction parameters per molecule-electrolyte pair (Ïca,m and Ïm,ca), commonly identified from regression of experimental data. However, quantifying the binary interaction parameters from regression may not always be practical since it requires availability of experimental data, substantial experience in data treatment, and necessity of properly defining the optimization problem. Recently, Ravichandran et al.  presented a novel methodology to predict the NRTL binary parameters from molecular dynamics (MD) simulations for nonelectrolyte solutions. Hence, it is of great interest to explore the feasibility of quantifying the binary interaction parameters of the eNRTL model from MD simulations for aqueous electrolytes.
In this study, we introduce a methodology to quantify the eNRTL binary interaction parameters for molecule-electrolyte pairs, (Ïca,m and Ïm,ca), by combining the statistical mechanic framework of two-liquid theory for electrolytes with MD simulations. Ï parameters are formulated as a function of species diameters (Ï), effective interaction strength (â), and the first neighbor shell radii around the center species (R), which are calculated from the MD simulations. Aqueous NaCl solution is selected for the case study due to availability of proven force field parameters  and its ubiquitous presence in nature. To investigate the composition dependence of the binary interaction parameters, the simulations are carried out in the concentration range of 0.11 to 16 mol.kg-1 NaCl, at 298.15 K and 1 bar. Furthermore, we obtain non-randomness factor parameter (Î±) from the MD simulations following the definition introduced by Brandani and Prausnitz  from the total coordination number (z) around a center species. The binary interaction parameters determined from the MD simulations are consistent with those obtained from regression. Future studies will investigate the applicability of the proposed approach for a broader range of electrolytes such as potassium and magnesium chlorides.
 S. Honarparvar, S. H. Saravi, D. Reible, and C.-C. Chen, "Comprehensive thermodynamic modeling of saline water with electrolyte NRTL model: A study on aqueous Ba2+-Na+-Clâ-SO42â quaternary system," Fluid Phase Equilibria, vol. 447, pp. 29-38, 2017.
 D. M. Sherman and M. D. Collings, "Ion association in concentrated NaCl brines from ambient to supercritical conditions: results from classical molecular dynamics simulations," Geochemical Transactions, vol. 3, pp. 102-107, 2002.
 W. Habermann and E. H. Pommer, "Biological fuel cells with sulphide storage capacity," Applied microbiology and biotechnology, vol. 35, pp. 128-133, 1991.
 K. S. Pitzer, âThermodynamics of electrolytes. I. Theoretical basis and general equations,â Journal of Physical Chemistry, vol. 77, pp. 268-277, 1973
 Y. Song and C.-C. Chen, "Symmetric electrolyte nonrandom two-liquid activity coefficient model," Industrial & Engineering Chemistry Research, vol. 48, pp. 7788-7797, 2009.
 N. Hossain, A. Ravichandran, R. Khare, and C.-C. Chen, âRevisiting Electrolyte Thermodynamic Models: Insight from Molecular Simulationsâ, submitted to AIChE Journal, 2018.
 A. Ravichandran, R. Khare, and C. C. Chen, "Predicting NRTL binary interaction parameters from molecular simulations," AIChE Journal, 2018. https://doi.org/10.1002/aic.16117
 V. Brandani and J. M. Prausnitz, "Two-fluid theory and thermodynamic properties of liquid mixtures: General theory," Proceedings of the National Academy of Sciences, vol. 79, pp. 4506-4509, 1982.