Next Article in Journal
The Effect of the Extraction Conditions on the Antioxidant Activity and Bioactive Compounds Content in Ethanolic Extracts of Scutellaria baicalensis Root
Next Article in Special Issue
Exploration of the Solubility Hyperspace of Selected Active Pharmaceutical Ingredients in Choline- and Betaine-Based Deep Eutectic Solvents: Machine Learning Modeling and Experimental Validation
Previous Article in Journal
Preparation and Thermal Properties of Magnetic PW@CaCO3@Fe3O4 Phase-Change Microcapsules and Their Application to Textile Fabrics
Previous Article in Special Issue
Tuning Ferulic Acid Solubility in Choline-Chloride- and Betaine-Based Deep Eutectic Solvents: Experimental Determination and Machine Learning Modeling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Selectivities of Carbon Dioxide over Ethane in Three Methylimidazolium-Based Ionic Liquids: Experimental Data and Modeling

by
Nadir Henni
,
Amr Henni
* and
Hussameldin Ibrahim
*
Clean Energy Technologies Research Institute (CETRi), University of Regina, Regina, SK S4S 0A2, Canada
*
Authors to whom correspondence should be addressed.
Submission received: 30 July 2024 / Revised: 27 August 2024 / Accepted: 29 August 2024 / Published: 1 September 2024

Abstract

:
This work focused on the solubility of ethane in three promising ionic liquids {1-Hexyl-3-methylimidazolium bis(trifluormethylsulfonyl) imide [HMIM][Tf2N], 1-Butyl-3-methyl-imidazolium dimethyl-phosphate [BMIM][DMP], and 1-Propyl-3-methylimidazolium bis(trifluoromethyl-sulfonyl)-imide [PMIM][Tf2N]}. The solubilities were measured at 303.15 K to 343.15 K and pressures up to 1.4 MPa using a gravimetric microbalance. The overall ranking of ethane solubility in the ionic liquids from highest to lowest is the following: [HMIM][Tf2N] > [PMIM][Tf2N] > [BMIM][DMP]. The Peng–Robinson equation of state was used to model the experimental data using three different mixing rules: van der Waals one, van der Waals two, and Wong–Sandler mixing rules combined with the Non-Random Two-Liquid model. The average absolute deviations for the three mixing rules for the ionic liquids at the three temperatures were 4.39, 2.45, and 2.45%, respectively. Henry’s Law constants for ethane in [BMIM] [DMP] were the highest (lowest solubility) amongst other ionic liquids studied in this work. The solubility ranking for the 3 ILs was confirmed by calculating their overall polarity parameter (N) using COSMO-RS. The selectivity of CO2 over C2H6 was estimated at three temperatures, and the overall ranking of the selectivity was in the following order: [PMIM][Tf2N] > [BMIM][DMP] > [HMIM][Tf2N] > Selexol. Selexol is an efficient and widely used physical solvent in gas sweetening. It has lower selectivity than the three ionic liquids studied. [PMIM][Tf2N], a promising solvent, has the highest selectivity among the three ILs studied and would, therefore, be the best choice if, in addition to carbon dioxide capture, ethane co-absorption was to be avoided. The enthalpy and entropy of solvation at infinite dilution were also estimated.

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 CO2 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 CO2 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 CO2 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 CO2 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 CO2 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 CO2 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

C2H6 solubility data were studied by many researchers who found the ionic liquids to have lower solubility than CO2. This work used three ionic liquids with different cations and anions to understand the factors affecting the solubility of C2H6 in ILs. C2H6 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 C2H6 in ionic liquids. As expected, it is observed that as the temperature rises, C2H6 solubility in the ionic liquid decreases. As pressure rises, C2H6 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 CO2, CH4, C2H6, and others. Anthony et al. [6] studied the anion effect on CO2 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 C2H6 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 C2H6 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 C2H6 solubility than other cation-based ionic liquids. The solubility of C2H6 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 CO2 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.
E r r = 100 n i = 1 n P E x p , i P C a l , i P E x p , i
The average absolute deviation 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 results for C2H6 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 C2H6 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 , H i , 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; f i L and f i V are the fugacity of solute (i) in the liquid phase and vapor phase, respectively; and x i and y i 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:
H i = lim x i 0 f i L T , P , x i x i = lim x i 0 f i V T , P , y i x i
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:
H = R ln H 1 / T P
S = R ln H ln T P
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 C2H6 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 CO2, and low energy requirement. Not much research has been conducted when it comes to the selectivity between gases. When CO2 is captured, it is not the only gas being absorbed. Other gases are absorbed, such as N2O, O2, SO2, H2S, CH4, 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 CO2 and C2H6. The equation used to estimate the selectivity is the following:
S C O 2 / C 2 H 6 = H C 2 H 6 H C O 2
The values of the selectivities are presented in Table 6. CO2/C2H6 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 CO2 was taken from previous research we performed measuring the CO2 solubility in [HMIM][Tf2N], [BMIM][DMP], and [PMIM][Tf2N] [13]. The study revealed that the three solvents were among the best in capturing CO2. The three ionic liquids followed this rank in terms of increasing CO2 solubility: [HMIM] [Tf2N] > [PMIM][Tf2N] >[BMIM][DMP], as reported in reference [13]. The study confirmed that [HMIM][Tf2N] is an excellent solvent for CO2 capture. [PMIM][Tf2N] was also found to have a high capacity for CO2 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 CO2/C2H6 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 C2H6 in ILs obtained in this work were correlated using the Peng–Robinson equation of state (PR-EoS), shown in Equation (6) [15], where a m is the intermolecular attractive force and b m 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).
P = R T v b m a m ( T ) v v + b m + b m ( v b m )

