1. Introduction
As reported by the Energy Information Agency (EIA) in 2023, approximately 83% of the total energy produced in the United States came from fossil fuels, i.e., oil, coal, and natural gas [
1]. Fossil fuels, when burned, produce large quantities of CO
2 and other gases. These gases will result in greenhouse gas effects, resulting in the absorption of the heat coming from the sun and causing the temperature of the earth’s atmosphere to rise. Global warming is consequently considered a detrimental and even catastrophic change in the climate. According to the Scripps Institution of Oceanography and the National Oceanic and Atmospheric Administration (NOAA), the CO
2 May monthly average in 2024 reached 426.7 ppm, an increase of 2.92 ppm over that recorded in 2023 [
2]. NOAA also reports that 2023 is the warmest year on record by far [
3].
There are multiple ways of slowing down the rate of global warming. One of those methods is through emission reduction with carbon capture. The most up-to-date technology in carbon capture is the process of absorbing CO2 generated by power generation, steel, and cement plants. This latest technology is similar to the removal of acid gases from natural gas sweetening, except that a mixture of amines is used as a solvent. Capturing CO2 using non-aqueous solvents is an advanced, next-generation post-combustion capture technology that seeks to remove CO2 emitted from power generation and industrial flue gas streams.
Natural gas is one of the types of energy from the diverse portfolio of fossil fuels. It is a naturally occurring hydrocarbon composed mainly of methane. In general, it also contains ethane, propane, butane, and pentane, as well as non-hydrocarbon gases such as nitrogen, carbon dioxide, and hydrogen sulfide. Natural gas is formed from the remains of plants and animals exposed to extreme heat and pressure under the earth’s surface hundreds of millions of years ago. It constitutes the major source of energy consumed today; the demand for natural gas has not ceased to increase year after year in the last decade. According to the International Energy Agency, the world’s natural gas consumption reached 4138 billion cubic meters (bcm) in 2022 and is expected to reach 4299 bcm in 2050 [
4].
To this day, carbon capture is performed by capturing CO
2 from both high-pressure natural gas streams and low-pressure gas streams with amines-based chemical solvents [
5]. However, this process is not very economical in the case of flue gases because of the large amount of energy required to regenerate the solvent and the serious solvent degradation that occurs. For natural gas applications, using physical rather than chemical solvents for carbon capture is more economical when the cleaning targets are not at the ppm level. Multiple processes have been used for carbon capture depending on the partial pressure of CO
2 in the feed and the required degree of treatment. Aside from amines and their blends, organic physical solvents are one of these processes. One of the main advantages is the low energy needed for solvent regeneration [
5]. Commonly, physical solvents have a low energy of regeneration and relatively low vapor pressure. They are non-corrosive, with high thermal stability. They could be used in treating natural gas as standalone or mixed with aqueous amines in the case of flue gases. Ionic liquids (ILs) represent a category of physical solvents that utilize intermolecular forces for CO
2 capture. ILs are normally defined as compounds composed of ions (a cation and an anion) with a melting point below 100 °C. In this study, three new and promising ionic liquids are studied to analyze their capacities for ethane at three different temperatures (303.15 K, 323.15 K, and 343.15 K). This range of temperatures is of great importance to the industry and the pressures used were up to 1.4 MPa. Using CO
2 solubility data, we recently published, selectivity values are reported for the three solvents at three temperatures. The experiments were conducted using a gravimetric microbalance (IGA-003).
2. Results
2.1. Verification of the C2H6 Solubility Measurement
A verification test was performed to ensure the IGA-003 is operating as it should. In this verification test, [BMIM][PF6] was used as the solvent, and the results were compared with Anthony et al. (2005) at 298.15 K [
6]. As shown in
Figure 1, the results between the measured data and those published by Anthony et al. are in excellent agreement.
2.2. Solubility of C2H6 in Ionic Liquids
C
2H
6 solubility data were studied by many researchers who found the ionic liquids to have lower solubility than CO
2. This work used three ionic liquids with different cations and anions to understand the factors affecting the solubility of C
2H
6 in ILs. C
2H
6 solubility in [HMIM][Tf2N], [BMIM][DMP], and [PMIM][Tf2N] were measured at 303.15 K, 323.15 K, and 343.15 K and pressures up to 1.4 MPa. The results were plotted in P-
x form in
Figure 2,
Figure 3 and
Figure 4 for the three solvents. The values of the experimental data can be found in
Table 1,
Table 2 and
Table 3.
3. Discussion
3.1. Factors That Affect the Solubility of C2H6
3.1.1. Pressure and Temperature
Figure 2,
Figure 3 and
Figure 4 show that temperature and pressure affect the solubility of C
2H
6 in ionic liquids. As expected, it is observed that as the temperature rises, C
2H
6 solubility in the ionic liquid decreases. As pressure rises, C
2H
6 solubility increases.
As presented in
Figure 5, a comparison between the solubility of ethane in the three ionic liquids at 303.15 K shows that the solvent that absorbs ethane the least is [BMIM][DMP].
3.1.2. Anion Effect
A significant amount of research has been performed to understand the mechanism of how different ionic liquids have different solubilities when it comes to capturing gases such as CO
2, CH
4, C
2H
6, and others. Anthony et al. [
6] studied the anion effect on CO
2 solubility in [BMIM]-based ionic liquids with three different anions. Experiments were performed between 283.15K and 323.15 K and up to 1.4 MPa for three ionic liquids with three different anions, namely [BMIM][Tf2N], [BMIM][BF4], and [BMIM][PF6]. They reported that ethane, along with other gases, such as ethylene, oxygen, and argon, have the lowest solubilities in ILs compared to other common solvents.
3.1.3. Alkyl Chain Length with [Tf2N]-Anion
As presented in
Figure 6, the alkyl chain effect will affect the solubility of C
2H
6 in ionic liquids. [HMIM] and [PMIM] are attached to an imidazolium-based cation with an n-methyl chain. The difference between their molecular structures is due to the hexyl and propyl alkyl chains. In their carbon chains, hexyl has six methyl groups, and propyl has three. The longer the chain, the more gas the ionic liquid will absorb. This observation was made also by Zoubeik (2014) [
7] and Tagiuri et al. (2014) [
8] based on several studies published in the literature. Note that the polarity parameter (N) for [PMIM] is lower than that of [HMIM], confirming its higher capacity in absorbing ethane. The polarity parameter is discussed further at the end of
Section 3.1.4 and also in reference [
9].
3.1.4. Cation with the Same [Tf2N] Anion
The effect of different cations with the same [Tf2N] anion is observed in
Figure 6. The [PMIM] cation allowed for less C
2H
6 solubility than the [HMIM] cation. The solubility in ILs with the same anions will depend mostly on the cation. Based on studies published in the literature, methyl imidazolium cation-based ionic liquid tends to usually have a higher C
2H
6 solubility than other cation-based ionic liquids. The solubility of C
2H
6 in [HMIM] is higher than that in [PMIM], most probably because of its longer carbon (alkyl) chain, as discussed earlier in the section entitled “Alkyl Chain Lengths”. As reported by Taguiri (2019) [
10] and others, this study confirms the fact that the anion of an ionic liquid has more of an impact on gas solubility than its cation.
For detailed information about the mechanism of gas solubility in ionic liquids, please refer to a comprehensive screening study we performed in 2011 [
9] using COSMO-RS. In it, molar volumes and polarity of the ionic liquids were used to quantify the solubility of gases in ILs. The study dealt with gases such as CO
2 and the selectivity of methane and nitrogen. A relative overall polarity parameter N was defined and computed for 2701 ionic liquids. N quantifies the degree of chemical affinity of water with ions due to vdW, misfit, and H-bond. The more negative the value of N, the more polar the IL is. N was found to be equal to 6.22 for [Tf2N], −62.04 for [DMP], 37.16 for [HMIM], 32.79 for [BMIM], and 32.70 for a cation similar to [PMIM]. Accordingly, [BMIM][DMP] with the lowest N should have the lowest ethane solubility, followed by [PMIM][Tf2N] and [HMIM][Tf2N]. This predicted trend is confirmed by the experimental data presented in
Figure 5. A comprehensive description of the structural variations of the ILs on the solubility can also be found. In addition to the alkyl chain, alkylation of ammonium and phosphonium cation, change in cation family for ring cation and non-ring cation, and presence of hydroxyl and ether groups were discussed in reference [
9].
3.2. Optimization of Binary Interaction Parameters
MATLAB was used for the thermodynamic model and was based on a bubble point technique for the three mixing rules [
5]. As shown in Equation (1) [
10], the binary interaction parameters were optimized by minimizing the objective function (Err). Because experimental data at lower pressures are more likely to be erroneous, the binary interaction parameters were optimized using a pressure range of 0.1 to 1.4 MPa.
The average absolute deviation for the vdW1, vdW2, and WS-NRTL mixing rules when used in the ionic liquids for C
2H
6 were 4.39%, 2.45%, and 2.45%, respectively. The results for C
2H
6 were best when using either the vdW2 or WS-NRTL mixing rule, as shown in
Table 4.
3.3. Henry’s Law Constants, Enthalpies, and Entropies
Using the solubility data, Henry’s Law constants were derived for the three ionic liquids. For each temperature, the fugacities of C
2H
6 were plotted against mole fractions. The limiting slope of the data’s second-order polynomial was used to determine Henry’s Law constants for each temperature. Henry’s Law constants
of the solute (i) in the solvent (j) is the ratio of the fugacity of the solute (i) to the mole fraction of the solute (i) in solvent (j) at infinite dilution at temperature T and pressure P;
and
are the fugacity of solute (i) in the liquid phase and vapor phase, respectively; and
and
are the mole fraction of solute (i) in the liquid phase and vapor phase, respectively, at temperature T and pressure P [
5]. It is expressed as the following:
Once Henry’s Law has been estimated from the solubility data, estimating the enthalpy of solvation (ΔH
∞) at infinite dilution and the entropy of solvation (ΔS
∞) at infinite dilution can also be performed using the following equations:
where ΔH
∞ is the enthalpy of solvation at infinite dilution, (ΔS
∞) is the entropy of solvation at infinite dilution, R is the universal gas constant, H is Henry’s Law constant, and T is temperature.
The comparison between the three ionic liquids is presented in
Figure 7 for the three temperatures studied. Henry’s Law constants for C
2H
6 in [HMIM][Tf2N] were published by C. Gomes et al. (2007) [
11] and are reported here for comparison. The values are very close, confirming the accuracy of the experimental data. The calculated values are presented in
Table 5. Henry’s Law constants of the three ionic liquids studied are compared to those found in the literature and presented in
Figure 7 for the three temperatures studied.
Overall, the ranking of the Henry’s Law constants for C2H6 is the following: [BMIM][DMP] > [PMIM][Tf2N] > [HMIM][Tf2N] at all temperatures. As the highest value is the best, and [BMIM][DMP] is therefore the IL that absorbs ethane the least. Compared to other ionic liquids, the three ionic liquids are not the best, but [BMIM][DMP] seems to be the best option.
The decreasing ranking in terms of the heat of absorption is the following: [BMIM][DMP] > [HMIM][Tf2N] > [PMIM][Tf2N]. The solubility of ethane in three ILs {[C3C1pip][FSI], [C3C1pyrr][FSI], and [N1223][FSI]} reported in
Figure 7 was published by Nath et al. (2017) [
12].
3.4. CO2/C2H6 Selectivity in the Three Ionic Liquid Systems
The major criteria for a potential solvent for the use of gas solubility are the following: high removal efficiency, high selectivity toward CO
2, and low energy requirement. Not much research has been conducted when it comes to the selectivity between gases. When CO
2 is captured, it is not the only gas being absorbed. Other gases are absorbed, such as N
2O, O
2, SO
2, H
2S, CH
4, and many other hydrocarbons. Therefore, with that in mind, it is crucial to study the effect of selectivity on ionic liquids. The selectivity of the ionic liquids used in this work was estimated as the ratio between Henry’s Law constant of CO
2 and C
2H
6. The equation used to estimate the selectivity is the following:
The values of the selectivities are presented in
Table 6. CO
2/C
2H
6 selectivity of the ionic liquids is not significantly high; however, they decrease as the temperature rises. The data obtained for the Henry’s Law constant of CO
2 was taken from previous research we performed measuring the CO
2 solubility in [HMIM][Tf2N], [BMIM][DMP], and [PMIM][Tf2N] [
13]. The study revealed that the three solvents were among the best in capturing CO
2. The three ionic liquids followed this rank in terms of increasing CO
2 solubility: [HMIM] [Tf2N] > [PMIM][Tf2N] >[BMIM][DMP], as reported in reference [
13]. The study confirmed that [HMIM][Tf2N] is an excellent solvent for CO
2 capture. [PMIM][Tf2N] was also found to have a high capacity for CO
2 among ILs published in the literature.
The selectivities of the three ionic liquids studied are compared to those found in the literature and presented in
Figure 8,
Figure 9 and
Figure 10 at the three temperatures studied.
For a comprehensive comparison, we have added in the table data related to Selexol, a widely used and efficient physical solvent (a mixture of polyethylene glycol dimethyl ethers). Selectivity values were calculated based on values reported by Rayer et al. (2012) [
14].
Overall, the ranking of the CO
2/C
2H
6 selectivity is the following: [PMIM][Tf2N] > [BMIM] [DMP] > [HMIM][Tf2N] > Selexol. As the highest selectivity is the one preferred, among the three ionic liquids studied, [PMIM][Tf2N] is considered the best. As shown below, these ILs do not possess the best selectivities among ionic liquids in the literature. Selexol seems to have the lowest selectivity compared to the ILs selected in this study. More details about the solvents in
Figure 8,
Figure 9 and
Figure 10 can be found in reference [
12].
4. Materials and Methods
4.1. Materials
Table 7 lists the structures, acronyms, suppliers, purities, and water content of all of the solvents and solutes used in this study. The manufacturers of the chemicals provided the information related to the purities of the solvents listed in the same table.
4.2. Solubility Measurement
The gas solubility experiments were performed with the use of an IGA-003 (Intelligent Gravimetric Analyzer, Hidden Isochema Ltd., Warrington, UK). The machine measures the weight change of a liquid sample once a gas/gas mixture has been injected into the system. The IGA-003 can run at a maximum pressure and temperature of up to 2.0 MPa and 773.15 K, respectively.
For each isotherm experiment, around 50 to 90 mg of ionic liquid sample was inserted into the sample container to determine the solubility of the gas in the ionic liquid. A Grant water bath was used to regulate the temperature of the sample. For each experiment, the sample was degassed for around 10 h at high temperatures and under a deep vacuum before beginning the isotherm process to remove any contaminants that may have been in the sample. Throughout the degassing process, the temperature of the sample was kept constant at 349.15 K.
Once the sample weight stabilized after being degassed for over 10 h, the operating temperature was reduced to the isotherm condition using a Grant water bath. The ionic liquid sample was kept under this condition for at least an hour to ensure that the sample weight was at equilibrium. Once the weight of the sample has been stable, the isotherm process can be initiated by choosing a list of pressures at which the gas solubility is to be measured. Once all the necessary parameters have been selected in the IGA computer software (Version 2010), C2H6 is injected into the chamber using the mass flow controller (MFC). The system’s pressure is kept constant at a pre-set pressure until the system reaches equilibrium. The IGA-003 records all real-time data, such as weight, temperature, and pressure.
Measurement of the density of the solvent is a required input in the experimental procedure. We reported the values of the densities of the 3 solvents in reference [
13]. The density values decreased in the following order: [PMIM][Tf2N] > [HMIM][Tf2N] > [BMIM][DMP]. [BMIM][DMP] was reported to have a high viscosity value of 409.88 mPa·s at 303.15 K [
13].
Figure 11 shows the setup for the IGA. Each instrument used in this work was numerically labelled. The setup consists of (1) a computer, (2) a microbalance, (3) a gas cylinder, (4) mass flow control, and (5) a water bath.
4.3. Thermodynamic Modeling
The solubility data of C
2H
6 in ILs obtained in this work were correlated using the Peng–Robinson equation of state (PR-EoS), shown in Equation (6) [
15], where
is the intermolecular attractive force and
is the van der Waals co-volume factor. The PR-EoS combined with three different mixing rules was used to build this model. The three mixing rules are the following:
Single binary interaction parameter based on van der Waals one (vdW1);
Two binary interaction parameters based on van der Waals two (vdW2);
Wong–Sandler mixing rules combined with the NRTL model (WS-NRTL).
4.3.1. Van der Waals Mixing Rule
For this case, van der Waals created two different mixing rules to estimate mixture parameters (
and
): van der Waals one (vdW1) and van der Waals two (vdW2) [
15,
16,
17]. In the vdW1 mixing rule, only one parameter (lij) is regressed, whereas, for vdW2, two binary interaction parameters are obtained (lij and kij). The parameter a
m for both vdW1 and vdW2 mixing rules was estimated using Equation (7). Using the temperature dependence binary interaction parameters, kij and aij can be estimated with Equation (10). The co-volume factor bm for both vdW1 and vdW2 is estimated with Equation (8) and Equation (9), respectively. Using the interaction parameter (lij),
can be estimated with Equation (11).
where
where
and
are the pure component intermolecular attractive force parameter and
and
are the pure components co-volume factor parameters for the components (
and
).
4.3.2. Wong–Sandler Mixing Rule
The mixture’s attractive force parameter, a, and co-volume parameter, b, can be calculated using Equations (12) and (13) [
15], respectively, for Wong–Sandler mixing rules. In this paper, the Non-Random Two-Liquid (NRTL) model is utilized to approximate the activity coefficient and Excess Gibbs Energy (Equations (18)–(24)) [
16]. The mixture parameters were estimated using three binary interaction parameters (
,
and
), where
and
are the NRTL parameters. The value of
in Equations (19) and (20) were set at 0.3 arbitrarily in this work. The energies
and
are the energies of interaction between molecules ‘
’ and ‘
’ [
17,
18,
19] and are calculated using Equations (21) and (22).
where
where
in Equation (14) calculated using Equation (17):
4.4. Critical Properties Calculations
The critical parameters of the gas and the solvent examined must be known to be used in the model. The modified Lydersen–Joback–Reid group contribution method was used to determine the critical temperature (Tc), critical pressure (Pc), and the acentric factor (
) of the ionic liquids in this study using the spreadsheet provided by Valderrama and Rojas [
20].
Table 8 lists the critical properties determined for this study’s three ionic liquids and C
2H
6.
5. Conclusions
The purpose of this work was to screen three new ionic liquids, measuring their capacity to absorb C2H6 at pressures up to 1.4 MPa and at temperatures ranging from 293.15 K to 343.15 K. Thermodynamic properties such as the enthalpies and entropies of absorption are reported and the data used to determine the selectivity values.
The ionic liquids selected were [HMIM][Tf2N], [BMIM][DMP], and [PMIM][Tf2N]. In an earlier publication, we reported that the density values of the solvents decreased in the following order: [PMIM][Tf2N] > [HMIM][Tf2N] > [BMIM][DMP].
As co-absorption of ethane is important in natural gas processing, Henry’s Law constants for C
2H
6 were estimated at 303.15 K, 323.15 K, and 343.15 K. As expected, they increased for all ILs as temperature increased. Overall, the ranking of the Henry’s Law constants for C
2H
6 is the following: [BMIM][DMP] > [PMIM][Tf2N] > [HMIM][Tf2N] at all temperatures. [BMIM][DMP] is, therefore, found to be the IL that absorbs ethane the least. Compared to other ionic liquids, the three ionic liquids are not the best, but [BMIM][DMP] seems to be the best option. The order in the selectivity can be predicted based on the overall polarity number N and its calculated values obtained from reference [
9].
The overall ranking of the CO2/C2H6 selectivity is the following: [PMIM][Tf2N] > [BMIM] [DMP] > [HMIM][Tf2N] > Selexol. The most promising ionic liquid in terms of selectivity is therefore [PMIM][Tf2N]. Note that the three ionic liquids studied had better selectivity than Selexol, a widely used physical solvent in gas sweetening plants.
Solubility data were correlated using the Peng–Robinson equation of state (PR-EoS), and the binary interaction coefficients were reported. The average absolute deviations for the vdW1, vdW2, and WS-NRTL mixing rules when used in the ionic liquids for C2H6 were 4.39%, 2.45%, and 2.45%, respectively. The enthalpy and entropies of adsorption were calculated. Negative enthalpies were reported, indicating an exothermic process of absorption.
Author Contributions
Conceptualization, A.H.; methodology, N.H.; software, N.H.; validation, N.H.; formal analysis, N.H.; investigation, N.H.; resources, A.H.; data curation, N.H.; writing—original draft preparation, N.H.; writing—review and editing, A.H. and H.I.; visualization, N.H.; supervision, H.I. and A.H.; project administration, A.H.; funding acquisition, A.H. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by the Natural Sciences and Engineering Research Council of Canada (NSERC) under a Discovery Grant number (RGPIN-06805-2018).
Data Availability Statement
All experimental data acquired are reported in the manuscript.
Acknowledgments
The first author would like to thank Devjyoti Nath for training him on the microbalance instrument and for his help in modeling the data.
Conflicts of Interest
The authors declare no conflicts of interest. The funding agency had no role in the design of the study; in the collection, analysis, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.
Nomenclature
AAD | Average absolute deviation |
[BF4] | Tetrafluoroborate |
[BMIM][DMP] | 1-Butyl-3-methyl-imidazolium dimethyl-phosphate |
[N1223][FSI] | N-propyl-n-methyl-n,n-dimethyl-ammonium bis(fluorosulfonyl)imide |
[C3C1PIP][FSI] | N-propyl-n-methyl-piperidinium bis(fluorosulfonyl)imide |
[C3C1Pyrr][FSI] | 1-Methyl-1-propylpyrrolidinium bis(fluorosulfonyl)imide |
[C4C1Im] [BF4] | 1-Butyl-3-methylimidazolium tetrafluoroborate |
CO2 | Carbon dioxide |
C2H6 | Ethane |
DCN | Dicyanamide |
EoS | Equation of state |
[FSI] | Bis(fluorosulfonyl)imide |
H | Henry’s Law constant |
[HMIM][Tf2N] | 1-Hexyl-3-methylimidazolium bis(trifluormethylsulfonyl)imide |
IGA | Intelligent Gravimetric Analyzer |
ILs | Ionic liquids |
K | Kelvin |
MDEA | Methyldiethanolamine |
MEA | Monoethanolamine |
Pc | Critical pressure |
PF6 | Hexafluoro-phosphate |
[PMIM][Tf2N] | 1-Propyl-3-methylimidazolium bis(trifluoromethyl-sulfonyl)-imide |
PR | Peng–Robinson |
RSO4 | Alkylsulfate |
Tb | Boiling point temperature |
Tc | Critical temperature |
vdW1 | van der Waals one |
vdW2 | van der Waals two |
WS-NRTL | Wong–Sandler mixing rules + Non-Random Two-Liquid model |
x | Liquid mole fraction (moles gas/moles gas + moles solVent) |
y | Vapor mole fraction |
Zc | Critical compressibility factor |
References
- Energy Information Administration (EIA). U.S. Primary Energy Consumption by Energy Source. 2024. Available online: https://www.eia.gov/energyexplained/us-energy-facts/ (accessed on 22 July 2024).
- National Oceanic and Atmospheric Administration. During a Year of Extremes, Carbon Dioxide Levels Surge Faster than Ever. 2024. Available online: https://www.noaa.gov/news-release/during-year-of-extremes-carbon-dioxide-levels-surge-faster-than-ever (accessed on 22 July 2024).
- National Oceanic and Atmospheric Administration. 2023 Was the World’s Warmest Year on Record, by Far. 2024. Available online: https://www.ncei.noaa.gov/access/monitoring/monthly-report/global/202313#:~:text=The%20year%202023%20was%20the,decade%20(2014%E2%80%932023) (accessed on 22 July 2024).
- International Energy Agency (IEA). World Energy Outlook. 2023. Available online: https://iea.blob.core.windows.net/assets/86ede39e-4436-42d7-ba2a-edf61467e070/WorldEnergyOutlook2023.pdf (accessed on 22 July 2024).
- Nath, D.; Henni, A. Solubility of carbon dioxide (CO2) in four bis (trifluoromethyl-sulfonyl)imide ([Tf2N]) based ionic liquids. Fluid Phase Equilibria 2020, 524, 112757. [Google Scholar] [CrossRef]
- Anthony, J.L.; Anderson, E.J.; Maginn, E.J.; Brennecke, J.F. Anion effects on gas solubility in ionic liquids. J. Phys. Chem. B 2005, 109, 6366–6374. [Google Scholar] [CrossRef] [PubMed]
- Zoubeik, M.; Henni, A. Experimental and thermodynamic study of CO2 solubility in promising [Tf2N and DCN] ionic liquids. Fluid Phase Equilibria 2014, 376, 22–30. [Google Scholar] [CrossRef]
- Tagiuri, A.; Sumon, K.Z.; Henni, A. Solubility of carbon dioxide in three [Tf2N] ionic liquids. Fluid Phase Equilibria 2014, 380, 39–47. [Google Scholar] [CrossRef]
- Sumon, Z.K.; Henni, A. Ionic liquids for CO2 capture using COSMO-RS: Effect of structure, properties and molecular interactions on solubility and selectivity. Fluid Phase Equilibria 2011, 310, 39–55. [Google Scholar] [CrossRef]
- Tagiuri, A. Studies of Solubility of CO2 in Ionic Liquids, Kinetics, and Heat of Reactions in Promising Cyclic Amines. Ph.D. Thesis, University of Regina, Regina, SK, Canada, 2019. [Google Scholar]
- Costa Gomes, M.F. Low-pressure solubility and thermodynamics of solvation of carbon dioxide, ethane, and hydrogen in 1-hexyl-3-methylimidazolium bis(trifluoromethylsulfonyl)amide between temperatures of 283 K and 343 K. J. Chem. Eng. Data 2007, 52, 472–475. [Google Scholar] [CrossRef]
- Nath, D.; Sumon, S.K.; Tagiuri, A.; Henni, A. Effect of cation on the solubility of ethane in three bis(fluorosulfonyl)imide ([FSI]) based low viscosity ionic liquids. Fluid Phase Equilibria 2017, 454, 78–90. [Google Scholar] [CrossRef]
- Henni, N.; Henni, A.; Ibrahim, H. Solubility of carbon dioxide in promising methylimidazolium-based ionic liquids. Fluid Phase Equilibria 2023, 565, 13619. [Google Scholar] [CrossRef]
- Rayer, A.V.; Henni, A.; Tontiwachwuthikul, P. High-Pressure Solubility of Methane (CH4) and Ethane (C2H6) in Mixed Polyethylene Glycol Dimethyl Ethers (Genosorb 1753) and Its Selectivity in Natural Gas Sweetening Operations. J. Chem. Eng. Data 2012, 57, 764–775. [Google Scholar] [CrossRef]
- Peng, D.Y.; Robinson, D.B. A New Two-Constant Equation of State. Ind. Eng. Chem. Fundam. 1976, 15, 59–64. [Google Scholar] [CrossRef]
- Renon, H.; Prausnitz, J.M. Local compositions in thermodynamic excess functions for liquid mixtures. AIChE J. 1968, 14, 135–144. [Google Scholar] [CrossRef]
- Orbey, H.; Sandler, S.I. Modeling Vapor-Liquid Equilibria: Cubic Equations of State and Their Mixing Rules; Cambridge University Press: Cambridge, UK, 1998. [Google Scholar]
- Eslamimanesh, A.; Mohammadi, A.H.; Salamat, Y.; Shojaei, M.-J.; Eskandari, S.; Richon, D. Phase behavior of mixture of supercritical CO2+ionic liquid: Thermodynamic consistency test of experimental data. AIChE J. 2013, 59, 3892–3913. [Google Scholar] [CrossRef]
- Lee, M.-T.; Lin, S.-T. Prediction of mixture vapor–liquid equilibrium from the combined use of Peng–Robinson equation of state and COSMO-SAC activity coefficient model through the Wong–Sandler mixing rule. Fluid Phase Equilibria 2007, 254, 28–34. [Google Scholar] [CrossRef]
- Valderrama, J.O.; Rojas, R.E. Critical properties of ionic liquids. Revisited. Ind. Eng. Chem. Res. 2009, 48, 6890–6900. [Google Scholar] [CrossRef]
Figure 1.
Comparison between the measured data for the solubility of C
2H
6 in [BMIM][PF6] and those reported by Anthony et al. (2005) [
6] at 298.15 K.
Figure 1.
Comparison between the measured data for the solubility of C
2H
6 in [BMIM][PF6] and those reported by Anthony et al. (2005) [
6] at 298.15 K.
Figure 2.
Solubility data for C2H6 in [HMIM][Tf2N] at different temperatures and pressures up to 1.4 MPa: ♦ 303.15 K; ■ 323.15 K; ▲: 343.15 K. Experimental VLE data were correlated with (a) PR + vdW1, (b) PR + vdW2, and (c) PR + WS + NRTL.
Figure 2.
Solubility data for C2H6 in [HMIM][Tf2N] at different temperatures and pressures up to 1.4 MPa: ♦ 303.15 K; ■ 323.15 K; ▲: 343.15 K. Experimental VLE data were correlated with (a) PR + vdW1, (b) PR + vdW2, and (c) PR + WS + NRTL.
Figure 3.
Solubility data of C2H6 in [BMIM][DMP] at different temperatures: ♦ 303.15 K; ■ 323.15 K; ▲ 343.15 K. Experimental VLE data were correlated with (a) PR + vdW1, (b) PR + vdW2, and (c) PR + WS + NRTL.
Figure 3.
Solubility data of C2H6 in [BMIM][DMP] at different temperatures: ♦ 303.15 K; ■ 323.15 K; ▲ 343.15 K. Experimental VLE data were correlated with (a) PR + vdW1, (b) PR + vdW2, and (c) PR + WS + NRTL.
Figure 4.
Solubility data of C2H6 in [PMIM][Tf2N] at different temperatures: ♦ 303.15 K; ■ 323.15 K; ▲ 343.15 K and model correlation: lines with (a) PR + vdW1, (b) PR + vdW2, and (c) PR + WS + NRTL.
Figure 4.
Solubility data of C2H6 in [PMIM][Tf2N] at different temperatures: ♦ 303.15 K; ■ 323.15 K; ▲ 343.15 K and model correlation: lines with (a) PR + vdW1, (b) PR + vdW2, and (c) PR + WS + NRTL.
Figure 5.
Comparison between the solubility of C2H6 in the 3 ILs at 303.15 K.
Figure 5.
Comparison between the solubility of C2H6 in the 3 ILs at 303.15 K.
Figure 6.
Solubility of C2H6 in the two ILs with [Tf2N]-anion at 303.15 K.
Figure 6.
Solubility of C2H6 in the two ILs with [Tf2N]-anion at 303.15 K.
Figure 7.
Comparison of Henry’s Law constants for C
2H
6 in ILs at (
a) 303.15 K, (
b) 323.15 K, and (
c) 343.15 K; red: ionic liquids used in this work and blue: other ILs obtained from the literature summarized by Nath et al. [
12].
Figure 7.
Comparison of Henry’s Law constants for C
2H
6 in ILs at (
a) 303.15 K, (
b) 323.15 K, and (
c) 343.15 K; red: ionic liquids used in this work and blue: other ILs obtained from the literature summarized by Nath et al. [
12].
Figure 8.
Comparison of CO
2/C
2H
6 selectivity in several ionic liquids reported by Nath et al. [
12] at 303.15 K.
Figure 8.
Comparison of CO
2/C
2H
6 selectivity in several ionic liquids reported by Nath et al. [
12] at 303.15 K.
Figure 9.
Comparison of CO
2/C
2H
6 selectivity in several ionic liquids reported by Nath et al. [
12] at 323.15 K.
Figure 9.
Comparison of CO
2/C
2H
6 selectivity in several ionic liquids reported by Nath et al. [
12] at 323.15 K.
Figure 10.
Comparison of CO
2/C
2H
6 selectivity with other ionic liquids reported by Nath et al. [
12] at 343.15 K.
Figure 10.
Comparison of CO
2/C
2H
6 selectivity with other ionic liquids reported by Nath et al. [
12] at 343.15 K.
Figure 11.
IGA-003 setup.
Figure 11.
IGA-003 setup.
Table 1.
Solubility of C2H6 in [HMIM][Tf2N] at 303.15 K, 323.15 K, and 343.15 K.
Table 1.
Solubility of C2H6 in [HMIM][Tf2N] at 303.15 K, 323.15 K, and 343.15 K.
[HMIM][Tf2N] |
---|
303.15 K | 323.15 K | 343.15 K |
---|
| Pressure (MPa) | | Pressure (MPa) | | Pressure (MPa) |
---|
0.003 | 0.01991 | 0.003 | 0.0199 | 0.003 | 0.0199 |
0.006 | 0.0499 | 0.005 | 0.0499 | 0.005 | 0.0499 |
0.013 | 0.09989 | 0.010 | 0.1000 | 0.008 | 0.0999 |
0.025 | 0.20004 | 0.021 | 0.1998 | 0.020 | 0.2002 |
0.049 | 0.3998 | 0.041 | 0.3997 | 0.037 | 0.3998 |
0.071 | 0.59988 | 0.059 | 0.5998 | 0.054 | 0.5998 |
0.092 | 0.7999 | 0.077 | 0.7998 | 0.070 | 0.7995 |
0.111 | 0.9994 | 0.094 | 0.9996 | 0.081 | 0.9999 |
0.120 | 1.0998 | 0.102 | 1.0998 | 0.091 | 1.1000 |
0.129 | 1.1988 | 0.110 | 1.2011 | 0.103 | 1.1994 |
0.135 | 1.2495 | 0.114 | 1.2503 | 0.107 | 1.2513 |
0.139 | 1.3007 | 0.119 | 1.3009 | - | - |
Table 2.
Solubility of C2H6 in [BMIM][DEP] at 303.15 K, 323.15 K, and 343.15 K.
Table 2.
Solubility of C2H6 in [BMIM][DEP] at 303.15 K, 323.15 K, and 343.15 K.
[BMIM][DMP] |
---|
303.15 K | 323.15 K | 343.15 K |
---|
| Pressure (MPa) | | Pressure (MPa) | | Pressure (MPa) |
---|
0.001 | 0.01978 | 0.001 | 0.0198 | 0.002 | 0.0196 |
0.002 | 0.0499 | 0.002 | 0.0499 | 0.002 | 0.0499 |
0.004 | 0.09989 | 0.004 | 0.0999 | 0.004 | 0.1001 |
0.008 | 0.19985 | 0.010 | 0.1999 | 0.010 | 0.1999 |
0.017 | 0.39991 | 0.018 | 0.4002 | 0.017 | 0.3996 |
0.026 | 0.59976 | 0.028 | 0.6006 | 0.023 | 0.5999 |
0.034 | 0.79991 | 0.038 | 0.7998 | 0.032 | 0.7999 |
0.0422 | 0.9998 | 0.045 | 0.9998 | 0.043 | 1.0006 |
0.050 | 1.2003 | 0.054 | 1.0997 | 0.047 | 1.0998 |
0.056 | 1.2996 | 0.059 | 1.2005 | 0.053 | 1.2012 |
0.059 | 1.3486 | - | - | 0.059 | 1.3011 |
0.063 | 1.4006 | - | - | 0.065 | 1.4006 |
Table 3.
Solubility of C2H6 in [PMIM][Tf2N] at 303.15K, 323.15K, and 343.15K.
Table 3.
Solubility of C2H6 in [PMIM][Tf2N] at 303.15K, 323.15K, and 343.15K.
[PMIM][Tf2N] |
---|
303.15 K | 323.15 K | 343.15 K |
---|
| Pressure (MPa) | | Pressure (MPa) | | Pressure (MPa) |
---|
0.002 | 0.0198 | 0.003 | 0.0199 | 0.002 | 0.0198 |
0.005 | 0.0499 | 0.005 | 0.0499 | 0.002 | 0.0500 |
0.009 | 0.0999 | 0.008 | 0.0999 | 0.005 | 0.0999 |
0.016 | 0.1999 | 0.015 | 0.1999 | 0.011 | 0.1999 |
0.032 | 0.4001 | 0.029 | 0.3999 | 0.023 | 0.4002 |
0.045 | 0.5999 | 0.042 | 0.5997 | 0.033 | 0.6011 |
0.058 | 0.7999 | 0.053 | 0.7998 | 0.044 | 0.8000 |
0.070 | 0.9999 | 0.065 | 0.9998 | 0.054 | 0.9999 |
0.083 | 1.1991 | - | - | 0.065 | 1.2012 |
0.088 | 1.30132 | - | - | 0.070 | 1.2998 |
- | - | - | - | 0.079 | 1.4000 |
Table 4.
Optimized binary interaction parameters for different models and calculated AAD%.
Table 4.
Optimized binary interaction parameters for different models and calculated AAD%.
Compound | Temperature (K) and Pressure (MPa) | Mixing Rule | a | b | | | AAD% |
---|
C2H6(1) + [HMIM][Tf2N] (2) | 303.15 and 0.1–1.4 | vdW1 | 0.0797 | - | - | - | 3.19 |
| | vdW2 | 0.1566 | 0.0163 | - | - | 0.47 |
| | WS + NRTL | 1.305 | - | 0.2151 | −0.1848 | 0.44 |
| 323.15 and 0.1–1.4 | vdW1 | 0.0694 | - | - | - | 2.11 |
| | vdW2 | 0.1116 | 0.0104 | - | - | 0.83 |
| | WS + NRTL | 1.2293 | - | 0.7031 | −0.5823 | 0.84 |
| 343.15 and 0.1–1.4 | vdW1 | 0.0452 | - | - | - | 3.37 |
| | vdW2 | 0.0557 | 0.0026 | - | - | 3.30 |
| | WS + NRTL | 1.0256 | | 0.0396 | 0.0123 | 3.30 |
C2H6 (1) + [BMIM][DMP] (2) | 303.15 and 0.1–1.4 | vdW1 | 0.1931 | - | - | - | 4.77 |
| | vdW2 | 0.052 | −0.0377 | | - | 2.35 |
| | WS + NRTL | 0.9997 | - | 1.1024 | 0.6362 | 2.34 |
| 323.15 and 0.1–1.4 | vdW1 | 0.1564 | - | - | - | 7.13 |
| | vdW2 | 0.0294 | −0.0346 | - | - | 6.65 |
| | WS + NRTL | 1.1384 | - | 0.6327 | 0.2915 | 6.66 |
| 343.15 and 0.1–1.4 | vdW1 | 0.1494 | - | - | - | 7.40 |
| | vdW2 | −0.2053 | −0.101 | - | - | 4.44 |
| | WS + NRTL | −0.1133 | - | 0.4517 | 1.8268 | 4.39 |
C2H6 (1) + [PMIM][Tf2N] (2) | 303.15 and 0.1–1.4 | vdW1 | 0.1240 | - | - | - | 5.23 |
| | vdW2 | 0.236 | 0.0308 | - | - | 1.47 |
| | WS + NRTL | 1.6977 | - | −0.0402 | −0.0385 | 1.53 |
| 323.15 and 0.1–1.4 | vdW1 | 0.0987 | - | - | - | 4.41 |
| | vdW2 | 0.2757 | 0.0467 | - | - | 0.66 |
| | WS + NRTL | 1.7468 | - | −0.5135 | 0.1939 | 0.66 |
| 343.15 and 0.1–1.4 | vdW1 | 0.1028 | - | - | - | 1.93 |
| | vdW2 | 0.0678 | −0.0093 | - | - | 1.86 |
| | WS + NRTL | 1.1086 | - | 0.0247 | 0.3612 | 1.87 |
Table 5.
Henry’s Law constants, enthalpies, and entropies for the solvation of C2H6 in ionic liquids at infinite dilution.
Table 5.
Henry’s Law constants, enthalpies, and entropies for the solvation of C2H6 in ionic liquids at infinite dilution.
Ionic Liquid | (MPa) | (KJ/kmol) | (KJ/mol·K) |
---|
| T = 303.15 K | T = 323.15 K | T = 343.15 K |
---|
[HMIM][Tf2N] | 7.63 | 9.43 | 10.69 | −1.16 | −3.61 |
[BMIM][DMP] | 23.60 | 22.07 | 25.93 | −0.81 | −2.63 |
[PMIM][Tf2N] | 12.85 | 13.61 | 17.35 | −6.40 | −20.02 |
Table 6.
CO2/C2H6 selectivity of the ionic liquids used in this work.
Table 6.
CO2/C2H6 selectivity of the ionic liquids used in this work.
| Temperature (K) |
---|
| 303.15 | 323.15 | 343.15 |
---|
Solvents | | | |
---|
[HMIM][Tf2N] | 2.47 | 2.32 | 2.07 |
[BMIM][DMP] | 3.19 | 2.30 | 2.01 |
[PMIM][Tf2N] | 3.87 | 2.98 | 2.79 |
Selexol | 1.71 | 1.77 | 1.88 |
Table 7.
List of all the solvents and solutes used in this study, structures, acronyms, suppliers, purities, and water content.
Table 7.
List of all the solvents and solutes used in this study, structures, acronyms, suppliers, purities, and water content.
Compounds | Structure | Abbreviation | Suppliers and Purities | Water Content |
---|
Ethane | | C2H6 | Linde (Regina, SK, Canada) (99.99%) | - |
1-Hexyl-3-methylimidazolium bis(trifluormethylsulfonyl)imide | | [HMIM][Tf2N] | Sigma-Aldrich (Oakville, ON, Canada) (≥98.0%) | - |
1-Butyl-3-methylimidazolium dimethylphosphate | | [BMIM][DMP] | IoLiTec (Heilbronn, Germany) (>98.0%) | <2500 ppm |
1-Propyl-3-methylimidazolium Bis(trifluoromethylsulfonyl)imide | | [PMIM][Tf2N] | Sigma-Aldrich (Oakville, ON, Canada) (≥99.0%) | ≤0.5 wt% |
1-Butyl-3-methylimidazolium hexafluorophosphate | | [BMIM][PF6] | Sigma-Aldrich (Oakville, ON, Canada) (≥97.0%) | - |
Table 8.
Molecular weights and critical properties of solutes and solvents used in this study.
Table 8.
Molecular weights and critical properties of solutes and solvents used in this study.
Components | MW (g/mol) | (K) | (MPa) | |
---|
[HMIM][Tf2N] | 447.4 | 1292.8 | 2.39 | 0.3893 |
[BMIM][DMP] | 264.26 | 851 | 2.15 | 0.961 |
[PMIM][Tf2N] | 405.3 | 1259.3 | 3.0 | 0.2575 |
C2H6 | 30.07 | 305.4 | 48.8 | 0.099 |
| Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |