Open Access Article
Rohan Sartape
a,
Rashmi Mishrab,
Vamsi V. Gandea,
Sophia Johnsona,
Bianca Berryc,
Ishaku Amosa,
Rohit Chauhana,
Amey Thorat
b,
Jindal K. Shah
b and
Meenesh R. Singh
*a
aDepartment of Chemical Engineering, University of Illinois Chicago, 929 W. Taylor St., Chicago, 60607, Illinois, USA. E-mail: mrsingh@uic.edu; Tel: +1 312-996-3424
bSchool of Chemical Engineering, Oklahoma State University, Stillwater, 74078, Oklahoma, USA
cLyons Township High School, Western Springs, 60558, IL, USA
First published on 27th March 2026
CO2 hydroxylation is the underlying reaction for a wide range of reactive CO2 capture (RCC) techniques. Despite being thermodynamically favorable, CO2 hydroxylation kinetics are sluggish due to a higher energy barrier associated with bending the linear CO2 molecule before nucleophilic attack by OH−. Here, we report a previously unrecognized catalytic effect in which ionic liquids (ILs) increase the rates of CO2 hydroxylation by pre-activating (or bending) CO2 locally around the IL. ILs are known for their higher CO2 solubilities due to their larger void fractions; however, their role in catalyzing the CO2 hydroxylation reaction, an important step in RCC, has not been reported previously. NMR and FTIR measurements, together with quantum-chemical calculations, are consistent with bicarbonate being the dominant stable species, while carbonate appears to revert to bicarbonate over time under the studied conditions. For instance, the addition of BMIM NTf2 to a KOH-ethylene glycol mixture accelerates CO2 hydroxylation by sixfold by lowering the activation energy (∼33%) without chemically binding CO2, highlighting its catalytic role. Isotope-labeled ATR-FTIR experiments confirm enhanced CO2 bending in ionic liquids through the appearance of a hot band, indicative of increased population of thermally accessible bent vibrational states. While low IL loadings enhance kinetics, higher concentrations hinder CO2 mass transfer due to reduced interfacial tension. The nature of IL cation and anion strongly influences rates, with shorter alkyl chains and I−/OTf− anions providing the highest activity. The system demonstrates reversibility and robustness under anhydrous flue gas conditions, offering a tunable pathway for efficient CO2 capture.
Broader contextThe escalating need to mitigate anthropogenic CO2 emissions has driven intense efforts to develop capture technologies that are both energy-efficient and scalable. Conventional amine-based processes, while widely deployed, suffer from high energy requirements for regeneration, volatility, and corrosivity, underscoring the need for sustainable alternatives. Ionic liquids (ILs) and deep eutectic solvents (DESs) offer promise due to their tunable properties, low volatility, and strong CO2 solubility, yet practical limitations such as high viscosity and slow mass transfer have hindered large-scale adoption. This study demonstrates a new concept: ILs can act as catalysts for CO2 hydroxylation, promoting bicarbonate formation through molecular pre-activation rather than direct binding. Using NMR, FTIR, and quantum-chemical simulations, we show that ILs bend and activate CO2, significantly lowering the activation energy for hydroxylation while maintaining reversibility. The catalytic effect is strongly influenced by IL structure, with cation alkyl chain length and anion chemistry dictating efficiency and stability, thereby offering rational design principles for optimized capture systems. By revealing IL-catalyzed hydroxylation as a tunable and scalable mechanism, this work provides critical molecular-level insights and positions IL-based systems as viable, energy-efficient alternatives to amines. The findings establish a foundation for advancing next-generation solvent systems that can overcome current bottlenecks in sustainable carbon capture. |
Ionic liquids (ILs) offer high CO2 capacity, tunable chemistry, and negligible vapor pressure.25–27 However, their high viscosity can severely limit CO2 mass transfer13 and challenges in synthesis and scalability further constrain their practical deployment.28–31 Although aprotic heterocyclic anions (AHAs) are available that can react with CO2, they are often described as “lower-viscosity” relative to many CO2-reactive ILs (∼3–20 mPa s for AHAs,32 ∼300–400 mPa s for acetate-based ILs at 313 K). They are not yet commercially available, making a comprehensive techno-economic analysis difficult, if not impossible.33–35 Furthermore, the viscosities of AHAs remain slightly higher than those of industrial aqueous amines, such as 30 wt% MEA (∼2–3 cP at 313 K), which can impose additional mass-transfer and pumping penalties. Deep eutectic solvents, an alternative to ILs for CO2 capture, are formed by combining a hydrogen bond acceptor (HBA) and a hydrogen bond donor (HBD), resulting in a mixture with a melting point lower than that of the individual components.36–38 At a certain molar ratio (typically 1
:
1, 1
:
2 or 2
:
1) of HBA to HBD, the melting point of the DES can be reduced significantly lower than that at the eutectic composition, enabling their applications at room temperature. DES inherits almost all the advantages of IL but offers lower production costs and ease of synthesis.39–41 Alkaline DES, composed of an alkaline salt as the HBA (e.g., potassium hydroxide) and an HBD (e.g., ethylene glycol), has demonstrated excellent CO2 capture performance.42–44
Recent studies from our group demonstrated that alkaline DES can achieve CO2 capture rates as high as 0.13 mmol m−2 s−1 with regeneration energies of 120 kJ mol−1 of captured CO2, using a pH-swing process.45–49 Mishra et al.50 investigated various intermolecular interactions involved in alkaline DES, comprised of ethylene glycol (EG) and KOH, to understand its enhanced performance in CO2 capture. Using FTIR and quantum chemical analyses, noncovalent interactions and structural changes were found to be critical for optimizing DES design for direct-air-capture applications. In parallel, mixing DESs with ILs (DES–IL) has been proposed to reduce energy losses in electrochemical CO2 capture.51,52 Our group observed a non-monotonic ionic conductivity behavior in various EG–IL mixtures as a function of concentration, attributed to ion pair aggregation that limits charge mobility53,54 and was linked to the extent of ion fluctuations in the first solvation shell of a given ion.55 Although the physicochemical properties of alkaline DES–IL mixtures have begun to be investigated, CO2 chemisorption kinetics in DES–IL mixtures remain unexplored and present an important avenue for future research.
In general, hydroxylation reactions with CO2 require prior CO2 activation and high activation energy, due to its linear, symmetric structure and chemical inertness. Increasing the temperature can facilitate activation but is often energy-intensive, making the use of activation agents/catalysts such as N-heterocyclic carbenes, polyoxometalates, and ILs essential.56,57 ILs, through electrostatic interactions, can interact with CO2.51,58–63 It is well established that the solvation effect and solvation energy in CO2 hydroxylation play a crucial role.64,65 Iida et al. demonstrated that CO2 hydroxylation in the gas phase is barrierless. However, in the aqueous phase, the reaction exhibits a significant activation barrier of 54 kJ mol−1. Their work further revealed that CO2 hydroxylation involves the desolvation of the hydroxide moiety, and the oxygen atom of CO2 increasingly attracts solvent molecules as the reaction coordinate decreases.64 In a similar vein, the presence of IL can modulate the solution energies around the hydroxide moiety and the O atom of CO2. This altered energy landscape facilitates CO2 bending and lowers the effective energy barrier, thereby catalyzing the CO2 capture process. A schematic illustration of the phenomena is also presented in Fig. 1. In an aqueous medium, the formation of bicarbonate or carbonate products from the hydroxylation of CO2 is strongly governed by how efficiently water can transfer protons and stabilize transition states.66 Since ILs and DES usually have trace amounts of water, products can remain in bicarbonate form, which is advantageous.67,68 Such a concept is also widely recognized in catalysis, where modulation of solvation energy plays a crucial role in lowering activation barriers. Various areas of chemistry leverage solvation-energy tuning to enhance reaction rates, including enzymatic catalysis69,70 and confined nanocage catalysis.71,72 A notable biological example is the enzyme carbonic anhydrase (CA), which catalyzes the hydration of CO2 in mammalian systems. The catalytic cycle of CA involves the deprotonation of a coordinated water molecule to generate a nucleophilic hydroxide species that subsequently attacks CO2 to form bicarbonate.73 Inspired by this mechanism, CA has also been explored in combination with amine solutions to enhance CO2 capture rates.74,75 However, the large-scale deployment of enzyme-assisted systems remains challenging due to the high cost and stability limitations associated with enzyme production and operation.76
Building on these principles, DESs have been explored as potential catalysts or activation agents for various reactions,77,78 offering sustainability and tunability.79–81 Similarly, IL–solvent mixtures, particularly at low ionic liquid mole fractions (xIL≤ 0.05), exhibit enhanced chemical activity compared to concentrated IL systems.82–92 For example, Man et al.82 reported maxima in substitution reaction rates (xIL = 0.02) for 1-butyl-3-methylimidazolium bis(trifluorosulfonyl)imide (BMIM NTf2), attributing this to microscopic interactions that influence the transition state. Zhang et al.93 demonstrated that free BMIM+ cations impact the hydroxylation rate of 2-methylcyclohexanone via ionic mechanisms. Similarly, Gilbert et al.89 reported enhanced nucleophilic substitution rates using small amounts of BMIM NTf2 without affecting stereoselectivity. Kochly et al.92 further highlighted that IL-induced rate enhancements in elimination reactions depend on the hydrogen-bonding properties of the solvent and the steric effects of the IL anion.
Herein, we demonstrate that an alkaline DES–IL mixture (specifically, KOH, ethylene glycol (EG), and IL) enhances CO2 capture by facilitating its conversion to a bicarbonate species. Due to the alkalinity of the system, an interplay between bicarbonate and carbonate species occurs, facilitated by in situ water formation during the CO2 capture process. To elucidate the reaction mechanism of IL-catalyzed CO2 capture, we systematically varied the IL concentration in DES–IL mixtures and observed the impact on the CO2 capture performance. Finally, to understand the role of IL structure, we screened various DES–IL combinations to evaluate the influence of both cation and anion chemistry on CO2 uptake. This study underscores the potential of DES–IL systems as a cost-effective, tunable platform for efficient CO2 capture and offers insight for the rational design of next-generation carbon capture technologies.
We also acknowledge that, in the field of CO2 capture, promoter species such as piperazine are widely used to enhance kinetics by increasing the rate of CO2 hydration in systems such as tertiary amines or K2CO3.94–97 In the present work, however, our spectroscopic observations and mechanistic analysis indicate that the enhancement observed includes catalytic contributions rather than solely conventional promoter behavior.
For the system containing DES, bubbling 13CO2 for 2 minutes, results in the appearance of new resonance peaks at 61, 67, 159, and 167 ppm. Furthermore, as time increases, these peaks appear to increase in intensity except for 167 ppm, which shifts upfield and decreases in intensity. Peaks at 159 and 167 ppm, labelled as α and β, can be ascribed to bicarbonate and carbonate species involved in the EG-KOH system. NMR signals from mixing isotope-labelled 0.1 M sodium bicarbonate and sodium carbonate with EG-KOH are also shown in Fig. S1 to validate the assigned peaks. The increase in the bicarbonate peak at 159 ppm over time, along with its stable chemical shift, coupled with the gradual decrease and downfield shift of the carbonate peak at 167 ppm, suggest that both species form at the beginning of the reaction. After 2 minutes of bubbling CO2 through the DES system, the carbonate peak begins to decrease in intensity, while the bicarbonate peak continues to increase. This indicates that carbonate progressively regenerates into bicarbonate. The process is analogous to the well-established pH-driven equilibrium dynamics of CO2 when dissolved in water.98 A tiny peak, labelled as δ, is also observed at 160.9 ppm after 8 minutes of 13CO2 bubbling. This peak is ascribed to the CO2 attached to the deprotonated ethylene glycol, as shown in Fig. 2(a). A similar set of new peaks at 61, 67, 159, and 167 ppm is again observed in the system containing BMIM NTf2 mixed with DES in Fig. 2(b). A comparison of normalized signal intensity at 159 ppm in Fig. S2 reveals an accelerated conversion to bicarbonate in the presence of BMIM NTf2. Furthermore, after 8 minutes of bubbling with 13CO2, a signal at 124 ppm corresponding to dissolved CO2 appears, indicating the completion of the chemisorption process.
In experiments with DES and DES–IL mixtures, resonance peaks at 61 and 67 ppm also increase with time and are ascribed to the methylene groups of the hydroxy ethyl carbonate (HEC) labelled as 2 and 3, as shown in Fig. 2(b).99 Such a characteristic peak indicates the attack of the deprotonated form of EG on CO2.99–101 Cui et al.99 observed that for a CO2 absorption system containing azolide-based ILs, instead of the formation of carbamate species by reacting with azolide ions, the CO2 reacts with the –OH group of ethylene glycol to form the carbonate species. Zhou et al.101 observed a similar reaction product when EG reacted with K2CO3. Based on their observations, the overall equilibrium reaction involves the reaction of carbonate and EG to form a hydroxide ion and HO–CH2–CH2–O–COO−. Based on their analysis, reaction pathways were proposed that involved the formation of bicarbonate, carbonate, water, and carbonic acid. Other published works also indicate that CO2 binds to the EG molecule following its deprotonation.36,102 These studies indicate that CO2 absorption in alkaline DES involves deprotonation of EG and formation of HO–CH2–CH2–O–COO− (HEC), which is also a feasible pathway other than the formation of free bicarbonate and carbonate species.
To gain further insight, the samples were further analyzed through 1H NMR, and the results are presented in Fig. 2(c). After bubbling 13CO2, a new peak at 3.7 ppm is observed in the system with DES. This peak is assigned to the proton associated with EG labelled as 3. A gradual shift in the peak position of –OH peaks at ∼5.1 ppm can also be observed in systems containing DES and DES–IL mixtures as 13CO2 is introduced, as shown in Fig. S3 and S4. An upfield shift in the –OH signal indicates the formation of water as a reaction product. Similar experiments were performed by Restrepo and Adams, who studied the effect of water content in glycols on 1H NMR shift and outlined a technique to estimate water content in the system.103,104
Additional experiments were performed to confirm that CO2 chemisorption occurred even in the presence of anhydrous flue gas (12.2% CO2, 0.16% SO2, 7% O2, and 80.2% N2). Fig. S5 and S6 show the presence of bicarbonate peaks via 13C NMR analysis. Systems with and without BMIM NTf2 in alkaline DES show the formation of bicarbonate. A comparison of normalized peak, at ∼159 ppm, intensities also shows that even in the presence of flue gas, CO2 chemisorption is accelerated in the system containing BMIM NTf2, as shown in Fig. S7. Presence of water in the flue gas can have detrimental impacts on CO2 capture process by reacting with the formed bicarbonate, yielding back to gaseous CO2 and carbonate. Interested readers can look into already published works.48,105 To validate the reversibility of the CO2 chemisorption, diluted H2SO4 was added to the DES. Results in Fig. S8 validate our observations and prove the reversibility of CO2 chemisorption in DES and DES–IL mixture for the pH-swing process.48,106 Typically, for the pH swing process, the carbonate and bicarbonate species remain in equilibrium, and a basic pH preserves bicarbonate and carbonate, while an acidic pH will release CO2 in gaseous form. The NMR shifts at 59 and 65 ppm become insignificant, indicating the reversible proton uptake by deprotonated EG and the dissolution of captured CO2 into gaseous CO2. Moreover, the bicarbonate peak at 159 ppm also disappears, indicating desorption of CO2 from the system.
The DES and DES–IL systems were also simulated using quantum-chemical calculations to validate the assigned peaks from NMR experiments. Details about the simulation are also provided in Section S1. The simulated 13C NMR spectra display with and without peak splitting, as shown in Fig. S9(a–c), with the corresponding chemical shift values provided in Table S2 in comparison to those of EG. This splitting of the EG peak is also observed in the experimental NMR spectra, confirming its deprotonation. Physical absorption of CO2 is generally characterized by 13C NMR signals in the range of 120–130 ppm,107 while chemisorption results in signals mainly between 155–165 ppm. As shown in the simulated 13C NMR spectrum in Figure S10(b), physisorption leads to a new carbon peak at approximately 130 ppm without any splitting of the EG peak. In contrast, chemisorption leads to a carbon peak near 163 ppm in Fig. S10(c), consistent with the experimentally observed peak at 160.9 ppm. The splitting of the EG peak in the chemisorbed state, Fig. S10(c), is assigned to EG deprotonation. The consistency between simulated and experimental data supports the conclusion that CO2 undergoes chemisorption with DES.
Calculated NMR shifts of CO2, bicarbonate, carbonate, and HEC are provided in Fig. S11(a)–(e) with the corresponding chemical shift values provided in Table S4. The simulated NMR spectra of deprotonated EG with CO2, shown in Fig. S11(e), closely match the experimental peaks. This indicates the formation of a CO2 adduct with deprotonated EG. The calculated chemical shifts at 67, 68, and 163 ppm correspond well with the experimental shifts at 61, 67, and 160.9 ppm, confirming chemisorption of CO2 by EG. Additionally, a peak at 167 ppm in the experimental spectra, which later disappears, can be attributed to carbonate. The calculated chemical shift for carbonate appears at 184 ppm in Fig. S11(c), but may be shifted upfield due to the experimental reaction environment. This supports the experimental observation that the peaks at 159, 160.9, and 167 ppm in the 13C NMR spectra correspond to bicarbonate, HEC, and carbonate species, respectively, indicating that both species are present at the onset of the reaction. The formation of a hydrogen peak at 5.9 ppm in the simulated 1H NMR spectrum of water contributed to the combination of the deprotonated hydrogen from EG and the hydroxyl ion from KOH, which correlates well with the 5.1 ppm peak observed in the experimental spectra. This assignment, as illustrated in Fig. S12 with the corresponding chemical shift values provided in Table S5, confirms the presence of water in the reaction system.
Quantum calculations were also carried out for BMIM+ and its interaction with CO2, as shown in Fig. S13(a and b), with the corresponding chemical shift values provided in Table S6. The simulated 13C NMR spectra for the DES–IL mixture in Fig. S13(b) reveal a new peak at around 162 ppm, as well as splitting of the C3 and C6 peaks upon CO2 interaction with BMIM+. However, this splitting is never observed in the experimental spectra, indicating that the BMIM+ cation does not interact with CO2. Theoretical NMR peaks for NTf2− and its interaction with CO2 are also shown in Fig. S14(a–c) with the corresponding chemical shift values provided in Table S7. Similar to BMIM+, if CO2 were to interact with the NTF2− anion, as depicted in Fig. S13(b and c), one would expect peak shifting and the appearance of an additional carbon signal in the 13C NMR spectra. Since no such peak shifting or additional peak is observed experimentally, it can be concluded that the anion does not chemically bond with CO2. Similarly, the absence of new peaks and chemical shifts corresponding to cation–anion interactions with CO2 in the experimental spectra suggests a plausible mechanism: BMIM NTF2 can function as a catalyst in the absorption of CO2 in the alkaline DES–IL system.
While promoter species are known to enhance kinetics in amine-based CO2 capture systems, the literature also recognizes that promoters may participate directly in CO2-binding equilibria, and promoter–CO2 reaction products can be present under capture conditions.108 In contrast, post-absorption NMR spectra of our system show no new peaks attributable to irreversible IL–CO2 adduct formation or IL degradation, indicating that the ionic liquid remains chemically conserved under the studied conditions. Taken together, these observations are consistent with an IL-enabled kinetic enhancement that occurs without measurable consumption of the ionic liquid, suggesting that the system exhibits catalytic contributions rather than solely conventional promoter behavior.
O stretching band of bicarbonate and carbonate species.4,109,110 The shoulder peak at 1300 cm−1 is ascribed to the stretching band of the C–O stretching vibration of HO–CH2–CH2–O–COO−. A similar band assignment was also made at 1296 cm−1 by Cheng et al. where CO2 absorption took place in a deep eutectic solvent comprising EG.111 Fig. 2(e) represents the FTIR spectra of 10% CO2 (vol%, balance N2) bubbling in the presence of BMIM NTf2 in alkaline DES. Three new peaks are observed at 1639, 1529, and 1300 cm−1, which are similar to the FTIR results obtained for the system without IL. These results also align with our previous NMR spectroscopy observations. Thus, corroborating our assumption that IL does not participate or form a chemical bond with CO2 in the reaction and acts like a catalyst for CO2 absorption in our system. Full spectra from the FTIR experiment with and without IL in alkaline DES have also been provided in Fig. S15 and S16. Additional information on experimental methods of FTIR is provided in Section S1.The interaction between CO2 and EG was further examined through simulated IR spectroscopy, as shown in Fig. S17. Comparison of the simulated IR spectra of hydroxyethyl carbonate (HEC) and EG reveals the appearance of new peaks at 1262, 1474, and 1649 cm−1 in the HEC spectrum. The peak at 1262 cm−1 in HEC closely corresponds to bands in the simulated spectra of bicarbonate (HCO3−) at 1278cm−1 and carbonate (CO32−) at 1311 cm−1 and is assigned to C
O and C–O stretching vibrations. The band at 1649 cm−1 in HEC matches the corresponding band in HCO3− and is assigned to the asymmetric stretching vibration of CO2. The peak at 1474 cm−1 is assigned to the symmetric stretching vibration of C
O. This band is absent in the simulated spectra of EG, HCO3−, and CO32−, where the molecular symmetry remains preserved. The assigned peaks in the simulated IR spectra agree well with the experimental bands at 1300, 1558, and 1638 cm−1. These results are consistent with the FTIR experimental observations, confirming the chemical interaction of CO2 with EG and the formation of bicarbonate and carbonate species.
Since an increased activity for CO2 chemisorption to bicarbonate was observed, experiments were conducted at higher temperatures to determine the activation energy (Ea) and to understand the applicability in a large-scale process.112–114 Experiments were conducted at different temperatures in a similar setup modified with a heating plate at the bottom of the glass vial to obtain Ea. Fig. 3(b) shows the comparison of the increase in reaction rates with temperature for DES without IL and with IL. Assuming a pseudo 1st order reaction is involved for the conversion of hydroxide to bicarbonate, the rate constants were determined to be 0.29, 0.43, 0.84, and 1.53 min−1 at 298 K, 313 K, 323 K, and 333 K, for the system without IL, respectively, as shown in Fig. 3(b) and Table S8. Furthermore, for the system containing BMIM NTf2, the rate constants were determined to be 0.73, 1.06, 1.48, and 2.3 min−1 at 298 K, 313 K, 323 K, and 333 K. Additional details regarding the pseudo 1st order reaction calculations are provided in Section S1. Timed data for CO2 hydroxylation at different temperatures are provided in Fig. S21. Furthermore, while the catalytic enhancement appears most pronounced at lower temperatures, this trend cannot be attributed to viscosity effects. A comparison of viscosity values for systems with and without IL shows a difference of only ∼1 cP at room temperature, which decreases further at elevated temperatures. Thus, changes in viscosity with temperature are insufficient to explain the diminished catalytic enhancement at higher temperatures. Instead, CO2 hydroxylation is primarily governed by the activation energy associated with CO2 activation and solvent reorganization. As reported by Iida et al.64 this barrier dominates the reaction kinetics. At higher temperatures, the intrinsic thermal energy increasingly overcomes this activation barrier even in the absence of IL, reducing the relative benefit of IL-induced solvent reorganization, making the catalytic effect less distinguishable. We also acknowledge that CO2 solubility decreases at elevated temperatures, particularly in organic systems. However, at the compositions investigated (5 vol% IL), the measured Henry constants for EG and EG–IL mixtures are comparable115 (∼485 bar vs. ∼415 bar), indicating only a modest increase in CO2 solubility upon IL addition. Such a difference corresponds to a relatively small change in dissolved CO2 concentration and is insufficient to account for the observed enhancement in capture rate. Therefore, while solubility effects may contribute to the overall behavior, they cannot alone account for the magnitude of the kinetic enhancement observed in the DES–IL system.
Fig. 3(c) shows the Arrhenius plot for estimating the Ea of CO2 chemisorption in the EG/KOH and EG/KOH/BMIM NTf2 system. Activation energy (Ea) values for CO2 chemisorption were determined by linear regression, yielding a satisfactory fit. It is evident that systems containing IL possess a reduced Ea for CO2 chemisorption, since the magnitude of the slope of the system containing IL is significantly lower than that of the system without IL. Fig. 3(d) shows a comparison of Ea with and without IL. The Ea required for the CO2 hydroxylation system without IL is 39.3 kJ mol−1, while the system with IL requires 26.6 kJ mol−1, which is a 33% reduction in Ea for CO2 chemisorption. Results from Eyring analysis of the same data shown in Table S9 revealed a similar drastic decrease in activation enthalpy and entropy (ΔH and ΔS), indicating that the rate constant enhancement is primarily governed by the ionic-liquid transition state. Such changes in the Ea were also observed by Hsieh et al.,84 who argued that the presence of ionic liquid can assist in the development of less charge at the transition state. In other words, the presence of ionic liquid modifies the nature of the transition state, which determines the extent of bond formation/breaking involved in bicarbonate/carbonate species.
Furthermore, Wang et al.116 have already argued that ILs reduce the activation barrier for CO2 hydroxylation by stabilizing polarized CO2 and key reaction intermediates through electrostatic interactions and hydrogen bonding. These interactions distort the linear CO2 molecule and facilitate nucleophilic attack, lowering the energetic cost of CO2 activation compared to non-ionic solvent environments. Another review by Wang et al.117 discusses the Arrhenius analysis reported for the synthesis of dimethyl carbonate (DMC) from CO2 and methanol in the presence of [GLY(mim)3][NTf2]3, which demonstrates a clear reduction in the apparent activation energy relative to the IL-free system. This decrease is attributed to strong electrostatic interactions and hydrogen bonding between the ionic liquid and CO2, which stabilize activated CO2 species and key transition states. The ionic liquid also enhances local reactant organization and polarizes CO2, thereby lowering the energetic cost of bond rearrangement. Collectively, these effects accelerate reaction kinetics, confirming the role of ionic liquids as catalysts. To validate, the enhancement in capture rate also occurs in the presence of IL under DAC-relevant conditions. We performed an additional experiment, details of which are provided in Section S1. As shown in Fig. S22, under 500 ppm CO2 feed conditions, the CO2 conversion rate is also enhanced by the presence of IL. The kinetics under 500 ppm CO2 conditions are sluggish due to a decrease in the effective driving force.118,119
Since it is evident that CO2 chemisorption is catalyzed, the next step is to identify the reaction catalyzed by the IL. To investigate this, the instantaneous reaction rates vs. concentration of OH− are plotted in Fig. 3(e). By comparison with the magnitude of the extent of the reaction as shown in Fig. S23, it can be further confirmed that reaction 1, i.e., bicarbonate formation, is catalyzed, as shown in the green region, since the maximum instantaneous reaction rate overlaps with the bicarbonate formation region. Furthermore, ζ1 (extent of reaction 1) becomes dominant after ∼25% conversion of hydroxide has taken place, as shown in Fig. 3(e) and Fig. S23. As previously discussed, the presence of IL lowers the activation energy for bicarbonate formation by bending CO2. Additionally, in Fig. S23 and S18(d), the ζ3 (extent of reaction 3) dominated region is small and negative, indicating that water formation occurs in this region during the initial time steps of CO2 bubbling. This is further validated by the significant NMR upfield shifts observed in OH− peaks at ∼5.1 ppm only at the initial time steps of CO2 bubbling, in Fig. S3 and S4.
Based on the extent-of-reaction analysis, the bicarbonate formation pathway is clearly catalyzed. As discussed in the Introduction, facilitated CO2 bending can lower the effective activation barrier, thereby catalyzing CO2 hydroxylation to bicarbonate.64 To experimentally probe whether enhanced CO2 bending occurs in the ionic liquid (IL) environment, ATR-FTIR spectroscopy was employed. The bending mode of 12CO2 (∼660 cm−1) lies at sufficiently low frequency that, at room temperature (kBT ≈ 200 cm−1), it is thermally populated for a small but non-negligible fraction of CO2 molecules.58,120 To eliminate interference from atmospheric CO2, isotopically labeled 13CO2 was used. As shown in Fig. 3(f), 13CO2 exhibits an asymmetric stretching band at 2277 cm−1. In the presence of BMIM NTf2, the intensity of this band increases with time, accompanied by the emergence of a weak shoulder approximately 11 cm−1 lower in frequency. This shoulder has been previously assigned to a hot band arising from coupling between the thermally populated bending mode and the asymmetric stretch, leading to an apparent red shift of the absorption feature.61,63 In contrast, for EG, no hot band is observed even after 10 minutes of 13CO2 bubbling, and the intensity of the 13CO2 asymmetric stretch at ∼2272 cm−1 remains unchanged. In ionic liquids, strong and spatially heterogeneous local electric field gradients disrupt the symmetric solvation shell around CO2, increasing polarization and enhancing the population of thermally accessible bent vibrational states. This manifests as increased coupling between the bending and asymmetric stretch modes, without requiring direct chemical coordination of CO2 to the ions.63 These observations indicate that enhanced CO2 bending occurs selectively in the presence of the IL, supporting our conclusion that the IL environment facilitates CO2 activation and thereby catalyzes bicarbonate formation.
Several alternative explanations could, in principle, be invoked to rationalize the observed enhancement in CO2 hydroxylation rates in the presence of ionic liquids (ILs). (a) CO2 solubility effects. BMIM NTf2 is indeed known for its high intrinsic CO2 solubility; however, the IL loading employed in this study is limited to 5 vol%, corresponding to a mol fraction of ∼0.05. At this low concentration, recent calculations from our group show that the Henry's constant (kH) decreases only modestly (from ∼500 bar to ∼270 bar), compared to the value for pure IL (∼29 bar).115,121,122 This relatively small change in kH cannot account for the drastic enhancement in hydroxylation rates observed, ruling out solubility-driven effects as the dominant mechanism. (b) Viscosity effects. The viscosity of pure BMIM NTf2 (∼48 cP) is higher than that of ethylene glycol (∼22 cP); however, applying a logarithmic mixing rule yields an estimated viscosity of ∼23 cP at a 0.05 mol fraction of IL, which is an increase of only ∼1 cP. Such a marginal change would be expected to slightly impede, rather than accelerate reaction kinetics, further weakening viscosity as a plausible explanation. (c) Mass-transfer limitations. All experiments were performed under conditions designed to ensure kinetic control. As shown in Fig. S24, increasing the CO2 flow rate beyond 10 mL min−1 does not further enhance hydroxylation rates, indicating saturation of the gas–liquid mass transfer rates. This behavior confirms that the system operates in a kinetically limited regime rather than being governed by external mass-transfer constraints. (d) Interfacial effects. While IL cations with long alkyl chains are known to be surface-active and may accumulate at the gas–liquid interface, this interfacial structuring does not explain the observed rate enhancement. In fact, as demonstrated in the subsequent section, increasing the IL concentration beyond an optimal value leads to slower CO2 hydroxylation, consistent with interfacial saturation and mass-transfer hindrance rather than acceleration. Taken together, these analyses demonstrate that changes in CO2 solubility, viscosity, mass transfer, or interfacial effects cannot rationalize the experimentally observed rate enhancement. Thus, the only explanation consistent with all observations is that the IL plays a catalytic role in facilitating CO2 hydroxylation.
To elucidate the impact of the interfacial area of the fluid on CO2 chemisorption phenomena, a mathematical model based on mass transfer in a bubble column was developed to describe the underlying reaction mechanisms and changes in interfacial property.129–131 Inspired by the Michaelis–Menten kinetics,132,133 the mathematical model approximates the reaction of IL with CO2 as a pseudo-steady state reaction, as shown in eqn (1), where the IL acts as a catalytic site that transiently binds CO2 to form an intermediate complex. After ILCO2 is formed, the reaction proceeds with the hydroxide dissolved in the liquid as shown in eqn (2). The rate law model, combined with the bubble column mass transfer model, was then fit to a system containing different volume percentages of IL, using experimental data on CO2 chemisorption with a DES–IL mixture to gain insight into the model parameters. Such an approach captures the saturation behavior observed at higher IL concentrations, where mass transfer or subsequent reaction steps, rather than CO2 binding, become rate-limiting. Additional details about the mathematical model, reactions, and assumptions are provided in Section S3.
![]() | (1) |
![]() | (2) |
Results from model fitting are in satisfactory agreement, as shown in Fig. S26. Importantly, the interfacial surface area of the bubbles was found to decrease as the volume percentage of IL increased. This is counterintuitive since the surface tension of bubbles decreases as the concentration of IL increases. The decline in interfacial surface area is likely caused by saturation of the liquid surface with IL molecules. At high IL concentrations, the CO2–IL interactions likely dominate, thereby reducing the interactions of CO2 with the OH− ions and ultimately slowing the chemisorption process. Similar behavior was also observed when BMIM NTf2 was replaced with HMIM NTf2, and the fitting results are provided in Fig. S27 and Table S10. This effect can be rationalized by the findings of Mezger et al.134 who reported that high IL concentrations lead to layering of alternate cations and anion-enriched regions at the liquid interface. Such layers can often decrease mass transfer coefficients and ultimately impact the diffusion coefficients of gas molecules.135
Fig. 5(b) shows the chemical conversion of KOH in the presence of BMIM-based ILs. Almost all BMIM-based ILs show enhanced CO2 chemisorption kinetics. The BMIM-based IL system could achieve 50% conversion at 0.68 min, 0.95 min, 1.6 min, and 2.8 min for BMIM OTf, BMIM I, BMIM DCA, BMIM NTf2, and BMIM PF6. For BMIM-based ILs, the rate constants are shown in Fig. S25 and follow the decreasing trend: I > OTf > DCA > NTf2 > PF6.
Fig. 5(c) shows the chemical conversion of KOH in the presence of HMIM-based ILs. Similarly, all HMIM-based ILs show enhanced CO2 chemisorption kinetics compared to chemisorption in alkaline DES alone. The HMIM-based IL system achieved 50% conversion at 0.94 min, 1.3 min, and 1.4 min for HMIM PF6, HMIM I, HMIM OTf, and HMIM NTf2, respectively. For HMIM-based ILs, the rate constant follows almost the same order as what has been observed before. The 1st-order rate constants follow the order HMIM OTf > HMIM NTf2 > HMIM I > HMIM PF6, as shown in Fig. S28. When compared with a fixed cation such as NTf2 and a varying cation, the 1st order rate constants follow a general order: EMIM > BMIM > HMIM, as observed in Table S8. Overall, the current study revealed that the nature of anion and cation also plays an important role in the enhancement of CO2 chemisorption rates. Furthermore, the choice of anion seems to become less significant towards the improvement of CO2 capture kinetics in ILs with cations with longer alkyl chain substitutions.
Although the scope of the present study does not include a comprehensive techno-economic analysis (TEA), detailed kinetic modeling, or process-scale simulations for DES–IL systems, if future studies establish process feasibility, the use of volume percentages may be practically relevant at an industrial scale. However, differences in rate enhancement arising from various anions and cations may also reflect differences in molar loading. To account for this, we performed experiments under equimolar conditions, maintaining a fixed ionic liquid concentration of 0.1 M and a 0.1 M hydroxide concentration in the DES capture solution. As shown in Fig. S29, DCA− anion-based ILs consistently exhibit slightly slower kinetics compared to OTf− and NTf2− based ILs, although they still provide measurable enhancement in CO2 hydroxylation relative to the IL-free system.
Based on Fig. S29, one might argue that at equimolar IL concentrations, CO2 solubility governs the observed enhancement. This interpretation is reasonable for DCA− based ILs, which are known to exhibit lower CO2 solubility and correspondingly reduced enhancement relative to OTf− and NTf2− based ILs.136,137 However, distinguishing between OTf− and NTf2− anions based solely on CO2 solubility becomes less straightforward. For EMIM-based ILs, NTf2− shows significant enhancement over OTf−, whereas for HMIM-based ILs, OTf− outperforms NTf2−. In contrast, BMIM-based ILs exhibit comparable performance for both OTf− and NTf2− anions. These trends suggest that CO2 solubility alone cannot consistently explain the observed rate enhancements across different cation–anion combinations.
Another possible explanation is viscosity effects. However, accounting for viscosity also does not support the observation in Fig. S29. For example, BMIM DCA, which has the least viscosity of approximately 28 cP, exhibits the least enhancement among the ILs studied. Interestingly, although the viscosity of BMIM OTf is roughly twice that of BMIM NTf2−, both ILs display similar enhancement in CO2 hydroxylation rates. While some differences are observed between EMIM+ and HMIM+ cation based-ILs, there is no systematic correlation between viscosity and rate enhancement. Therefore, the argument that viscosity tracks with the enhancement observed in these experiments is not supported by the data.
Collectively, these results indicate that the observed rate enhancement cannot be attributed solely to CO2 solubility or viscosity effects. Instead, we emphasize that the enhancement arises from a combination of catalytic, interfacial, electrostatic, and solubility effects introduced by ILs. A more detailed molecular-level investigation, potentially involving DFT or molecular mechanics simulations of the gas–liquid interfacial structure, would be required to fully deconvolute these contributions, but such an analysis lies beyond the scope of the present work.
| This journal is © The Royal Society of Chemistry 2026 |