Mean Activity Coefficients of NaCl in the Mixture of ... - ACS Publications

Aug 1, 2017 - Department of Chemistry, University of Pamplona, (57 + 7) 5685303−5685304, IBEAR FJ-207 Biofuels Laboratory, Pamplona,. Colombia...
0 downloads 0 Views 1MB Size
Article pubs.acs.org/jced

Mean Activity Coefficients of NaCl in the Mixture of 2‑Hydroxyethylammonium Butyrate + H2O at 298.15 K Eliseo Amado-González,* Irina Lupita González-Gutierrez, and Wilfred Goméz-Jaramillo Department of Chemistry, University of Pamplona, (57 + 7) 5685303−5685304, IBEAR FJ-207 Biofuels Laboratory, Pamplona, Colombia S Supporting Information *

ABSTRACT: The main objective of this work is to continue with a series of electrochemistry studies of mixtures of ionic liquids and water as a solvent to provide accurate data for future particular applications [Amado-Gonzalez et al. J. Chem. Eng. Data 2017, 62, 752]. The mean activity coefficients for NaCl in [2-hydroxyethylammonium butyrate (2-HEAB) + H2O] as a solvent mixture were determined by cell potential measurements: Na−ion selective electrode (ISE)|NaCl (m), 2-HEAB (w), H2O (1 − w)|Cl−ion selective electrode (ISE) at molalities from 0.10 to 3.20 mol·kg−1 at 298.15 K. Different weight fractions (w) of 2-HEAB with w = 0.01, 0.05, 0.1, 0.2, 0.3, and 0.4 were used. At higher concentrations of w = 0.40, NaCl is salting out by 2HEAB in water. The Pitzer ion interaction parameters β0, β1, and Cγ were used to find the values of osmotic coefficients, solvent activity, and the excess Gibbs free energy for the mixed electrolyte system. The results may be interpreted by the clathrate-like formation of 2-HEAB + water. A qualitative description of the relation between water and the 2-HEAB was done using the general AMBER force field (GAFF).

1. INTRODUCTION Thermodynamic properties of mixtures of [ionic liquids (ILs) + H2O] are important to evaluate future applications.1 Ionic liquids (ILs) or the future new solvents can be structurally designed as proposed by Lowe and Rendall.2 Due to the physicochemical properties of ILs,3−5 they may have many different industrial uses.6 Activity coefficients and osmotic coefficients of electrolyte solutions are extremely useful to test new electrolyte solution theories.7−15 Even though activity coefficients of ammonium ILs have been reported in aqueous solutions by the isopiestic method,16−23 to our knowledge, research papers about the activity coefficients of NaCl in ILs + water are still scarce, and studies of mean ionic activity coefficients of NaCl reported in the literature for ternary liquid mixtures are still limited.24−26 Because of the high costs of ILs and import restrictions, the synthesis of new ILs should be a goal. Protic ionic liquids offer a good possibility because they can be produced by a simple acid− base neutralization reaction.27 2-Hydroxyethylammonium butyrate (2-HEAB) is therefore evaluated for its potential abilities as a solvent. In this sense, our goal is to evaluate how the thermodynamic properties of salts like NaCl would be affected by the (2-HEAB + H2O) solvent mixture. In this work, we found the activity coefficients of NaCl in the mixture (2-HEAB + H2O) by cell potential measurements at 298.15 K. The results were fitted to the Pitzer model. Then, the osmotic coefficients, solvent mixture activities, and excess Gibbs free energies of these systems were calculated. A calculation of hydrogen interaction at w = 0.4 was completed using the general AMBER force field (GAFF). © 2017 American Chemical Society

