Screening Deep Eutectic Solvents for CO2 Capture With COSMO-RS

GRAPHICAL ABSTRACT

Graphical Abstract. Corrected COSMO-RS model for screening CO₂ solubility in DES.

Introduction

With the rapid development of industry, the demand for energy is growing. Fossil fuels currently account for the majority of energy supply (Rahmanifard and Plaksina, 2019). The use of fossil fuels poses a range of environmental problems. For example, the use of fossil fuels emits a large amount of carbon dioxide (CO₂) (Hanif et al., 2019; Hosseini et al., 2019). This, in turn, leads to the serious greenhouse-gas effect. It was reported that the global CO₂ emissions reached an all-time high of 33.143 billion tons in 2018 by the International Energy Agency (IEA) (Qu, 2018). This calls for CO₂ capture.

There are several technologies for CO₂ capture, for example, physical, or chemical solvent scrubbing (Abdeen et al., 2016) and pressure swing adsorption (Riboldi and Bolland, 2016). Because of the complexity of the gas components, most technologies still suffer from high energy demand, high cost, and serious secondary pollutions. Developing new capture technologies, including new solvents and novel processes, is the key point to CO₂ capture.

Recently, ionic liquids (ILs) have been proposed as potential candidates for CO₂ capture due to their unique properties. However, the conventional ILs are expensive mainly due to the complex synthesis process. The newly emerged DESs possess analogous properties to ILs and share many common properties with ILs, whereas DESs have additional merits of low cost, biodegradability, low toxicity, easy preparation, and no purification requirement (Sarmad et al., 2017a). The most fascinating property of DESs is the structural diversity, and it can be prepared by mixing a hydrogen bond donor (HBD) with a hydrogen bond acceptor (HBA) in appropriate mole ratios (Hayyan et al., 2013).

To develop DESs for CO₂ capture, the CO₂ absorption capacity (i.e., CO₂ solubility) in DESs is one of the important properties. It is reported that the CO₂ solubility in DESs depends on pressure, temperature, and the mole ratio of HBD to HBA, and it increases with increasing pressure and decreasing temperature (Sarmad et al., 2017a). Based on the studies of Kamps et al. (2003) and Aki et al. (2004), it was found that the CO₂ solubilities in DESs were comparable with the imidazolium-based ILs. Leron and Li measured the CO₂ solubilities in choline chloride (ChCl)-glycerol (GLY) 1:2 at 313.15 K and 0–6 MPa (Leron and Li, 2013). The result indicated that the measured CO₂ solubilities (0.1103–3.0718 mol/kg) in this DES are greater than those in the 1-ethyl-3-methylimidazolium-based ILs, such as 1-ethyl-3-methylimidazolium tetrafluoroborate, 1-ethyl-3-methylimidazolium dicyanamide, and 1-ethyl-3-methylimidazolium ethyl sulfate. In the work by Sarmad et al., 35 DESs with 209 data points for CO₂ solubility at 298.15 K and below 2 MPa were reported (Sarmad et al., 2017b). Among the synthesized DESs, 15 samples exhibit higher CO₂ solubilities than the conventional ILs, such as 1-octyl-3-methylimidazolium tetrafluoroborate, 1-octyl-3-methylimidazolium hexafluorophosphate, and 1-ethyl-3-methylimidazolium ethylsulfate. Sarmad et al. reviewed the CO₂ solubilities in 45 DESs in wide ranges of mole ratio (HBA:HBD), temperature, and pressure (Sarmad et al., 2017a). Based on this review, top three DESs with respect to CO₂ solubilities were acquired, i.e., 3.6929, 3.5592, and 3.1265 mol/kg for ChCl-GLY (1:2, 303.15 K, 5.863 MPa), ChCl-Urea (1:2, 303.15 K, 5.654 MPa), and ChCl-ethylene glycerol (EG) (1:2, 303.15 K, 0.5774 MPa), respectively.

Numerous possible DESs can be synthesized. It is nevertheless a challenge to identify and suggest the best DESs for capturing CO₂ based on the available experimental measurements only. A rapid and priori screening method to predict the CO₂ absorption capacity in DESs is needed. COSMO-RS (Conductor-like Screening Model for Real Solvents) is recommended for predicting some thermodynamic properties (Gonzalez-Miquel et al., 2011; Liu Y.-R., et al., 2016), such as activity coefficients, solubilities, and Henry's constants, and it can also be used for calculating the properties of mixtures at various temperatures and pressures, rendering it an effective tool to predict and develop task-specific DESs for a specific application. Previous work has demonstrated that Henry's constants can be used as one of the criteria to screen ILs for CO₂ capture (Manan et al., 2009; Palomar et al., 2011), and the Henry's constants of CO₂ in ILs have been successfully estimated with COSMO-RS (Gonzalez-Miquel et al., 2011). However, using COSMO-RS to predict the CO₂ solubility or the Henry's constant of CO₂ in DESs is still scarce. To the best of our knowledge, only Kamgar et al. (2017b) predicted CO₂ solubility in the DES of ChCl-Urea (1:2) with COSMO-RS, demonstrating a reliable prediction only at low pressures and high temperatures, where the gas could be assumed to be an ideal gas. No report is available to predict the Henry's constants of CO₂ in DESs and compare with the experimental results for verification.