4.3.1. Van der Waals Mixing Rule

For this case, van der Waals created two different mixing rules to estimate mixture parameters ( a m and b m ): 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 am 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), b i j can be estimated with Equation (11).
a m = i j x i x j a i j
b m = i x i b i
b m = i j x i x j b i j
where
a i j = a i i a j j ( 1 k i j )
b i j = b i + b j 2 ( 1 l i j )
where a i i and a j j are the pure component intermolecular attractive force parameter and b i and b j are the pure components co-volume factor parameters for the components ( i and j ).

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 ( k i j , τ i j , and τ j i ), where τ i j and τ j i are the NRTL parameters. The value of α i j in Equations (19) and (20) were set at 0.3 arbitrarily in this work. The energies g i j and g j j are the energies of interaction between molecules ‘ i ’ and ‘ j ’ [17,18,19] and are calculated using Equations (21) and (22).
a = R T Q D 1 D
b = Q 1 D
where
Q = i j x i x j b a R T i j
D = i x i a i b i R T + G e x C R T
C = ln 1 + 2 / 2
where b a R T i j in Equation (14) calculated using Equation (17):
b a R T i j = 1 2 b i a i R T + b j a j R T 1 k i j
  G e x R T = x i x j τ j i G j i x i + x j G j i + τ i j G i j x j + x i G i j
G i j = e x p α i j τ i j
G j i = e x p α i j τ j i
τ i j = g i j g j j / R T
τ j i = g j i g i i / R T
ln γ i = x j 2 τ j i G j i x i + x j G j i 2 + τ i j G i j x j + x i G i j 2
ln γ j = x i 2 τ i j G i j x j + x i G i j 2 + τ j i G j i x i + x j G j i 2

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 C2H6.

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 C2H6 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 C2H6 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

AADAverage 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
CO2Carbon dioxide
C2H6Ethane
DCNDicyanamide
EoSEquation of state
[FSI]Bis(fluorosulfonyl)imide
HHenry’s Law constant
[HMIM][Tf2N]1-Hexyl-3-methylimidazolium bis(trifluormethylsulfonyl)imide
IGAIntelligent Gravimetric Analyzer
ILsIonic liquids
KKelvin
MDEAMethyldiethanolamine
MEAMonoethanolamine
PcCritical pressure
PF6Hexafluoro-phosphate
[PMIM][Tf2N]1-Propyl-3-methylimidazolium bis(trifluoromethyl-sulfonyl)-imide
PR Peng–Robinson
RSO4Alkylsulfate
TbBoiling point temperature
TcCritical temperature
vdW1 van der Waals one
vdW2van der Waals two
WS-NRTLWong–Sandler mixing rules + Non-Random Two-Liquid model
xLiquid mole fraction (moles gas/moles gas + moles solVent)
yVapor mole fraction
ZcCritical compressibility factor

References

  1. 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).
  2. 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).
  3. 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).
  4. 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).
  5. 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]
  6. 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]
  7. 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]
  8. 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]
  9. 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]
  10. 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]
  11. 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]
  12. 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]
  13. 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]
  14. 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]
  15. Peng, D.Y.; Robinson, D.B. A New Two-Constant Equation of State. Ind. Eng. Chem. Fundam. 1976, 15, 59–64. [Google Scholar] [CrossRef]
  16. Renon, H.; Prausnitz, J.M. Local compositions in thermodynamic excess functions for liquid mixtures. AIChE J. 1968, 14, 135–144. [Google Scholar] [CrossRef]
  17. 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]
  18. 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]
  19. 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]
  20. 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 C2H6 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 C2H6 in [BMIM][PF6] and those reported by Anthony et al. (2005) [6] at 298.15 K.