2. EXPERIMENTAL SECTION 2.1. Materials, Synthesis, and Measurements. In Table 1 the chemicals, suppliers, and stated purity are shown. The IL, 2HEAB was synthesized by acid−base neutralization reactions between ethanolamine and the organic acid as previously described.26 The reagents monoethanolamine and butyric acid (Sigma-Aldrich, 99%) are analytical (AR grade) reagents. The synthesis of 2-HEAB was done as reported by literature.28−30 Strong agitation was performed at heating at 323 K. The IL was dried under a vacuum of 20 kPa during 8 h until ΔpH (pH at the equivalent point − pH experimental) < 0.07. An ATR-FTIR spectra for 2-HEAB was done. The ammonium structure was found at 3200−2400 cm−1. The OH stretching vibration was embedded in this band. The carbonyl stretching and N−H inplane bending vibrations were at 1600 cm−1. Water content was verified by a volumetric titroline KF (Schott instruments) as 6 × 10−4 mass fraction. No phase change of ILs was observed. The analysis of 1H NMR gave the next results for the 2-HEAB: (400 MHz, D2O) δ 4.72 (s, 4H, H3b y H4b), 3.75−3.66 (m, 2H, H2b), 3.06−2.97 (m, 2H, H1b), 2.05 (t, J = 7.3 Hz, 2H, H1a), 1.52−1.40 (m, 2H, H2a), 0.79 (t, J = 7.4 Hz, 3H, H3a) and a purity higher that 99% mass fraction (details available in the Supporting Information). Received: March 18, 2017 Accepted: July 10, 2017 Published: August 1, 2017 2384

DOI: 10.1021/acs.jced.7b00278 J. Chem. Eng. Data 2017, 62, 2384−2391

Journal of Chemical & Engineering Data

Article

Table 1. Chemicals, Suppliers, and Method To Check the Stated Puritya chemical name NaCl butyric acid ethanolamine 2-HEABb

source

purity

Sigma-Aldrich, Riedel-de Haën, p.a., Reag. ACS CAS no. 7647-14-5 Merck, CAS no.107-92-6 Merck, CAS no.141-43-5 synthesis

99.8%

no further purification

99.0% 99.5% 97.0%

no further purification no further purification Fischer-Karl titration 1H RMN, ATR-FTIR spectra double distillation

water a

purification method

mass water contenta

electrical conductivity (S·cm−1)

≤0.02% ≤0.09% ≤0.06% 2-HEAB/water 0.20 > 2-HEAB/water 0.30 > 2-HEAB/water 0.05 > 2-HEAB/water 0.01 > 2-HEAB/ water 0.10. The decrease of the γ± (NaCl) may be the result of the variation in the dielectric constant of the (2-HEAB + H2O) mixture solvent from 78.08 to 72.69 where the addition of NaCl increases the electrostatic interactions. In (NaCl + 2-HEAB + H2O) mixtures, structures clathrate-like may be present for the minimum of the curves. Head-Gordon43 found the solvation of nonpolar groups by hydrogen-bonded. In Figure 3 the osmotic coefficients show a minimum of the curves varying between 0.32 and 0.34 m at w = 0.4. Then we found that the osmotic coefficient increased with the NaCl increment. One of the problems at higher w of IL is NaCl salting out. It may be possible that, after clathrate-like formation, this structure disappears as the NaCl concentration increase due to competition between the ions of NaCl and 2-HEAB. In Table 4, β0, β1, and Cγ with the standard deviation (σ) are listed. β0 and β1 are nearly dependent on w. However, Cγ is nearly independent of w. 3.3. Activity of Solvent (as) and Excess Gibbs Free Energy (GE) by the Pitzer Model. In this work, the solvent is the mixture of (2-HEAB + H2O) in weight fraction (w = 0.01, 0.05, 0.10, 0.20, 0.30, and 0.40). Following the methodology,25 the molar weight of (2-HEAB + H2O) was calculated for each weigh fraction as Ms = X1M1 + X2M2 and M1= 149.1855 g/mol and M2 = 18.0151 g/mol. In Table 5, GE and as are calculated by

Figure 4. as for the (2-HEAB + H2O) solvent mixture at T = 298.15 K.

Figure 5. Excess free Gibbs energy GE of NaCl in the (2-HEAB + H2O) solvent mixture at T = 298.15 K.

Figure 6. Hydrogen bond number along 0.5 ns of classical molecular dynamic simulations for w = 0.4.

GE = RT[vAmA (1 − ϕ + ln γA )]