The aim of this work was to predict CO₂ solubilities, the Henry's constants of CO₂, and the interactions for CO₂ capture in DESs with COSMO-RS. A comprehensive survey of the published experimental results of CO₂ solubility and Henry's constant of CO₂ was firstly carried out. COSMO-RS was used to predict these two properties, and the predictions were compared with the experimental results. COSMO-RS was further developed with a correction based on the experimental CO₂ solubility to improve the model performance. The σ-profiles predicted by COSMO-RS were used to reflect the strength of molecular interactions between an HBA (or HBD) and CO₂, and the calculated excess enthalpy was applied to acquire the dominant interactions for CO₂ capture in DESs.

Computational Detail

COSMO-RS Computation Details

COSMO-RS calculations were performed using the software COSMOtherm (version C3.0, release 14.01, applied with parameterization BP_TZVP_C30_1401, COSMOlogic, Leverkusen, Germany) (Larriba et al., 2017). Following the standard method, first, the quantum chemical Gaussian09 package was used to optimize the structure of the studied compounds at the B3LYP/6-31++G (d, p) level. Second, the COSMOfiles of the optimized structures were opened by Gaussian03, and the COSMO continuum solution models were obtained using the BVP86/TZVP/DGA1 level theory. Third, CO₂ solubility, the Henry's constant of CO₂, σ-profiles, and excess enthalpy were determined with COSMO-RS (Liu Y.-R., et al., 2016). In computation, the temperature and pressure were set to be the same values as the experimental conditions.

In COSMO-RS, all DESs were implemented in COSMOtherm software following the electroneutral approach, where each DES was treated as three different compounds in a stoichiometric mixture (Larriba et al., 2017). For HBA, the mole fractions of the cation and anion were treated as equal, i.e., n_cation = n_anion = n_HBA. The solubility and Henry's constant of CO₂ were calculated based on the equations described by Klamt et al. (2001) and Loschen and Klamt (2014). The total and contribution of the excess enthalpies were calculated with the method reported by Casas et al. (2012).

In COSMO-RS, the solubility (X) (i.e., the mole fraction) of CO₂ in DESs was obtained with the following equation (Equation 1) (Li et al., 2014).

$\begin{array}{ll}{X}_{{CO}_2}=\frac{n_{{CO}_2}}{n_{{CO}_2}+n_{DES}}& (1)\end{array}$

where n_DES is the mole quantity of absorbent, and it was obtained according to the mass weight and the mole mass of DESs.

