N2O Analogy Using Molecular Simulation

For water, the three-site rigid model SPC/E,(42) the four-site rigid model TIP4P ...... and Critical Properties as Predicted by a Simple Molecular Mod...
2 downloads 0 Views 634KB Size
Article pubs.acs.org/IECR

Validation of the CO2/N2O Analogy Using Molecular Simulation Qu Chen,† Sayee Prasaad Balaji,† Mahinder Ramdin,† Juan José Gutiérrez-Sevillano,† André Bardow,‡ Earl Goetheer,§ and Thijs J. H. Vlugt*,† †

Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands ‡ Lehrstuhl für Technische Thermodynamik, RWTH Aachen University, Schinkelstrasse 8, 52062 Aachen, Germany § TNO, Leeghwaterstraat 46, 2628 CA Delft, The Netherlands S Supporting Information *

ABSTRACT: CO2 readily reacts in aqueous amine solutions. The properties of free CO2 in amine solutions are therefore difficult to obtain directly and are often predicted from the nonreacting molecule N2O due to the similarities in mass and structure. This often-used empirical “CO2/N2O analogy” is verified in this work using molecular simulation. Continuous fractional component Monte Carlo (CFCMC) simulations in the osmotic ensemble were used to compute the Henry coefficients of CO2 and N2O in the solvents water, ethanol, n-heptane, and a 30% aqueous MEA solution at a temperature of 303 K. Molecular dynamics (MD) simulations were performed to compute the self-diffusivities of CO2 and N2O in the aforementioned solvents at 303 K. Different force fields for CO2 and water were used. The computed Henry coefficients and self-diffusivities of CO2 and N2O in the solvents are in good agreement with available experimental data. The simulation results indicate that the CO2/N2O analogy is valid for aqueous MEA solution at 303 K. The Henry coefficient and self-diffusivity ratios of CO2 to N2O in water and 30% MEA solution are approximately 0.77 and 1.1, respectively. Additional simulations where all the amines have reacted with CO2 confirm that reactions have little impact on the physical absorption and diffusion properties of CO2.

1. INTRODUCTION As a byproduct of burning fossil fuels, carbon dioxide (CO2) accounts for a large portion of the greenhouse gas emissions.1,2 CO2 capture and storage (CCS) could play a significant role in reducing CO2 emissions from power plants.3−5 The removal of CO2 from natural gas or in the postcombustion process is typically achieved in an absorber using an amine solution.6−8 For example, an aqueous solution containing 30% monoethanolamine (MEA) is a frequently used solvent for CO2 capture.6 For a proper design of a gas absorption unit, it is very important to understand the mass-transfer process of CO2. This requires knowledge of the solubility and diffusivity of CO2 in aqueous amine solutions. Due to reactions between CO2 and amines, the solubility of CO2 measured in experiments is the total solubility, which combines both physical and chemical absorption. However, the physical solubility is needed to obtain the rate of the CO2 absorption.9 When the CO2 molecules absorb into an amine solvent, diffusion and reaction will take place simultaneously. As there is almost no unreacted CO2, it is practically impossible to obtain the solubility and diffusivity of CO2 in amine solutions directly from experiments.10,11 In practice, this problem is usually overcome by substituting CO2 with a similar but inert gas. Nitrous oxide (N2O) can be considered as an ideal candidate since the gases have similar molecular weights, volumes, and structures. Substituting CO2 for N2O to obtain physical properties of CO2 is known as “the CO2/N2O analogy”. According to the CO2/N2O analogy, the solubility and diffusivity of both gases are related by12 © XXXX American Chemical Society

⎡ HCO ⎤ ⎡ HCO ⎤ 2 2 ⎢ ⎥ ⎥ =⎢ ⎢⎣ HN2O ⎥⎦ ⎢ HN2O ⎥⎦ ⎣ in water in amine solution

(1)

⎡ DCO ⎤ ⎡ DCO ⎤ 2 2 ⎢ ⎥ ⎥ =⎢ ⎢⎣ DN2O ⎥⎦ ⎢ DN2O ⎥⎦ ⎣ in water in amine solution

(2)