(6)

N ⎡ mi ⎤ ⎥ as = exp⎢ −ϕMs ∑ ⎢⎣ 1000 ⎥⎦ i=1

(7)

Figure 4 shows that solvent activity is reduced by increasing molalities of NaCl in the ternary system. Schröder et al.44 considered that the values of GE are produced by the electrostatic interaction among the IL, anions, cations, and water. In Figure 5, as shown on Table 5, the negative value to GE increases with water at w = 0.4; this may be produced by the interaction of the OH group in the anion. Figure 6 shows that increments of the negative values of GE depend on the increase of 2-HEAB content

In eq 2, the ion-interaction parameters were evaluated by using multiple linear regression technique. The calculations of the osmotic coefficient38 (ϕ) were given by ϕ − 1 = f ϕ + mBϕ + m2C γ

∑i = 1 xivP i i

(3)

The dielectric constant of (2-HEAB + H2O) solvent mixtures were calculated by Wang and Anderko procedure.39,40 Equation 2388

DOI: 10.1021/acs.jced.7b00278 J. Chem. Eng. Data 2017, 62, 2384−2391

Journal of Chemical & Engineering Data

Article

Figure 7. Molecular dynamic snapshot along the NpT ensemble.

salting out by 2-HEAB in water. This suggests that 2-HEAB may have future uses in water desalting. The parameters of the Pitzer model β0, β1, and Cγ were calculated by the multiple linear regression method. The values of GE for NaCl in mixtures of (2HEAB + H2O) are negative. This suggests that the electrostatic interactions between the 2-HEAB and NaCl are dominant. The behavior of γ± (NaCl) and ϕ was discussed by the clathrate-like formation of (2-HEAB + H2O). At w = 0.4 the maximum number of hydrogen bonds calculated by molecular dynamic simulations is around 12 between the IL and water.

in the (2-HEAB + H2O) solvent mixture. In a previous paper, we found that ILs like 1-ethyl-3-methyl-imidazolium methanesulfonate, [Emim][MeSO3], or 1-ethyl-3-methyl-imidazolium ethyl sulfate, ([Emim][EtSO4] + H2O), have an decreasing effect of the NaCl activity coefficient.1 As proposed by the structural and electrostatic model,45 the values of GE of the system NaCl in the (2-HEAB + H2O) solvent mixture may be mainly controlled by electrostatic competitive interaction.46 3.4. Molecular Dynamics. All classical molecular dynamic simulations were performed using GROMACS 5.0.7 at w = 0.4 along 0.5 ns. Electrostatic interactions were calculated using the smooth particle mesh Ewald summation with a real space cutoff of 0.6 nm. Lennard−Jones interactions were calculated using a 12−6 potential function with a cutoff radii of 0.6 nm. Simulations were carried out with a single ionic liquid molecule surrounded by 127 TIP3 water molecules into a cubic box with an edge length of approximately 1.5 nm. For all cases, standard periodic boundary conditions in all directions were considered. At the beginning, 90 000 steps of steepest descent minimization were performed to get an energy minimum and elsewise to remove any bad intra−intermolecular bad contact. Molecular simulations were driven in three steps: (i) initial short 50 ps NVT ensemble simulation with a constant temperature of 298.15 K using the modified Berendsen thermostat with update frequencies of 0.4 ps; (ii) short 50 ps isobaric−isothermal (NpT) ensemble simulation at 298.15 K and 1 bar. The pressure control for NpT was achieved using a Parrinello−Rahman barostat with coupling time of 2 ps; (iii) long 500 ps NpT ensemble using the same above condition to final analysis. All bonds and angles were kept fixed using the LINCS algorithm.47−51 The variations of hydrogen bond numbers for IL + water suggest that the number of hydrogen bonds are higher in anion (butyrate) than in the cations as shown in Figure 6 for w = 0.40. A molecular dynamic snapshop at 125, 250, 375, and 500 ps along NpT ensemble that the IL (2-HEAB) has an effect over the structure of water as shown in Figure 7.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jced.7b00278. Data for calibration of Na+ and Cl− selective electrode pair at 298.15 K, ATR-FTIR spectra for 2-HEAB, and the 1H NMR for the 2-HEAB (PDF)