Molecules 29 04152 g001
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.
Molecules 29 04152 g002
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.
Molecules 29 04152 g003aMolecules 29 04152 g003b
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.
Molecules 29 04152 g004aMolecules 29 04152 g004b
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.
Molecules 29 04152 g005
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.
Molecules 29 04152 g006
Figure 7. Comparison of Henry’s Law constants for C2H6 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 C2H6 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].
Molecules 29 04152 g007
Figure 8. Comparison of CO2/C2H6 selectivity in several ionic liquids reported by Nath et al. [12] at 303.15 K.
Figure 8. Comparison of CO2/C2H6 selectivity in several ionic liquids reported by Nath et al. [12] at 303.15 K.
Molecules 29 04152 g008
Figure 9. Comparison of CO2/C2H6 selectivity in several ionic liquids reported by Nath et al. [12] at 323.15 K.
Figure 9. Comparison of CO2/C2H6 selectivity in several ionic liquids reported by Nath et al. [12] at 323.15 K.
Molecules 29 04152 g009
Figure 10. Comparison of CO2/C2H6 selectivity with other ionic liquids reported by Nath et al. [12] at 343.15 K.
Figure 10. Comparison of CO2/C2H6 selectivity with other ionic liquids reported by Nath et al. [12] at 343.15 K.
Molecules 29 04152 g010
Figure 11. IGA-003 setup.
Figure 11. IGA-003 setup.
Molecules 29 04152 g011
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 K323.15 K343.15 K
x C 2 H 6 Pressure (MPa) x C 2 H 6 Pressure (MPa) x C 2 H 6 Pressure (MPa)
0.0030.019910.0030.01990.0030.0199
0.0060.04990.0050.04990.0050.0499
0.0130.099890.0100.10000.0080.0999
0.0250.200040.0210.19980.0200.2002
0.0490.39980.0410.39970.0370.3998
0.0710.599880.0590.59980.0540.5998
0.0920.79990.0770.79980.0700.7995
0.1110.99940.0940.99960.0810.9999
0.1201.09980.1021.09980.0911.1000
0.1291.19880.1101.20110.1031.1994
0.1351.24950.1141.25030.1071.2513
0.1391.30070.1191.3009--
Standard uncertainty u(x) = 0.006; standard uncertainty u(T) = 0.1 K; standard uncertainty u(P) = 0.0008 MPa.
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 K323.15 K343.15 K
x C 2 H 6 Pressure (MPa) x C 2 H 6 Pressure (MPa) x C 2 H 6 Pressure (MPa)
0.0010.019780.0010.01980.0020.0196
0.0020.04990.0020.04990.0020.0499
0.0040.099890.0040.09990.0040.1001
0.0080.199850.0100.19990.0100.1999
0.0170.399910.0180.40020.0170.3996
0.0260.599760.0280.60060.0230.5999
0.0340.799910.0380.79980.0320.7999
0.04220.99980.0450.99980.0431.0006
0.0501.20030.0541.09970.0471.0998
0.0561.29960.0591.20050.0531.2012
0.0591.3486--0.0591.3011
0.0631.4006--0.0651.4006
Standard uncertainty u(x) = 0.006; standard uncertainty u(T) = 0.1 K; standard uncertainty u(P) = 0.0008 MPa.
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 K323.15 K343.15 K
x C 2 H 6 Pressure (MPa) x C 2 H 6 Pressure (MPa) x C 2 H 6 Pressure (MPa)
0.0020.01980.0030.01990.0020.0198
0.0050.04990.0050.04990.0020.0500
0.0090.09990.0080.09990.0050.0999
0.0160.19990.0150.19990.0110.1999
0.0320.40010.0290.39990.0230.4002
0.0450.59990.0420.59970.0330.6011
0.0580.79990.0530.79980.0440.8000
0.0700.99990.0650.99980.0540.9999
0.0831.1991--0.0651.2012
0.0881.30132--0.0701.2998
----0.0791.4000
Standard uncertainty u(x) = 0.006; standard uncertainty u(T) = 0.1 K; standard uncertainty u(P) = 0.0008 MPa.
Table 4. Optimized binary interaction parameters for different models and calculated AAD%.
Table 4. Optimized binary interaction parameters for different models and calculated AAD%.
CompoundTemperature (K)
and Pressure (MPa)
Mixing Rule k 12  a l 12  b τ 12 τ 21 AAD%
C2H6(1) + [HMIM][Tf2N] (2)303.15 and 0.1–1.4vdW10.0797---3.19
vdW20.15660.0163--0.47
WS + NRTL1.305-0.2151−0.18480.44
323.15 and 0.1–1.4vdW10.0694---2.11
vdW20.11160.0104--0.83
WS + NRTL1.2293-0.7031−0.58230.84
343.15 and 0.1–1.4vdW10.0452---3.37
vdW20.05570.0026--3.30
WS + NRTL1.0256 0.03960.01233.30
C2H6 (1) + [BMIM][DMP] (2)303.15 and 0.1–1.4vdW10.1931---4.77
vdW20.052−0.0377 -2.35
WS + NRTL0.9997-1.10240.63622.34
323.15 and 0.1–1.4vdW10.1564---7.13
vdW20.0294−0.0346--6.65
WS + NRTL1.1384-0.63270.29156.66
343.15 and 0.1–1.4vdW10.1494---7.40
vdW2−0.2053−0.101--4.44
WS + NRTL−0.1133-0.45171.82684.39
C2H6 (1) + [PMIM][Tf2N] (2)303.15 and 0.1–1.4vdW10.1240---5.23
vdW20.2360.0308--1.47
WS + NRTL1.6977-−0.0402−0.03851.53
323.15 and 0.1–1.4vdW10.0987---4.41
vdW20.27570.0467--0.66
WS + NRTL1.7468-−0.51350.19390.66
343.15 and 0.1–1.4vdW10.1028---1.93
vdW20.0678−0.0093--1.86
WS + NRTL1.1086-0.02470.36121.87
a  k 12 = k 21 ;   b  l 12 = l 21 .
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 H e n r y s   L a w   C o n s t a n t (MPa) H
(KJ/kmol)
S
(KJ/mol·K)
T = 303.15 KT = 323.15 KT = 343.15 K
[HMIM][Tf2N]7.639.4310.69−1.16−3.61
[BMIM][DMP]23.6022.0725.93−0.81−2.63
[PMIM][Tf2N]12.8513.6117.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.15323.15343.15
Solvents
[HMIM][Tf2N]2.472.322.07
[BMIM][DMP]3.192.302.01
[PMIM][Tf2N]3.872.982.79
Selexol1.711.771.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.
CompoundsStructureAbbreviationSuppliers and PuritiesWater Content
EthaneMolecules 29 04152 i001C2H6Linde (Regina, SK, Canada) (99.99%)-
1-Hexyl-3-methylimidazolium bis(trifluormethylsulfonyl)imideMolecules 29 04152 i002[HMIM][Tf2N]Sigma-Aldrich (Oakville, ON, Canada) (≥98.0%)-
1-Butyl-3-methylimidazolium dimethylphosphateMolecules 29 04152 i003[BMIM][DMP]IoLiTec (Heilbronn, Germany) (>98.0%)<2500 ppm
1-Propyl-3-methylimidazolium Bis(trifluoromethylsulfonyl)imideMolecules 29 04152 i004[PMIM][Tf2N]Sigma-Aldrich (Oakville, ON, Canada) (≥99.0%)≤0.5 wt%
1-Butyl-3-methylimidazolium hexafluorophosphateMolecules 29 04152 i005[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.
ComponentsMW (g/mol) T c (K) P c (MPa) ω
[HMIM][Tf2N]447.41292.82.390.3893
[BMIM][DMP]264.26 8512.150.961
[PMIM][Tf2N]405.31259.33.00.2575
C2H630.07305.448.80.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.

Share and Cite

MDPI and ACS Style

Henni, N.; Henni, A.; Ibrahim, H. Selectivities of Carbon Dioxide over Ethane in Three Methylimidazolium-Based Ionic Liquids: Experimental Data and Modeling. Molecules 2024, 29, 4152. https://doi.org/10.3390/molecules29174152

AMA Style

Henni N, Henni A, Ibrahim H. Selectivities of Carbon Dioxide over Ethane in Three Methylimidazolium-Based Ionic Liquids: Experimental Data and Modeling. Molecules. 2024; 29(17):4152. https://doi.org/10.3390/molecules29174152

Chicago/Turabian Style

Henni, Nadir, Amr Henni, and Hussameldin Ibrahim. 2024. "Selectivities of Carbon Dioxide over Ethane in Three Methylimidazolium-Based Ionic Liquids: Experimental Data and Modeling" Molecules 29, no. 17: 4152. https://doi.org/10.3390/molecules29174152

Article Metrics

Back to TopTop