To evaluate the model performance, the discrepancies between the results (i.e., CO₂ solubility X, Henry's constant H) estimated with COSMO-RS and the corresponding experimental data points were quantified with the absolute relative error (ARD) as defined by Equation 2 (Kamgar et al., 2017a).

$\begin{array}{ll}\text{ARD}\%=\frac{100}{NP}{{\displaystyle \sum }}_{i=1}^{NP}\left|\frac{Y^{Exp.}-Y^{\text{COSMO}-\text{RS}}}{Y^{Exp.}}\right|& (2)\end{array}$

where NP is the total number of data points. Y^COSMO−RS (i.e., $\text{ln }{\text{X}}_{\text{C}{\text{O}}_2}^{\text{COSMO}-\text{RS}}$ or H^COSMO−RS) is the result predicated with COSMO-RS (original or corrected) at a given temperature and pressure, and Y^Exp. is the corresponding experimental result (i.e., $\text{ln }{X}_{C{O}_2}^{Exp.}$ or H^Exp.).

DESs-Database and COSMOfiles

A literature survey was conducted, and the DESs with CO₂ solubility and the Henry's constant of CO₂ in DESs were summarized and used as databases for verifying and further developing COSMO-RS. According to the survey, for CO₂ solubility, the DESs with HBA:HBD at the mole ratios of 1:2, 1:3, 1:4, and 1:5 have more experimental results compared to those with other mole ratios, and they were selected for developing COSMO-RS model in this work, while the limited experimental results at the mole ratio of 1:6 were used for verification of the developed COSMO-RS. For Henry's constants, the DESs (1:1.5, 1:2, 1:3, 1:4, 1:5) were chosen due to that sufficient experimental results are available for these DESs in literatures. The DESs databases together with the experimental measurement conditions are summarized in Tables 1, 2, respectively, and the detailed CO₂ solubilities and the Henry's constants under different conditions are provided in Tables S1–S6.

TABLE 1

Table 1. DESs database studied in this work for CO₂ solubility.

TABLE 2

Table 2. DESs database studied in this work for Henry's constants of CO₂.

To use COSMO-RS, the COSMOfiles for all the studied HBAs and HBDs are needed. In this work, those for the HBAs of ATPP⁺, AC⁺, BTEA⁺, BTMA⁺, and MTPP⁺, and for the HBDs of lactic acid (LA), ethylene glycol (EG), levulinic acid (LEV), furfuryl alcohol (FA), triethylene glycol (TEG), guaiacol (GC), decanoic acid (DecA), 1,4-butanediol (1,4-BUT), and 2,3-butanediol (2,3-BUT) were calculated based on the procedures described in the computation details. The optimized structures of HBAs and HBDs as well as the COSMOfiles from Gaussian are provided in Supplementary Information. The COSMOfiles for other HBAs and HBDs studied in this work as listed in Tables 1, 2 were directly taken from the COSMO-RS database.

Model Results and Discussion

CO₂ Solubility With COSMO-RS

For conducting COSMO-RS prediction, 502 experimental data points of CO₂ solubilities in four types of DESs (1:2, 1:3, 1:4, 1:5) under different conditions summarized from the published work were used as the input to predict the CO₂ solubility. The model predictions were then compared with the experimental results. Together with the selected 502 experimental data points of CO₂ solubilities, the predicted results and the corresponding ARDs are listed in Tables S1–S4.

Taking the results listed in Table S1 as examples, it can be observed that COSMO-RS is capable of predicting CO₂ solubilities in DESs, but qualitatively. For example, the experimental ln X_CO2 of TMACl-LA (1:2) is −5.1743 at 308 K, 194 kPa, and −5.2419 at 318 K, 194 kPa. The increase of temperature from 308 to 318 K leads to a decrease of CO₂ solubility with a difference of (Δln X_CO2 = −0.0676). The corresponding COSMO-RS predicted values of ln X_CO2 are −2.9257 and −3.1451, respectively, indicating a temperature increase results in a decreased CO₂ solubility, however, the difference is (Δln X_CO2 = −0.2194) that are much lower than the experimental observation. For different kinds of DESs, the predicted CO₂ solubility also provides the same trend with the experiment results, for example, the experimental ln X_CO2 for TBACl-LA (1:2) is −4.3513 at 308 K, 194 kPa, which is greater than that of TMACl-LA (1:2) under the same condition. The prediction of COSMO-RS provides a value of ln X_CO2 = −2.2662 for TBACl-LA (1:2) at 308 K, 194 kPa, which is also higher than that for TMACl-LA (1:2) (i.e., ln X_CO2 = −2.9257). The results for other DESs listed in Tables S2–S4 show similar observations as those listed in Table S1. Hence, COSMO-RS can be used to screen DESs qualitatively, reducing the amount of unnecessary experimental measurements.

Further comparisons were performed based on ln X_CO2 between the experimental and COSMO-RS results. As shown in Figure 1, systematic deviations can be observed, i.e., all the COSMO-RS predictions are higher than the experimental results. With increasing pressure and decreasing temperature, the discrepancies become larger, which is consistent with the observation by Kamgar et al. (2017b). This indicates that using the COSMO-RS with the parameters obtained from the conventional compound systems to predict the thermodynamic properties of DESs will lead to large deviations, which was also pointed out by others (Han et al., 2018). One of the reasons can be attributed to the formation of nanoscopic structures in DESs (Cerajewski et al., 2018), which is completely different from the solvation properties of molecular solutes in the conventional solvents at the microscopic scale.

FIGURE 1

Figure 1. Experimental ln X_CO2 vs. the ln X_CO2 predicted with the original COSMO-RS. (A) HBA:HBD 1:2; (B) HBA:HBD 1:3; (C) HBA:HBD 1:4; (D) HBA:HBD 1:5.

Henry's Constant With COSMO-RS

The Henry's constants of CO₂ in DESs (1:1.5, 1:2, 1:3, 1:4, 1:5) with 132 data points were predicated with COSMO-RS and compared with the experimental results (Figure 2). The specific values of the experimental and predicted Henry's constants are listed in Table S6, and the ARDs are summarized in Table 3. From Table S6, it can be found that the predicted Henry's constants are in agreement with the experimental data, and they increase with increasing temperature. As shown in Figure 2, most of the predicted Henry's constants are lower than the experimental values, which is consistent with the observation for the CO₂ solubility because the Henry's constant is inversely proportional to the CO₂ solubility (Liu X. Y., et al., 2016). A 13.7–36.3% deviation can be observed as listed in Table 3, indicating that the predicated Henry's constants are more accurate than the predicted CO₂ solubility with the original COSMO-RS (Table 4).

FIGURE 2

Figure 2. Experimental vs. predicted Henry's constants with the original COSMO-RS.

TABLE 3

Table 3. ARDs for the Henry's constants determined experimentally with respect to those predicted with the original COSMO-RS.

TABLE 4

Table 4. Corrected COSMO-RS for predicting ln X_CO2 in the DESs of four mole ratios.

COSMO-RS Correction for CO₂ Solubility

Although, from one side, it is unfortunate to observe a large deviation between the experimental CO₂ solubility and the COSMO-RS prediction, from the other side, it illustrates that it is possible to develop a systematic correction of COSMO-RS for improving the model performance of CO₂ solubility in DESs. Liu et al. reported a corrected COSMO-RS for predicting the activity coefficient of CO₂ in ILs and acquired a good agreement between the experimental and predicted results after correction (Liu et al., 2018). Following this idea, in this work, a temperature-pressure-dependent correction was firstly proposed as summarized in Equations 3–6 (Table 4) for the DESs with the same mole ratio. As displayed in Table 4, the corrected COSMO-RS includes three parameters, i.e., A (K⁻¹), B (KPa⁻¹), and C, and these parameters were adjusted based on the experimental CO₂ solubility at different temperatures and pressures for each group of DESs, i.e., the DESs with the same mole ratio. The predicted results with the corrected COSMO-RS are given in Tables S1–S4. The performance of the corrected COSMO-RS is further illustrated in Figure 3, and the deviations in ARDs between the experimental and modeling results are reported in Table 4. It can be found that, with the corrected COSMO-RS, i.e., the COSMO-RS with a temperature-pressure-dependent parameter, the predicted logarithmic CO₂ solubilities (square symbol in Figure 3) are in agreement with the experimental results, with much smaller ARDs of 6.5, 4.8, 6.5, and 4.5% for these four groups of DESs compared to the ARDs with the original COSMO-RS. According to the results listed in Tables S1–S4, the ARD decreases with increasing temperature, which agrees with the observation by Kamgar et al. (2017b). As shown in Table 4, the adjustable parameters A (K), B (KPa⁻¹), and C in four corrected Equations 3–6 are almost linearized with the mole ratios. To further improve the model prediction capacity, the adjustable parameters A, B, and C were set to be mole-ratio-dependent which can be described as follows:

$\begin{array}{ll}A=k_1\times w+k_2& (7)\end{array}$

$\begin{array}{ll}B=k_3\times w+k_4& (8)\end{array}$

$\begin{array}{ll}C=k_5\times w+k_6& (9)\end{array}$

where w is the mole ratio, being 2, 3, 4, 5, etc.

FIGURE 3

Figure 3. Experimental ln X_CO2 vs. the ln X_CO2 calculated with the corrected COSMO-RS models using Equations 3–6 (square symbol) and Equation 10 (circular symbol). (A) HBA:HBD 1:2; (B) HBA:HBD 1:3; (C) HBA:HBD 1:4; (D) HBA:HBD 1:5.

To obtain k₁-k₆, in this work, the experimental results of CO₂ solubilities in the DESs with four mole-ratios of 1:2, 1:3, 1:4, and 1:5 were used, and, in total, 502 experimental data points were included in fitting with a linear least square method. The fitted parameters of k₁-k₆ together with the corresponding ARDs are listed in Table 5, and the fitted CO₂ solubilities are given in Tables S1–S4.

TABLE 5

Table 5. Adjustable parameters of Equation 10 and the corresponding ARDs.

To illustrate the fitting performance with this set of universal parameters (k₁-k₆), the CO₂ solubilities obtained with the corrected COSMO-RS (circular symbol) were compared with the experimental results as well as those with the corrected COSMO-RS but using the individual parameters at each mole-ratio. As shown in Figure 4, the predicted $\text{ln }{X}_{C{O}_2}^{corr.}$ with the corrected COSMO-RS using the universal parameters (i.e., Equation 10) agrees with the experimental data, and only a slight deviation can be observed between these two corrected COSMO-RS models, i.e., with Equation 10 and Equations 3–6, respectively. Additionally, from Table 5, it can be found that the ARDs with Equation 10 are 6.8, 5.2, 6.6, and 4.7% for the DESs at mole ratios of 1:2, 1:3, 1:4, and 1:5, respectively, which are almost the same as the ARDs calculated individually with Equations 3–6 (Table 4). All these indicate that it is reasonable to use the mole-ratio-dependent parameters.

$\begin{array}{l}\text{ln }{X}_{C{O}_2}^{corr.}=\text{ln }{X}_{C{O}_2}^{COSMO-RS}+\left(k_1\times w+k_2\right)\times \frac{1}{T}+\left(k_3\times w+k_4\right)\\ \times P+(k_5\times w+k_6)& (10)\end{array}$

FIGURE 4

Figure 4. Experimental ln X_CO2 vs. the ln X_CO2 predicted with the original and corrected COSMO-RS (Equation 10) at the mole ratio of 1:6.

Further Verification of the Corrected COSMO-RS

COSMO-RS with the universal parameters k₁–k₆ can be used to predict CO₂ solubilities in the DESs with any mole ratios according to Equation 10_. In order to further investigate the model prediction capacity, it was used to predict CO₂ solubilities in the DESs at the mole ratio of 1:6. This type of DES was excluded in parameter fitting for correcting COSMO-RS due to the limited number of available experimental data points. The comparison of the prediction with the experimental data as well as those predicted with the original COSMO-RS is displayed in Table S5 and Figure 4. From Figure 4, it can be seen that the ARD between $\text{ln }{X}_{C{O}_2}^{COSMO-RS}$ and $\text{ln }{X}_{C{O}_2}^{Exp.}$is 78.2%, while it is only 10.3% between $\text{ln }{X}_{C{O}_2}^{corr.}$ and $\text{ln }{X}_{C{O}_2}^{Exp.}$, indicating that the corrected COSMO-RS is applicable for different kinds of DESs.

In summary, COSMO-RS is a predictive model, but the performance is not always satisfactory for the DESs with CO₂. The model performance of CO₂ solubility in DESs can be improved significantly with the corrected COSMO-RS.

Further Investigation With COSMO-RS and Discussion

σ-Profiles Prediction

The molecular interactions between two compounds can be linked to their σ-profiles, the wider the complementary of their σ-profiles (i.e., in the same region, the σ-profile for one compound increases, while that for the other compound decreases), the stronger the molecular interactions between these two compounds. In this work, in order to study the effects of HBAs and HBDs on the interactions with CO₂, the HBAs of TMACl, TBACl, TBABr, and BTMACl as well as the HBDs of AC, GLY, and LA were selected to predict the σ-profiles. The predicted results are listed in Table S7 and illustrated in Figure 5. Following the previous work, the σ-profiles can be divided into three regions: H-bond donor region (σ < −0.0082 e/Ų), non-polar region (−0.0082 < σ < +0.0082 e/Ų), and H-bond acceptor region (σ > +0.0082 e/Ų) (Liu Y. R. et al., 2017). In the H-bond acceptor region, the σ-profiles of HBAs for TMACl and TBACl are almost the same. However, the σ-profile of TBACl has a wider region complemented with CO₂ than that for TBACl in both H-bond donor and non-polar regions, indicating that TBACl has a strong interaction with CO₂ compared to TMACl. Meanwhile, the location of the σ-profile peak in the H-bond acceptor region also reflects the interaction strength with CO₂. TBACl and TBABr have the same σ-profile curves in the H-bond donor and non-polar regions, while in the H-bond acceptor region, the peak of TBACl (0.019 e/Ų) is located to the right of TBABr (0.017 e/Ų), evidencing that TBACl has a stronger interaction with CO₂ than TBABr. The σ-profile prediction shows that both the alkyl chain length of cations and the different types of anions of HBAs can affect the CO₂ solubility, agreeing with the experimental results of X_CO2 (TBACl-LA 1:2) > X_CO2 (TMACl-LA 1:2) (Zubeir et al., 2014) and X_CO2 (TBACl-AC 1:2) > X_CO2 (TBABr-AC 1:2) (Sarmad et al., 2017b) at the same temperature and pressure. Additionally, as shown in Figure 5, it can be seen that AC is more complementary with CO₂ than GLY, and the σ-profile curve of AC is lower than that of GLY. This implies that AC as the HDB in DES has a strong interaction with CO₂ with respect to GLY, being consistent with the experimental results of X_CO2 (BTMACl-AC 1:2) > X_CO2 (BTMACl-GLY 1:2) at the same temperature and pressure (Sarmad et al., 2017b). Therefore, the σ-profile can be used to reflect the interaction strength of a DES with CO₂.

FIGURE 5

Figure 5. σ-profiles of HBAs, HBDs, and CO₂.

Excess Enthalpy Prediction

The total excess enthalpy of six kinds of DESs at infinitely dilute CO₂ concentrations predicted by COSMO-RS was illustrated in Figure 6. These DESs can be divided into three types: (1) TMACl-LA 1:2 and TBACl-LA 1:2 have the same HBD and molar ratio but different HBA; (2) ChCl-DEG 1:3 and ChCl-TEG 1:3 have the same HBA and molar ratio but different HBD; (3) ChCl-FA 1:3 and ChCl-FA 1:5 have the same HBA and HBD but different molar ratios. The higher the absolute values of excess enthalpy, the stronger the interaction between DES and CO₂, i.e., the higher the capacity for CO₂ capture. As shown in Figure 6, the absolute excess enthalpy of TBACl-LA 1:2+CO₂ is higher than that for TMACl-LA 1:2+CO₂, indicating that TBACl-LA 1:2 has a higher capacity for CO₂ capture. The comparison of ChCl-DEG 1:3+CO₂ and ChCl-TEG 1:3+CO₂ shows that ChCl-TEG 1:3+CO₂ has a high absolute excess enthalpy compared to ChCl-DEG 1:3+CO₂ system, being in agreement with the experimental results obtained by Sarmad et al. (2017b) and Ghaedi et al. (2017), i.e., increasing the alkyl chain length in HBA and HBD results in an increased CO₂ solubility. Moreover, by increasing the molar ratio of HBD in DES, the absolute excess enthalpy of ChCl-FA 1:5+CO₂ is increased compared to ChCl-FA 1:3+CO₂ system, which agrees with the experimental results (Lu et al., 2015). Therefore, predicting excess enthalpy may be an efficient option for designing the potential DES for CO₂ capture.

FIGURE 6

Figure 6. Excess enthalpy of DESs+CO₂ system.

Figure 7 displays the contribution of each interaction [H-bond (HB), misfit (MF), and Van der Waals force (VdW)] to the total excess enthalpy of DESs+CO₂ systems. It is observable in Figure 7 that HB and VdW are consistently the dominant interactions for the DES+CO₂ systems, followed by MF. Shukla et al. reported that the CO₂ solubility in DES depends on the HB interactions between HBA and HBD (Shukla and Mikkola, 2018). Cao et al. indicated that the formation of HB between HBA and HBD enhanced the CO₂ solubility in DESs (Cao et al., 2015). In addition, Atilhan et al. investigated the interactions between DES and SO₂ by the quantum chemistry, confirming a dominant VdW interaction between DES and SO₂ (Atilhan et al., 2019). These findings support our results.

FIGURE 7

Figure 7. Contribution to the excess enthalpy of DESs+CO₂ system.

Conclusion

This work established a database containing 502 experimental data points for CO₂ solubility and 132 for the Henry's constant of CO₂ in DESs. This database was used for further verification and development of COSMO-RS.

The logarithmic CO₂ solubility predicted with the original COSMO-RS shows ARDs of 62.2, 59.6, 63.0, and 59.1% for the DESs with the HBD:HBA at 1:2, 1:3, 1:4, and 1:5, respectively. The Henry's constants of CO₂ in the DESs (1:1.5, 1:2, 1:3, 1:4, 1:5) predicted with the original COSMO-RS with ARDs of 13.7–36.3% are more accurate compared to the predicted CO₂ solubilities. To improve the performance, COSMO-RS was corrected based on 502 data points of CO₂ solubility in the DESs (1:2, 1:3, 1:4, 1:5). It shows that the adjustable parameters in the corrected COSMO-RS can be universal, the corrected COSMO-RS with the universal parameters can be used to reliably predict the CO₂ solubility in DESs, and the ARDs for the logarithmic CO₂ solubility in the DESs (1:2, 1:3, 1:4, 1:5) are of 6.8, 5.2, 6.6 and 4.7%, respectively. The corrected COSMO-RS with the universal parameters was further used to predict CO₂ solubility in the DESs (1:6), showing that a much lower ARD (10.3%) compared to that with the original COSMO-RS (ARD, 78.2%). Additional, the σ-profiles can reflect the strength of molecular interactions between an HBA (or HBD) and CO₂, and the dominant interactions for CO₂ capture in DESs are the H-bond and Van der Waals force, followed by the misfit, according to the results of excess enthalpies.

This work provides a reliable tool for DESs screening and the corrected COSMO-RS can be used to quantitatively predict CO₂ solubilities in DESs.

Data Availability Statement

All datasets generated for this study are included in the article/Supplementary Material.

Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

Funding

This work was financially supported by the National Key R&D Program of China (No. 2018YFB0605802), Carl Tryggers Stiftelse foundation, Swedish Energy Agency (P47500-1). HY and XJ thank for the financial support from the Joint Research Fund for Overseas Chinese Scholars and Scholars in Hong Kong and Macao Young Scholars (No. 21729601). YN thanks for the financial support from the National Natural Science Foundation of China (No. 21776276) and Zhengzhou High Level Talent (No. 2018030045). SZe thanks for the financial support from the National Natural Science Foundation of China (No. 21890764).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2020.00082/full#supplementary-material

References

Abdeen, F. R. H., Mel, M., Jami, M. S., Ihsan, S. I., and Ismail, A. F. (2016). A review of chemical absorption of carbon dioxide for biogas upgrading. Chinese J. Chem. Eng. 24, 693–702. doi: 10.1016/j.cjche.2016.05.006

CrossRef Full Text | Google Scholar

Aki, S. N. V. K., Mellern, B. R., Saurer, E. M., and Brennecke, J. F. (2004). High-pressure phase behavior of carbon dioxide with imidazolium-based ionic liquids. J. Phys. Chem. B 108, 20355–20365. doi: 10.1021/jp046895+

CrossRef Full Text | Google Scholar

Atilhan, M., Altamash, T., and Aparicio, S. (2019). Quantum chemistry insight into the interactions between deep eutectic solvents and SO₂. Molecules 24:2963. doi: 10.3390/molecules24162963

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, L. D., Huang, J. H., Zhang, X. P., Zhang, S. J., Gao, J. B., and Zeng, S. J. (2015). Imidazole tailored deep eutectic solvents for CO₂ capture enhanced by hydrogen bonds. Phys. Chem. Chem. Phys. 17, 27306–27316. doi: 10.1039/c5cp04050g

PubMed Abstract | CrossRef Full Text | Google Scholar

Casas, A., Palomar, J., Alonso, M. V., Oliet, M., Omar, S., and Rodriguez, F. (2012). Comparison of lignin and cellulose solubilities in ionic liquids by COSMO-RS analysis and experimental validation. Ind. Crop. Prod. 37, 155–163. doi: 10.1016/j.indcrop.2011.11.032

CrossRef Full Text | Google Scholar

Cerajewski, U., Trager, J., Henkel, S., Roos, A. H., Brehm, M., and Hinderberger, D. (2018). Nanoscopic structures and molecular interactions leading to a dystectic and two eutectic points in [EMIm][Cl]/urea mixtures. Phy. Chem. Chem. Phy. 20, 29591–29600. doi: 10.1039/C8CP04912B

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y. F., Ai, N., Li, G. H., Shan, H. F., Cui, Y. H., and Deng, D. S. (2014). Solubilities of carbon dioxide in eutectic mixtures of choline chloride and dihydric alcohols. J. Chem. Eng. Data 59, 1247–1253. doi: 10.1021/je400884v

CrossRef Full Text | Google Scholar

Ghaedi, H., Ayoub, M., Sufian, S., Shariff, A. M., Hailegiorgis, S. M., and Khan, S. N. (2017). CO₂ capture with the help of Phosphonium-based deep eutectic solvents. J. Mol. Liq. 243, 564–571. doi: 10.1016/j.molliq.2017.08.046

CrossRef Full Text | Google Scholar

Gonzalez-Miquel, M., Palomar, J., Omar, S., and Rodriguez, F. (2011). CO₂/N₂ selectivity prediction in supported ionic liquid membranes (SILMs) by COSMO-RS. Ind. Eng. Chem. Res. 50, 5739–5748. doi: 10.1021/ie102450x

CrossRef Full Text | Google Scholar

Haider, M. B., Jha, D., Sivagnanam, B. M., and Kumar, R. (2018). Thermodynamic and kinetic studies of CO₂ capture by glycol and amine-based deep eutectic solvents. J. Chem. Eng. Data 63, 2671–2680. doi: 10.1021/acs.jced.8b00015

CrossRef Full Text | Google Scholar

Han, J. L., Dai, C. N., Yu, G. Q., and Lei, Z. G. (2018). Parameterization of COSMO-RS model for ionic liquids. Green Energy Environ. 3, 247–265. doi: 10.1016/j.gee.2018.01.001

CrossRef Full Text | Google Scholar

Hanif, I., Aziz, B., and Chaudhry, I. S. (2019). Carbon emissions across the spectrum of renewable and nonrenewable energy use in developing economies of Asia. Renew. Energ. 143, 586–595. doi: 10.1016/j.renene.2019.05.032

CrossRef Full Text | Google Scholar

Hayyan, M., Hashim, M. A., Hayyan, A., Al-Saadi, M. A., AlNashef, I. M., Mirghani, M. E., et al. (2013). Are deep eutectic solvents benign or toxic? Chemosphere 90, 2193–2195. doi: 10.1016/j.chemosphere.2012.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Hosseini, E., Stevens, G. W., and Scholes, C. A. (2019). Membrane gas-solvent contactors undergoing oscillating solvent flow for enhanced carbon dioxide capture. Sep. Purif. Technol. 227:115653. doi: 10.1016/j.seppur.2019.05.095

CrossRef Full Text | Google Scholar

Kamgar, A., Bakhtyari, A., Mohsenpour, S., D'Agostino, C., Moggridge, G. D., and Rahimpour, M. R. (2017a). Mutual diffusion in concentrated liquid solutions: a new model based on cluster theory. J. Mol. Liq. 232, 516–521. doi: 10.1016/j.molliq.2017.02.079

CrossRef Full Text | Google Scholar

Kamgar, A., Mohsenpour, S., and Esmaeilzadeh, F. (2017b). Solubility prediction of CO₂, CH₄, H₂, CO and N₂ in Choline Chloride/Urea as a eutectic solvent using NRTL and COSMO-RS models. J. Mol. Liq. 247, 70–74. doi: 10.1016/j.molliq.2017.09.101

CrossRef Full Text | Google Scholar

Kamps, Á. P.-S., Tuma, D., Xia, J. Z., and Maurer, G. (2003). Solubility of CO₂ in the ionic liquid [bmim][PF6]. J. Chem. Eng. Data 48, 746–749. doi: 10.1021/je034023f

CrossRef Full Text | Google Scholar

Klamt, A., Eckert, F., and Hornig, M. (2001). COSMO-RS: a novel view to physiological solvation and partition questions. J. Comput. Aid. Mol. Des. 15, 355–365. doi: 10.1023/a:1011111506388

PubMed Abstract | CrossRef Full Text | Google Scholar

Larriba, M., Ayuso, M., Navarro, P., Delgado-Mellado, N., Gonzalez-Miquel, M., García, J., et al. (2017). Choline chloride-based deep eutectic solvents in the dearomatization of gasolines. ACS Sustain. Chem. Eng. 6, 1039–1047. doi: 10.1021/acssuschemeng.7b03362

CrossRef Full Text | Google Scholar

Leron, R. B., and Li, M.-H. (2013). Solubility of carbon dioxide in a eutectic mixture of choline chloride and glycerol at moderate pressures. J. Chem. Thermodyn. 57, 131–136. doi: 10.1016/j.jct.2012.08.025

CrossRef Full Text | Google Scholar

Li, G. H., Deng, D. S., Chen, Y. F., Shan, H. F., and Ai, N. (2014). Solubilities and thermodynamic properties of CO₂ in choline-chloride based deep eutectic solvents. J. Chem. Thermodyn. 75, 58–62. doi: 10.1016/j.jct.2014.04.012

CrossRef Full Text | Google Scholar

Li, X. Y., Hou, M. H., Han, B. X., Wang, X. L., and Zou, L. Z. (2008). Solubility of CO₂ in a choline chloride + urea eutectic mixture. J. Chem. Eng. Data 53, 548–550. doi: 10.1021/je700638u

CrossRef Full Text | Google Scholar

Liu, X. B., Gao, B., Jiang, Y. T., Ai, N., and Deng, D. S. (2017). Solubilities and thermodynamic properties of carbon dioxide in guaiacol-based deep eutectic solvents. J. Chem. Eng. Data 62, 1448–1455. doi: 10.1021/acs.jced.6b01013

CrossRef Full Text | Google Scholar

Liu, X. Y., Huang, Y., Zhao, Y. S., Gani, R., Zhang, X. P., and Zhang, S. J. (2016). Ionic liquid design and process simulation for decarbonization of shale gas. Ind. Eng. Chem. Res. 55, 5931–5944. doi: 10.1021/acs.iecr.6b00029

CrossRef Full Text | Google Scholar

Liu, X. Y., Zhou, T., Zhang, X. P., Zhang, S. J., Liang, X. D., Gani, R., et al. (2018). Application of COSMO-RS and UNIFAC for ionic liquids based gas separation. Chem. Eng. Sci. 192, 816–828. doi: 10.1016/j.ces.2018.08.002

CrossRef Full Text | Google Scholar

Liu, Y.-R., Thomsen, K., Nie, Y., Zhang, S.-J., and Meyer, A. S. (2016). Predictive screening of ionic liquids for dissolving cellulose and experimental verification. Green Chem. 18, 6246–6254. doi: 10.1039/C6GC01827K

CrossRef Full Text | Google Scholar

Liu, Y. R., Meyer, A. S., Nie, Y., Zhang, S. J., Zhao, Y. S., Fosbøl, P. L., et al. (2017). Freezing point determination of water–ionic liquid mixtures. J. Chem. Eng. Data 62, 2374–2383. doi: 10.1021/acs.jced.7b00274

CrossRef Full Text | Google Scholar

Loschen, C., and Klamt, A. (2014). Prediction of solubilities and partition coefficients in polymers using COSMO-RS. Ind. Eng. Chem. Res. 53, 11478–11487. doi: 10.1021/ie501669z

CrossRef Full Text | Google Scholar

Lu, M. Z., Han, G. Q., Jiang, Y. T., Zhang, X. D., Deng, D. S., and Ai, N. (2015). Solubilities of carbon dioxide in the eutectic mixture of levulinicacid (or furfuryl alcohol) and choline chloride. J. Chem. Thermodyn. 88, 72–77. doi: 10.1016/j.jct.2015.04.021

CrossRef Full Text | Google Scholar

Manan, N. A., Hardare, C., Jacquemin, J., Rooney, D. W., and Youngs, T. G. A. (2009). Evaluation of gas solubility prediction in ionic liquids using COSMOthermX. J. Chem. Eng. Data 54, 2005–2022. doi: 10.1021/je800857x

CrossRef Full Text | Google Scholar

Palomar, J., Gonzalez-Miquel, M., Polo, A., and Rodriguez, F. (2011). Understanding the physical absorption of CO₂ in ionic liquids using the COSMO-RS method. Ind. Eng. Chem. Res. 50, 3452–3463. doi: 10.1021/ie101572m

CrossRef Full Text | Google Scholar

Qu, T. F. (2018). The International Energy Agency Report: Global CO2 Emissions Researched an All-Time High of 33.143 Billion Tons in 2018. Available online at: http://world.huanqiu.com/exclusive/2019-03/14612504.html?agt=15438 (accessed October 1, 2019).

Rahmanifard, H., and Plaksina, T. (2019). Hybrid compressed air energy storage, wind and geothermal energy systems in Alberta: Feasibility simulation and economic assessment. Renew. Energ. 143, 453–470. doi: 10.1016/j.renene.2019.05.001

CrossRef Full Text | Google Scholar

Riboldi, L., and Bolland, O. (2016). Pressure swing adsorption for coproduction of power and ultrapure H₂ in an IGCC plant with CO₂ capture. Int. J. Hydrog. Energ. 41, 10646–10660. doi: 10.1016/j.ijhydene.2016.04.089

CrossRef Full Text | Google Scholar

Sarmad, S., Mikkola, J. P., and Ji, X. Y. (2017a). Carbon dioxide capture with ionic liquids and deep eutectic solvents: a new generation of sorbents. ChemSusChem 10, 324–352. doi: 10.1002/cssc.201600987

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarmad, S., Xie, Y. J., Mikkola, J.-P., and Ji, X. Y. (2017b). Screening of deep eutectic solvents (DESs) as green CO₂ sorbents: from solubility to viscosity. New J. Chem. 41, 290–301. doi: 10.1039/C6NJ03140D

CrossRef Full Text | Google Scholar

Shukla, S. K., and Mikkola, J.-P. (2018). Intermolecular interactions upon carbon dioxide capture in deep-eutectic solvents. Phys. Chem. Chem. Phys. 20, 24591–24601. doi: 10.1039/c8cp03724h

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, Y. J., Dong, H. F., Zhang, S. J., Lu, X. H., and Ji, X. Y. (2013). “Thermophysical properties and gas solubilities in choline chloride/urea for CO2 separation,” in ICAE 2013 (Pretoria, South Africa).

Google Scholar

Zubeir, L. F., Lacroix, M. H. M., and Kroon, M. C. (2014). Low transition temperature mixtures as innovative and sustainable CO₂ capture solvents. J. Phy. Chem. B 118, 14429–14441. doi: 10.1021/jp5089004

PubMed Abstract | CrossRef Full Text | Google Scholar

Zubeir, L. F., van Osch, D. J. G. P., Rocha, M. A. A., Banat, F., and Kroon, M. C. (2018). Carbon dioxide solubilities in decanoic acid-based hydrophobic deep eutectic solvents. J. Chem. Eng. Data 63, 913–919. doi: 10.1021/acs.jced.7b00534

PubMed Abstract | CrossRef Full Text | Google Scholar