AUTHOR INFORMATION

Corresponding Author

*Tel.: +51 7 3114621948. E-mail: [email protected]. ORCID

Eliseo Amado-González: 0000-0003-4523-1323 Wilfred Goméz-Jaramillo: 0000-0003-3645-4007 Funding

This project was supported by University of Pamplona. The first author is grateful to Escola d’Engyneria, University of Santiago de Compostela, for a scholarship and Dr. Miguel Angel Iglesias for acceptance in the research group. Financial support from University of Pamplona PR400-156.012-008 (GA313-BP2015) is acknowledged. Notes

The authors declare no competing financial interest.



4. CONCLUSION The γ± of NaCl in mixtures of the (2-HEAB + H2O) solvent mixture was calculated by studied by the electrode potential method using Na-ISE and Cl-ISE at 298.15 K. The activity coefficients of NaCl in pure aqueous solution and in the (2HEAB + H2O) solvent mixture were measured by galvanic cells without liquid junction. The measurements were made at different weight fractions (w) of 2-HEAB. The addition of 2HEAB into the solvent (water/water + 2-HEAB) increases the electrostatic attraction between [2-HEA]+ and [B]− in the solvent mixture due to the decrease in the dielectric constant. The presence of NaCl has a competitive effect with 2-HEAB for water molecules. At the higher concentration of w = 0.40, NaCl is

LIST OF SYMBOLS E:cell potential GE:excess Gibbs free energy m:molality M:molecular weight Pm:polarization w:molar fraction xi:mole fraction of component

Greek Letters

αs:solvent activity ρ:density ϕ:osmotic coefficient γ:activity coefficient εr:dielectric constant

2389

DOI: 10.1021/acs.jced.7b00278 J. Chem. Eng. Data 2017, 62, 2384−2391

Journal of Chemical & Engineering Data

Article