where H is the Henry coefficient and D is the self-diffusivity. The idea of the CO2/N2O analogy was first conceived by Clarke,13 who studied the absorption kinetics of CO2 in aqueous MEA solutions at short contact times using the laminar liquid jet technique. In this study, it was suggested that the solubility and diffusivity of CO2 in MEA solutions are comparable to those of N2O due to the similarity in their mass and molecular interactions. Joosten and Danckwerts14 presented experimental data on the solubility and diffusivity of N2O in equimolar potassium−potassium bicarbonate solutions of various ionic strengths. These authors indicated that the corresponding properties of CO2 can be reliably predicted by these N2O data. Joosten and Danckwerts showed that the diffusivities of CO2 and N2O in a pure liquid and an electrolyte solution are correlated by ⎡ Dsoln ⎤ ⎡D ⎤ = ⎢ soln ⎥ ⎢ ⎥ ⎣ D0 ⎦CO ⎣ D0 ⎦N O 2

2

(3)

Received: September 3, 2014 Revised: October 23, 2014 Accepted: October 28, 2014

A

dx.doi.org/10.1021/ie503488n | Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Industrial & Engineering Chemistry Research

Article

molecules fluctuates.25 This can be realized by adding/ removing molecules to/from a so-called “open system”. A system where the number of molecules fluctuates is called an open system. The conventional Monte Carlo schemes are inefficient since they suffer from the low insertion/deletion acceptance probabilities at high densities.26,27 Therefore, methods that can handle open systems at high densities are ́ et required to efficiently compute the Henry coefficient. Lisal al.28 obtained the Henry coefficient for CO2 in water via the Widom test-particle insertion method and the staged freeenergy perturbation method. These authors also compared different interaction models of water and CO2. Cichowski et al.29 determined the Henry coefficient of CO2 in ethanol by transition matrix Monte Carlo simulations coupled with the expanded ensemble technique. A similar study on the CO2 solubility in ethanol was carried out by Zhang and Siepmann,30 who used Gibbs ensemble Monte Carlo simulations to directly calculate the Henry coefficient. The continuous fractional component Monte Carlo (CFCMC) method was recently developed by Shi and Maginn.26,31−33 This method overcomes the difficulties encountered in conventional Monte Carlo (MC) methods by utilizing gradual insertions/deletions of molecules using a continuous coupling parameter and an adaptive bias potential. The CFCMC method in the osmotic ensemble was applied to obtain the Henry coefficient of CO2 in ethanol, and the results are in agreement with previous studies.31 Self-diffusivities of the absorbed CO2 and N2O can be computed by equilibrium molecular dynamics (MD) simulations.34 This approach requires interaction potential models capable of accurately describing the intermolecular and intramolecular interactions of the molecules. Some efforts were made to investigate the influence of different force fields on the diffusion coefficient. In het Panhuis et al.35 used the SPC/E model for water and adjusted the Lennard-Jones size parameters and quadruple moments (partial charges) of the Gromos CO2 potential. This work concluded that, whereas the Lennard-Jones interaction of the CO2 potential has a large impact on the CO2 diffusion behavior, the diffusivity of CO2 is almost independent of the Coulombic interaction of the CO2 potential. Vlcek et al.36 proposed and applied a coupling parameter scheme to optimize the Lennard-Jones cross interactions for the SPC/E−EPM2 models. This successfully improved the accuracy of the predicted mutual diffusion coefficients of CO2 and water in simulations. Zeebe37 presented studies on the diffusivity of dissolved CO2 in water over the temperature range from 273 to 373 K and provided a general power-law equation to predict the diffusivity of CO2 as a function of temperature. Recently, Moultos et al.38 compared different combinations of CO2 and water models to investigate the diffusivity of CO2 over a broad range of temperatures and pressures. These authors found that at normal pressure the TIP4P/2005−EPM2 combination was best to calculate the diffusivity of CO2 at lower temperatures ( 30% MEA aqueous solution > water. To validate the CO2/N2O analogy, the ratios of the Henry coefficients of CO2 and N2O in all investigated solvents were obtained. As listed in Table 2, the CO2/N2O ratios for water Table 2. Ratios of Henry Coefficients of CO2 to N2O in Water, Ethanol, n-Heptane, and a 30% MEA Solution at a Temperature of 303 Ka

Figure 3. Fugacity of CO2 or N2O in the gas phase plotted as a function of mole fraction of CO2/N2O in ethanol and n-heptane at 303 K. The EPM2 and TraPPE models were used for CO2.

solvent water (TIP4P) water (TIP5P) ethanol n-heptane 30% MEA (TIP4P) 30% MEA (TIP5P)

and N2O in the gas phase are plotted as a function of the solute mole fraction in solution. By extrapolating to low pressures, the Henry coefficient can be obtained. As shown in Table 1, all the simulation results for the Henry coefficients of CO2 and N2O in ethanol and n-heptane underestimate the experimental values. To the best of our knowledge, experimental data for the Henry coefficients of CO2 and N2O in ethanol have been reported only at 293 K as 16.28 and 13.72 MPa, respectively.19 The two models of CO2 (EPM2 and TraPPE) in ethanol and heptane, unlike in water, result in Henry coefficients that are very similar. The Henry coefficient of N2O is always lower than that of CO2 in both ethanol and heptane. Note that the solubility is higher when the Henry coefficient is lower. This trend for ethanol was also observed in the work by Kierzkowska-Pawlak et al.,19 although their experiments were conducted at 293 K. Finally, we performed simulations for the CO2 and N2O solubilities in a 30% (w/w) MEA solution. Note that at this stage the reactions between CO2 and MEA have been “switched off” in the MC simulations. In section 3.2 a case will be presented where the reactions between CO2 and MEA have been taken into consideration. On the basis of our results of CO2 solubility in water, we only used the TIP4P and TIP5P models for water and the TraPPE model for CO2. The results are shown in Figure 4 and Table 1. Although a certain deviation from the experiment can be observed, the simulation data show that 30%

[HCO2/HN2O]sim

[HCO2/HN2O]exp

diff (%)

± ± ± ± ± ±

0.7712 0.7712 n/a 1.4760,61 n/a n/a

−7.8 +1.3 n/a −29 n/a n/a

0.71 0.78 1.11 1.05 0.82 0.84

0.02 0.02 0.05 0.02 0.02 0.04

a

Simulation results are compared with the experimental data reported with the difference presented. We only used the TraPPE model for CO2 and the TIP4P/TIP5P model for H2O for comparison.

and MEA solution are all below unity while the CO2/N2O ratios for ethanol and heptane are both above unity. The differences between simulations and experiments for the CO2/ N2O Henry coefficient are below 10% for the TIP4P and TIP5P water models, but for heptane the difference is nearly 30%. The explanation could be that the experimental ratio for heptane is in fact calculated from data of two different studies.60,61 For solvents with larger dipole moments (e.g., water), more CO2 molecules are likely to be absorbed. Solvents with smaller dipole moments (e.g., heptane) tend to absorb more N2O molecules. The validity of eq 1 can be verified by comparing the CO2/N2O ratios for the Henry coefficients in water and in 30% MEA solution. Li and Lai12 reported the ratio HCO2/HN2O in water to be 0.77 at a temperature of 303 K and E

dx.doi.org/10.1021/ie503488n | Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Industrial & Engineering Chemistry Research

Article

the other solvents was also investigated, and again no notable differences could be observed for the different system sizes (results not shown). As for the solubility, a comparison of the performance for the different models for H2O with respect to the diffusivity is provided. Three water models (i.e., SPC/E, TIP4P, and TIP5P) were tested using the EPM2 CO2 model. From the diffusivity regime up to 20 ns, the diffusivity of CO2 in the SPC/E, TIP4P, and TIP5P water was obtained as (2.727 ± 0.071) × 10−9, (2.428 ± 0.080) × 10−9, and (3.500 ± 0.079) × 10−9 m2/s, respectively, compared to the experimental value of 2.120 × 10−9 m2/s reported by Li and Lai.12 It is obvious that the combination EPM2−TIP5P performs poorly, and therefore we used the SPC/E and TIP4P models for further research. Simulations were then performed in an ensemble containing 1 CO2 solvated in 500 water molecules, 500 ethanol molecules, 500 heptane molecules, and a 30% MEA solution (56 MEA plus 444 H2O), respectively, at 303 K. Similar simulations were performed with 1 N2O molecule in the aforementioned solvents as well. Two models for H2O (SPC/E and TIP4P) and two models for CO2 (EPM2 and TraPPE) lead to four combinations. Table 3 shows the results for all simulations at

applied this ratio for calculating the Henry coefficient of CO2 in amine solutions. We found that this ratio is in agreement with our simulation results, since all the data points for water and MEA solution fall within 10% of 0.77. The CO2/N2O analogy for the Henry coefficient is therefore supported at a temperature of 303 K. So far, all the MC simulations assume that CO2 does not react with the MEA molecules; hence the solubility of CO2 in the MEA solution was purely due to physisorption. In reality, CO2 will react with MEA and it is therefore very important to understand how reactions or reaction products affect the Henry coefficients of CO2 and N2O. Assuming that all MEA molecules have been reacted with CO2, a simulation of a fully reacted s y s t e m w i t h 2 8 H O C H 2 C H 2 N H CO O − i o n s , 2 8 HOCH2CH2NH2+ ions, and 444 water molecules was performed. In experiments, CO2 can fully react with MEA molecules even at low temperatures (e.g., 293 K).62 At 303 K, the MC simulations were performed to compute the Henry coefficient of CO2. The best performing combination for CO2− H2O (TraPPE−TIP5P) was used. The Henry coefficient of free CO2 in this fully reacted MEA solution is calculated to be 107.2 ± 5.6 MPa, while the Henry coefficient in the nonreacted MEA solution is 92.3 ± 1.9 MPa. This indicates that the Henry coefficient of the MEA solution increases when the reaction products are formed. In other words, the Henry coefficient of CO2 is higher in an unreacted MEA solution than in a solution containing their reaction products. However, there is no huge impact of reaction products on the physical solubility of CO2. As a result, the CO2/N2O analogy for the physical solubility may also hold during the whole chemical reaction process. 3.2. Diffusivity. The self-diffusivity was calculated from Einstein’s equation in three dimensions: 1 1 Di = lim 6Ni m →∞ mΔt

Table 3. Self-Diffusivities of CO2 in Water, Ethanol, nHeptane, and 30% MEA Solution at 303 Ka solute CO2 (EPM2) CO2 (TraPPE) N2O CO2 (EPM2) CO2 (TraPPE) N2O CO2 (EPM2) CO2 (TraPPE) N2O CO2 (EPM2) CO2 (TraPPE) N2O CO2 (EPM2)

Ni

∑ ⟨[rl ,i(t + mΔt ) − rl ,i(t )]2 ⟩ l=1

(6)

where Ni is the total number of molecules of component i, m is the number of time steps, and Δt is the time step used in MD simulations. rl,i(t) is the position of the lth molecule of type i at time t.63 Moultos et al.38 showed that the error bar of the CO2 diffusivity in water is considerably large (>10%). This accuracy is not suitable for the purpose of this study, because the diffusivity ratio of CO2 and N2O will then have an even larger error bar and it will not be possible to obtain a value statistically different from unity. In this work, we carried out at least 50 parallel simulations, from which the statistical uncertainty was determined. An order-n algorithm64 was used to compute the mean squared displacement (MSD). The average MSDs of CO2 as a function of time for each simulation were plotted. In a log−log plot, the different distinct regimes can be observed. Before molecules collide with each other, they undergo ballistic diffusion, which occurs in a very short time ( ethanol > 30% MEA > water. Molecular dynamics simulations also correctly predict the observed CO2/ N2O diffusivities in the solvents, which obey the following order: heptane > ethanol > water > 30% MEA. The simulation results indicate that the CO2/N2O analogy is valid for water and 30% MEA solution at 303 K. The ratio of the Henry coefficients CO2/N2O is approximately 0.77 (±10%), and the ratio of the self-diffusivities CO2/N2O is around 1.1 (±10%). Additional simulations were performed to investigate the effect of a chemical reaction between CO2 and MEA on the physical solubility and self-diffusivity of CO2. The simulation results indicate that both the physical solubility and the self-diffusivity of CO2 are lower in the system containing reaction products of CO2 and MEA, which is a result of the reaction between CO2 and MEA.



ASSOCIATED CONTENT

S Supporting Information *

Force field parameters of all components; comparison of densities for 30% MEA solution at different CO2 loadings. This material is available free of charge via the Internet at http:// pubs.acs.org.



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was performed as a part of the CATO-2 program, the Dutch national R&D program on CO2 capture, transport and storage, funded by the Dutch Ministry of Economic Affairs (Q.C. and S.P.B.). Financial support by ADEM, a green Deal in Energy Materials, a program of the Dutch Ministry of Economic Affairs, is also acknowledged (M.R.). This work was also supported by NWO Exacte Wetenschappen (Physical Sciences) for the use of supercomputer facilities, with financial support from the Nederlandse Organisatie voor Wetenschappelijk Onderzoek (Netherlands Organization for Scientific Research, NWO).



REFERENCES

(1) Meinshausen, M.; Meinshausen, N.; Hare, W.; Raper, S. C. B.; Frieler, K.; Knutti, R.; Frame, D. J.; Allen, M. R. Greenhouse-gas Emission Targets for Limiting Global Warming to 2°C. Nature 2009, 458, 1156−1162. (2) Solomon, S. Irreversible Climate Change due to Carbon Dioxide Emissions. Proc. Natl. Acad. Sci. U. S. A. 2009, 106, 1704−1709. (3) Metz, B.; Davidson, O.; de Coninck, H.; Loos, M.; Meyer, L. Carbon Dioxide Capture and Storage: Special Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, U.K., 2005. H

dx.doi.org/10.1021/ie503488n | Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Industrial & Engineering Chemistry Research

Article

(27) Dubbeldam, D.; Torres-Knoop, A.; Walton, K. S. On the Inner Workings of Monte Carlo Codes. Mol. Simul. 2013, 39, 1253−1292. (28) Lísal, M.; Smith, W. R.; Aim, K. Analysis of Henry’s Constant for Carbon Dioxide in Water via Monte Carlo Simulation. Fluid Phase Equilib. 2004, 226, 161−172. (29) Cichowski, E. C.; Schmidt, T. R.; Errington, J. R. Determination of Henry’s law Constants through Transition Matrix Monte Carlo Simulation. Fluid Phase Equilib. 2005, 236, 58−65. (30) Zhang, L.; Siepmann, J. I. Direct Calculation of Henry’s Law Constants from Gibbs Ensemble Monte Carlo Simulations: Nitrogen, Oxygen, Carbon Dioxide and Methane in Ethanol. Theor. Chem. Acc. 2006, 115, 391−397. (31) Shi, W.; Maginn, E. J. Continuous Fractional Component Monte Carlo: An Adaptive Biasing Method for Open System Atomistic Simulations. J. Chem. Theory Comput. 2007, 3, 1451−1463. (32) Shi, W.; Maginn, E. J. Improvement in Molecule Exchange Efficiency in Gibbs Ensemble Monte Carlo: Development and Implementation of the Continuous Fractional Component Move. J. Comput. Chem. 2008, 29, 2520−2530. (33) Liu, H.; Maginn, E. A Molecular Dynamics Investigation of the Structural and Dynamic Properties of the Ionic Liquid 1-n-Butyl-3methylimidazolium Bis(trifluoromethanesulfonyl)imide. J. Chem. Phys. 2011, 135, 124507. (34) Liu, X.; Schnell, S. K.; Simon, J. M.; Krüger, P.; Bedeaux, D.; Kjelstrup, S.; Bardow, A.; Vlugt, T. J. H. Diffusion Coefficients from Molecular Dynamics Simulations in Binary and Ternary Mixtures. Int. J. Thermophys. 2013, 34, 1169−1196. (35) In het Panhuis, M.; Patterson, C. H.; Lynden-Bell, R. M. A Molecular Dynamics Study of Carbon Dioxide in Water: Diffusion, Structure and Thermodynamics. Mol. Phys. 1998, 94, 963−972. (36) Vlcek, L.; Chialvo, A. A.; Cole, D. R. Optimized Unlike-pair Interactions for Water-Carbon Dioxide Mixtures Described by the SPC/E and EPM2 Models. J. Phys. Chem. B 2011, 115, 8775−8784. (37) Zeebe, R. E. On the Molecular Diffusion Coefficients of Dissolved CO2, HCO3−, and CO32− and their Dependence on Isotopic Mass. Geochim. Cosmochim. Acta 2011, 75, 2483−2498. (38) Moultos, O. A.; Tsimpanogiannis, I. N.; Panagiotopoulos, A. Z.; Economou, I. G. Atomistic Molecular Dynamics Simulations of CO2 Diffusivity in H2O for a Wide Range of Temperatures and Pressures. J. Phys. Chem. B 2014, 118, 5532−5541. (39) Dubbeldam, D.; Calero, S.; Ellis, D. E.; Snurr, R. Q. RASPA 1.0: Molecular Software Package for Adsorption and Diffusion in Flexible Nanoporous Materials; 2013. (40) Peng, D.; Robinson, D. B. A New Two-Constant Equation of State. Ind. Eng. Chem. Fundam. 1976, 15, 59−64. (41) Wang, F.; Landau, D. P. Efficient, Multiple-Range Random Walk Algorithm to Calculate the Density of States. Phys. Rev. Lett. 2001, 86, 2050−2053. (42) Berendsen, H. J. C.; Grigera, J. R.; Straatsma, T. P. The Missing Term in Effective Pair Potentials. J. Phys. Chem. 1987, 91, 6269−6271. (43) Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79, 926−935. (44) Mahoney, M. W.; Jorgensen, W. L. A Five-Site Model for Liquid Water and the Reproduction of the Density Anomaly by Rigid, Nonpolarizable Potential Functions. J. Chem. Phys. 2000, 112, 8910− 8922. (45) Chen, B.; Potoff, J. J.; Siepmann, J. I. Monte Carlo Calculations for Alcohols and their Mixtures with Alkanes. Transferable Potentials for Phase Equilibria. 5. United-Atom Description of Primary, Secondary, and Tertiary Alcohols. J. Phys. Chem. B 2001, 105, 3093−3104. (46) Martin, M. G.; Siepmann, J. I. Transferable Potentials for Phase Equilibria. 1. United-Atom Description of n-Alkanes. J. Phys. Chem. B 1998, 102, 2569−2577. (47) Jorgensen, W. L.; Maxwell, D. S.; Tirado-Rives, J. Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids. J. Am. Chem. Soc. 1996, 118, 11225−11236.

(48) Potoff, J. J.; Siepmann, J. I. Vapor-Liquid Equilibria of Mixtures Containing Alkanes, Carbon Dioxide, and Nitrogen. AIChE J. 2001, 47, 1676−1682. (49) Harris, J. G.; Yung, K. H. Carbon Dioxide’s Liquid-vapor Coexistence Curve and Critical Properties as Predicted by a Simple Molecular Model. J. Phys. Chem. 1995, 99, 12021−12024. (50) Lachet, V.; Creton, B.; de Bruin, T.; Bourasseau, E.; Desbiens, N.; Wilhelmsen, O.; Hammer, M. Equilibrium and Transport Properties of CO2 + N2O and CO2 + NO Mixtures: A Molecular Simulation and Equation of State Modeling Study. Fluid Phase Equilib. 2012, 322, 66−78. (51) Hirschfelder, J. O.; Curtiss, C. F.; Brid, R. B. Molecular Theory of Gases and Liquids; John Wiley and Sons: New York, 1954. (52) Darden, T.; York, D.; Pedersen, L. Particle mesh Ewald: An N · log(N) Method for Ewald Sums in Large Systems. J. Chem. Phys. 1993, 98, 10089−10092. (53) Apol, E.; Apostolov, R.; Berendsen, H. J. C.; van Buuren, A.; Bjelkmar, P.; van Drunen, R.; Feenstra, A.; Groenhof, G.; Kasson, P.; Larsson, P.; Meulenhoff, P.; Murtola, T.; Pall, S.; Pronk, S.; Schulz, R.; Shirts, M.; Sijbers, A.; Tieleman, P.; Hess, B.; van der Spoel, D.; Lindahl, E. GROMACS User Manual, version 4.5.4; 2010. (54) Hanwell, M. D.; Curtis, D. E.; Lonie, D. C.; van der Meersch, T.; Zurek, E.; Hutchison, G. R. Avogadro: an Advanced Semantic Chemical Editor, Visualization, and Analysis Platform. J. Cheminf. 2012, 4, 17. (55) Vaidya, P. D.; Kenig, E. Y. CO2-alkanolamine Reaction Kinetics: A Review of Recent Studies. Chem. Eng. Technol. 2007, 30, 1467− 1474. (56) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark; M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, J. M.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, Ö ; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J. Gaussian 09, revision D.01; Gaussian, Inc.: Wallingford, CT, 2009. (57) Carroll, J. J.; Slupsky, J. D.; Mather, A. E. The Solubility of Carbon Dioxide in Water at Low Pressure. J. Phys. Chem. Ref. Data 1991, 20, 1201−1209. (58) Gwinner, B.; Roizard, D.; Lapicque, F.; Favre, E.; Cadours, R.; Boucot, P.; Carrette, P. L. CO2 Capture in Flue Gas: Semiempirical Approach to Select a Potential Physical Solvent. Ind. Eng. Chem. Res. 2006, 45, 5044−5049. (59) Wagner, W.; Pruss, A. The IAPWS Formulation 1995 for the Thermodynamic Properties of Ordinary Water Substance for General and Scientific Use. J. Phys. Chem. Ref. Data 2002, 31, 387−535. (60) King, M. B.; Al-Najjar, H. The Solubilities of Carbon Dioxide, Hydrogen Sulphide and Propane in some normal Alkane SolventsI: Experimental Determinations in the Range 15−70°C and Comparison with Ideal Solution Values. Chem. Eng. Sci. 1977, 32, 1241−1246. (61) Yen, L. C.; McKetta, J. J., Jr. Solubility of Nitrous Oxide in some Nonpolar Solvents. J. Chem. Eng. Data 1962, 7, 288−289. (62) Bö ttinger, W.; Maiwald, M.; Hasse, H. Online NMR Spectroscopic Study of Species Distribution in MEA-H2O-CO2 and DEA-H2O-CO2. Fluid Phase Equilib. 2008, 263, 131−143. (63) Liu, X.; Martín-Calvo, A.; McGarrity, E.; Schnell, S. K.; Calero, S.; Simon, J.-M.; Bedeaux, D.; Kjelstrup, S.; Bardow, A.; Vlugt, T. J. H. Fick Diffusion Coefficients in Ternary Liquid Systems from Equilibrium Molecular Dynamics Simulations. Ind. Eng. Chem. Res. 2012, 51, 10247−10258. I

dx.doi.org/10.1021/ie503488n | Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX

Industrial & Engineering Chemistry Research

Article

(64) Dubbeldam, D.; Ford, D. C.; Ellis, D. E.; Snurr, R. Q. A New Perspective on the order-n Algorithm for Computing Correlation Functions. Mol. Simul. 2009, 35, 1084−1097. (65) Snijder, E. D.; te Riele, M. J. M.; Versteeg, G. F.; van Swaaij, W. P. M. Diffusion Coefficients of CO, CO2, N2O, and N2 in Ethanol and Toluene. J. Chem. Eng. Data 1995, 40, 37−39. (66) Wilke, C. R.; Chang, P. Correlation of Diffusion Coefficients in Dilute Solutions. AIChE J. 1955, 1, 264−270. (67) Poling, B. E.; Prausnitz, J. M.; O’Connell, J. P. The Properties of Gases and Liquids, 5th ed.; McGraw-Hill: New York, 2001.

J

dx.doi.org/10.1021/ie503488n | Ind. Eng. Chem. Res. XXXX, XXX, XXX−XXX