(21) Blanco, L. H.; Amado, E.; Calvo, J. Osmotic and activity coefficients of dilute aqueous solutions of the series Me4NI to MeBu3NI at 298.15 K. Fluid Phase Equilib. 2008, 268, 90−94. (22) Amado, E.; Blanco, L. H. Osmotic and Activity Coefficients of Dilute Aqueous Solutions of Unsymmetrical Tetraalkylammonium Iodides at 298.15 K. J. Chem. Eng. Data 2012, 57, 1044−1049. (23) Golabiazar, R.; Sadeghi, R. Vapor Pressure Osmometry Determination of the Osmotic and Activity Coefficients of Dilute Aqueous Solutions of Symmetrical Tetraalkyl Ammonium Halides at 308.15 K. J. Chem. Eng. Data 2014, 59, 76−81. (24) Sirbu, F.; Iulian, O.; Catrinel Ion, A.; Ion, I. Activity Coefficients of Electrolytes in the NaCl + Na2SO4 + H2O Ternary System from Potential Difference Measurements at (298.15, 303.15, and 308.15) K. J. Chem. Eng. Data 2011, 56, 4935−4943. (25) Ghalami-Choobar, B. Thermodynamic study of the ternary mixed electrolyte (NaCl + NiCl2 + H2O) system: Application of Pitzer model with higher-order electrostatic effects. J. Chem. Thermodyn. 2011, 43, 901−907. (26) Galleguillos-Castro, H.; Hernandez-Luis, F.; Fernandez-Merida, L.; Esteso, M. Thermodynamic Study of the NaCl + Na2SO4 + H2O System by Emf Measurements at Four Temperatures. J. Solution Chem. 1999, 28, 791−806. (27) Bicak, N. A new ionic liquid: 2-hydroxy ethylammonium formate. J. Mol. Liq. 2004, 116, 37−44. (28) Pinto, R.; Silvana Mattedi, S.; Aznar, M. Synthesis and Physical Properties of Three Protic Ionic Liquids with the Ethyl-ammonium Cation. J. Chem. Eng. Trans. 2015, 43, 1165−1170. (29) Alvarez, V. H.; Dosil, N.; González-Cabaleiro, R.; Mattedi, S.; Martin-Pastor, M.; Iglesias, M.; Navaza, J. M. Brønsted ionic liquids for sustainable process: synthesis and physical properties. J. Chem. Eng. Data 2010, 55, 625−632. (30) Rocha-Pinto, R.; Santos, D.; Mattedi, S.; Aznar, M. Density, refractive index, apparent volumes and excess molar volumes of four protic ionic liquids + water at T = 298.15 and 323.15 K. Braz. J. Chem. Eng. 2015, 32, 671−682. (31) Huang, M.; Jiang, Y.; Sasisanker, P.; Driver, G.; Weingärtner, H. Static Relative Dielectric Permittivities of Ionic Liquids at 25 °C. J. Chem. Eng. Data 2011, 56, 1494−1499. (32) Chirico, R. D.; Frenkel, M.; Magee, J. W.; Diky, V.; Muzny, C. D.; Kazakov, A. F.; Kroenlein, K.; Abdulagatov, I.; Hardin, G. R.; Acree, W. E., Jr.; Brenneke, J. F.; Brown, P. L.; Cummings, P. T.; de Loos, T. W.; Friend, D. G.; Goodwin, A. R. H.; Hansen, L. D.; Haynes, W. M.; Koga, N.; Mandelis, A.; Marsh, K. N.; Mathias, P. M.; McCabe, C.; O’Connell, J. P.; Padua, A.; Rives, V.; Schick, C.; Trusler, J. P. M.; Vyazovkin, S.; Weir, R. D.; Wu, J. Improvement of quality in publication of experimental thermophysical property data: challenges, assessment tools, global implementation, and online support. J. Chem. Eng. Data 2013, 58, 2699−2716. (33) Hernández-Luis, F.; Amado-González, E.; Esteso, M. Activity coefficients of NaCl in trehalose-water and maltose-wáter mixtures at 298.15 K. Carbohydr. Res. 2003, 338, 1415−1424. (34) Case, D. A.; Cerutti, D. S.; Cheatham, T. E., III; Darden, T. A.; Duke, R. E.; Giese, T. J.; Gohlke, H.; Goetz, A. W.; Greene, D.; Homeyer, N.; Izadi, S.; Kovalenko, A.; Lee, T. S.; LeGrand, S.; Li, P.; Lin, C.; Liu, J.; Luchko, T.; Luo, R.; Mermelstein, D.; Merz, K. M.; Monard, G.; Nguyen, H.; Omelyan, I.; Onufriev, A.; Pan, F.; Qi, R.; Roe, D. R.; Roitberg, A.; Sagui, C.; Simmerling, C. L.; Botello-Smith, W. M.; Swails, J.; Walker, R. C.; Wang, J.; Wolf, R. M.; Wu, X.; Xiao, L.; York, D. M.; Kollman, P. A. AMBER 2017; University of California: San Francisco, 2017. (35) Mobley, D. L.; Chodera, J. D.; Dill, K. A. On the Use of Orientational Restraints and Symmetry Corrections in Alchemical Free Energy Calculations. J. Chem. Phys. 2006, 125, 084902. (36) Abraham, M. J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J. C.; Hess, B.; Lindahl, E. GROMACS: High Performance Molecular Simulations through Multi-Level Parallelism from Laptops to Supercomputers. SoftwareX 2015, 1−2, 19−25.

vi:molar volume



REFERENCES

(1) Amado-Gonzalez, E.; Esteso, M.; Gomez-Jaramillo, W. Mean Activity Coefficients for NaCl in the Mixtures Containing Ionic Liquids [Emim][MeSO3] + H2O and [Emim][EtSO4]+H2O at 298.15 K. J. Chem. Eng. Data 2017, 62, 752−761. (2) Lowe, B.; Rendall, H. Aqueous solutions of unsymmetrical quaternary ammonium iodides. Trans. Faraday Soc. 1971, 67, 2318− 2327. (3) Deetlefs, M.; Seddon, K. R.; Shara, M. Predicting Physical Properties of Ionic Liquids. Phys. Chem. Chem. Phys. 2006, 8, 642−649. (4) Earle, M. J.; Esperança, J. M. S. S.; Gilea, M. A.; Canongia Lopes, J. N.; Rebelo, L. P. N.; Magee, J. W.; Seddon, K. R.; Widegren, J. A. The Distillation and Volatility of Ionic Liquids. Nature 2006, 439, 831−834. (5) Kosmulski, M.; Gustafsson, J.; Rosenholm, J. B. Thermal Stability of Low Temperature Ionic Liquids Revisited. Thermochim. Acta 2004, 412, 47−53. (6) Binnemans, K. Ionic Liquid Crystals. Chem. Rev. 2005, 105, 4148− 4204. (7) Pitzer, K. S. Thermodynamics of electrolytes. I. Theoretical basis and general equations. J. Phys. Chem. 1973, 77, 268−277. (8) Harned, H. S.; Robinson, R. A. Multicomponent Electrolyte Solutions; Pergamon Press: London, 1968. (9) Zhong, S. Y.; Sang, S. H.; Zhang, J. J.; Huang, W. Y.; Hu, J. X. Mean Activity Coefficients of NaBr in NaBr + Na2B4O7 + H2O Ternary System at 298.15 K Determined by Potential Difference Measurements. J. Chem. Eng. Data 2014, 59, 1603−1608. (10) Zhang, J. J.; Sang, S. H. Studies on mean activity coefficients of NaBr in NaBr + Na2B4O7 + H2O system at 298.15 K by potential difference measurements. J. Sichuan Univ. (Eng. Sci. Ed.) 2012, 44, 240− 243. (11) Zhou, M. F.; Sang, S. H.; Zhang, J. J.; Hu, J. X.; Zhong, S. Y. Studies on Mean Activity Coefficients of NaBr in NaBr−SrBr2−H2O Ternary System at 298.15 K by EMF Method. J. Chem. Eng. Data 2014, 59, 3779−3784. (12) Zhong, S. Y.; Sang, S. H.; Zhang, J. J.; Wei, C. Mean activity coefficients of KBr in KBr + K2B4O7 + H2O ternary system at 298.15 K determined by the electromotive force method. J. Chem. Eng. Data 2014, 59, 455−460. (13) Zhang, J. J.; Sang, S. H.; Zhong, S. Y. Mean activity coefficients of KBr in the KBr+K2SO4+H2O ternary system at 298.15 K by a potential difference method. J. Chem. Eng. Data 2012, 57, 2677−2680. (14) Clegg, S. L.; Rard, J. A.; Miller, D. G. Isopiestic Determination of the Osmotic and Activity Coefficients of NaCl + SrCl2 + H2O at 298.15 K and Representation with an Extended Ion-Interaction Model. J. Chem. Eng. Data 2005, 50, 1162−1170. (15) Zhou, M. F.; Sang, S. H.; Liu, Q. Z.; Wang, D.; Fu, Ch. Mean Activity Coefficients of NaCl in NaCl + SrCl2 + H2O Ternary. J. Chem. Eng. Data 2015, 60, 3209−3214. (16) Lindenbaum, S.; Boyd, G. E. Osmotic and activity coefficients for the symmetrical tetralkyl ammonium halides in aqueous solutions at 25 °C. J. Phys. Chem. 1964, 68, 911−917. (17) Amado, E.; Blanco, L. H. Isopiestic Determination of the Osmotic and Activity Coefficients of Aqueous Solutions of Symmetrical and Unsymmetrical Quaternary Ammonium Bromides at T = (283.15 and 288.15) K. J. Chem. Eng. Data 2009, 54, 2696−2700. (18) Amado, E. G.; Blanco, L. H. Osmotic and activity coefficients of dilute aqueous solutions of symmetrical and unsymmetrical quaternary ammonium bromides at 293.15 K. Fluid Phase Equilib. 2006, 243, 166− 170. (19) Amado, E.; Blanco, L. H. Isopiestic determination of the osmotic and activity coefficients of dilute aqueous solutions of symmetrical and unsymmetrical quaternary ammonium bromides with a new isopiestic cell at 298.15 K. Fluid Phase Equilib. 2005, 233, 230−233. (20) Blanco, L. H.; Amado, E.; Avellaneda, J. Isopiestic determination of the osmotic and activity coefficients of dilute aqueous solutions of the series MeEt3NI to HepEt3NI at 298.15 K. Fluid Phase Equilib. 2006, 249, 147−152. 2390

DOI: 10.1021/acs.jced.7b00278 J. Chem. Eng. Data 2017, 62, 2384−2391

Journal of Chemical & Engineering Data

Article

(37) Archer, D. thermodynamic properties of the NaCl + H2O System ll. Thermodynamic properties of NaCl (aq), NaCl.2 H2O(cr) and Phase equilibria. J. Phys. Chem. Ref. Data 1992, 21, 793−829. (38) Pitzer, K. S.; Simonson, J. M. Thermodynamics of Multicomponent, Miscible, Ionic Systems: Theory and Equations. J. Phys. Chem. 1986, 90, 3005−3009. (39) Zhuo, K.; Chen, Y.; Kang, L. Dielectric Constants for Binary Amino Acid-Water Solutions from (278.15 to 313.15) K. J. Chem. Eng. Data 2009, 54, 137−141. (40) Tang, J.; Ma, Y.; Li, S.; Zhai, Q.; Jiang, Y.; Hu, M. Activity Coefficients of RbCl in Ethylene Glycol + Water and Glycerol + Water Mixed Solvents at 298.15 K. J. Chem. Eng. Data 2011, 56, 2356−2361. (41) Oster, G. The dielectric properties of mixture liquids. J. Am. Chem. Soc. 1946, 68, 2036−2040. (42) Wang, P.; Anderko, A. Computation of dielectric constants of solvent mixtures and electrolyte solutions. Fluid Phase Equilib. 2001, 186, 103−122. (43) Head-Gordon, T. Is water structure around hydrophobic groups clathrate-like? Proc. Natl. Acad. Sci. U. S. A. 1995, 92, 8308−8312. (44) Schröder, C.; Neumayr, G.; Steinhauser, O. On the collective network of ionic liquid/water mixtures. III. Structural analysis of ionic liquids on the basis of Voronoi decomposition. J. Chem. Phys. 2009, 130, 194503−194511. (45) Zhuo, K. L.; Liu, H. X.; Zhang, H. H.; Liu, Y. H.; Wang, J. Activity Coefficients of CaCl in (Maltose + water) and (lactose + water) Mixtures at 298.15 K. J. Chem. Thermodyn. 2008, 40, 889−896. (46) Zhuo, K.; Ren, H.; Wei, Y.; Chen, Y.; Ma, J. Activity Coefficients of [Cmim]Br (n = 3 to 8) Ionic Liquids in Aqueous Fructose Solution at T = 298.15 K. J. Chem. Eng. Data 2014, 59, 640−648. (47) Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. G. A Smooth Particle Mesh Ewald Method. J. Chem. Phys. 1995, 103, 8577. (48) Bussi, G.; Donadio, D.; Parrinello, M. Canonical Sampling through Velocity Rescaling. J. Chem. Phys. 2007, 126, 014101. (49) Sprenger, K. G.; Jaeger, V. W.; Pfaendtner, J. The General AMBER Force Field (GAFF) Can Accurately Predict Thermodynamic and Transport Properties of Many Ionic Liquids. J. Phys. Chem. B 2015, 119, 5882−5895. (50) Parrinello, M. Polymorphic Transitions in Single Crystals: A New Molecular Dynamics Method. J. Appl. Phys. 1981, 52, 7182. (51) Hess, B.; Bekker, H.; Berendsen, H. J. C.; Fraaije, J. G. E. M. LINCS: A Linear Constraint Solver for Molecular Simulations. J. Comput. Chem. 1997, 18, 1463−1472.

2391

DOI: 10.1021/acs.jced.7b00278 J. Chem. Eng. Data 2017, 62, 2384−2391