Abstract
The intimidating level of anthropogenic CO2 in the atmosphere responsible for global warming and erratic weather conditions needs to be addressed on a priority basis. Different kinds of materials were used to capture CO2 to curtail the alarming and drastic effects of global warming. An ionic liquid (IL) 1-butyl-3-methylimidazolium methanesulfonate [C4mim][CH3SO3] was chosen, owing to its unique and efficient characteristics required for CO2 capture. Thermos-physical characteristics such as sigma surface, sigma profile, and sigma potential are calculated from the COSMO-RS model independent of any kind of experimental or coefficient data as an input. The mandatory information required for the interaction of IL with CO2 was obtained from this model. The COSMO-RS model depends upon unimolecular quantum chemical analysis associated with statistical thermodynamics, molecular structure, and conformation. The structural confirmation of [C4mim][CH3SO3] IL was performed by FTIR, 1H NMR, and 13C NMR spectroscopic methods. Spectrochemical properties are calculated by FTIR, NMR, UV-visible, and fluorescence. Maximum CO2 solubility performed at room temperature (RT) and 45 bar was found to be ~2.7 mmol/g. The uptake of CO2 indicates the presence of sulphur-functionalized anions and bulky alkyl groups in IL’s significant affinity towards CO2. According to hysteresis-based classification, CO2 sorption and desorption follows type H3 classification, which indicates the presence of microporous and mesoporous in the IL sample. The effect of functionalized anions and alkyl groups on CO2 capture is highlighted in this study. The present study is aimed at providing a detailed overview related to theoretical and experimental study and application in terms of CO2 capture of IL.
1. Introduction
An alarming use of fossil fuels to produce energy for different purposes releases massive scale of CO2 primarily responsible for climatic fluctuations all over the world due to global warming [1–3]. This huge CO2 emission to the atmosphere is increasing rapidly due to enormous development taking place everywhere. A substantial part of CO2 emission is coming from power plants (41%) followed by transport (23%), industries (20%), construction (4%), and other sectors (12%) [4]. Excessive CO2 emission is produced by combustion to generate electricity and heat [5, 6]. The expansion of CO2 emission and its harmful effects such as global warming and erratic climatic conditions concern researchers, scientists, and academicians all over the world. To minimize the intimidating effects of global warming, a simple, quick, and efficient method is required to decrease the emission of CO2 into the atmosphere [7, 8]. The precombustion, postcombustion, and oxycombustion methods are currently being implemented to reduce the emission of CO2 into the atmosphere [9–11]. In addition, to attain methane with low CO2 content emission, impurities are eliminated to use methane as a biomethane, a useful replacement for natural gas [12–14].
Three important types of amines including monoethanolamine (MEA), diethanolamine (DEA), and methyldiethanolamine (MDA) are utilized to capture CO2 also known as early amine technology for CO2 capture [15, 16]. MEA (30%) together with water was used in industry for CO2 capture owing to its excellent properties [17]. In addition, chemical reactions utilized higher CO2 enthalpy to show the high-energy amount in the transformation process [18]. The amine-based technology utilized a CO2 capture process in the current power area to lower energy by ~25-40% [19]. Besides several advantages of amine-based technology, many disadvantages such as the requirement of excessive energy, instrument disintegration, adsorbent volatility, and low CO2 capture were also found [20]. Furthermore, this technology uses massive investment and postcombustion of CO2 uptake, which makes it economically not affordable [21]. To overcome the drawbacks of amine-based technology, efficient and economical kind of materials such as stable ionic liquids (ILs) were currently being used for CO2 capture [3, 22, 23]. The low energy requirement in the regeneration stage of ILs forces researchers worldwide to shift research from amine-based technology to IL-based technology [24]. For the CO2 capture process, stable room temperature ionic liquids (RTILs) also known as green solvents are used owing to their unique characteristics. Recently, a cost target of $40/tCO2 for efficient, innovative, and stable generation solvents was set by the Department of Energy (DOE) [25, 26]. Recently, different kinds of ILs such as 1-n-butyl-3-methylimi-dazolium hexafluorophosphate ([bmim][PF6]), 1-n-butyl-2,3-dimethylimidazolium hexafluorophosphate ([bmmim][PF6]), 3-butyl-1-methylimidazolium bis(trifluoromethane-sulfonyl)amide ([bmim][Tf2N]), 1-n-butyl-3-methylimidazolium tetrafluoroborate ([bmim][BF4]), 1-n-butyl-2,3-dimethylimidazolium tetrafluoroborate ([bm-mim][BF4]), 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([emim][Tf2N]), 1-butylpyridinium bis(trifluoromethanesulfonyl)amide ([bpy][Tf2N]), 1-propyl-2,3-dimethylimidazolium bis(trifluoromethane-sulfonyl)amide ([pdmim][Tf2N]), and 1-ethyl-2,3-dimethylimidazolium bis(trifluoromethylsulfonyl)imide ([emmim][Tf2N]) were utilized in order to understand the mechanism of CO2 solubility in imidazolium-based ILs at different temperatures and pressures. This kind of imidazolium-based ILs is considered important for CO2 capture owing to its unique characteristics including lower melting points, low viscosity, facile synthesis, and significant stability towards oxidative and reductive behaviors. Furthermore, it was discovered that anion plays the primary role to interact with the CO2 and cations dominate the secondary role in imidazolium-based ILs.
Screening model COSMO-RS based on the approach of density functional theory (DFT) was utilized to screen different types of solvents for CO2 capture [27–29]. Molecular structure in terms of conformation collected from electronic calculations of the desired compound is the isolated parameter mandatory for COSMO-RS calculations [30–32]. The sigma potentials, sigma surfaces, sigma profiles, phase equilibrium, and activity coefficient parameters were obtained from this screening model [33–35]. The nominal type of data required for different types of calculations makes the COSMO-RS type model interesting and ambitious as compared to thermodynamic calculations where several kinds of parameters are essential to fix the data to accomplish proper results. To perceive the solubility process based on the analysis of quantum calculations, the screening model uses statistical thermodynamics. In short, the advantages mentioned in this section cited from different literature sources make this model an absolute replacement to calculate parameters particularly the sigma surface, sigma potential, and sigma profile of the ILs.
Ionic liquid 1-butyl-3-methylimidazolium methanesulfonate [C4mim][CH3SO3] was investigated both theoretically and experimentally. Detailed examination to study charge density distribution, sigma profile, and sigma potential was performed by the computational-based conductivity model COSMO-RS. The sigma profile, potential, and surfaces are derived from quantum mechanical calculation with molecular and confirmation data only. The physical and chemical properties, quantitative and qualitative measurements, and identification of compounds were analyzed in this research study. FTIR, 1H NMR, 13C NMR, fluorescence, and UV-visible spectroscopic methods performed these measurements. Another aim of this research study was to perform CO2 capture at room temperature (RT) of IL, which may help in the future to minimize the harmful effects of global warming and climatic fluctuations. The COSMO-RS helps us to measure the sigma surface, sigma profile, and sigma potential which provides us virtual conductor environment required for the prediction of the interaction of IL with CO2. Different types of characterizations such as 1H NMR, 13C NMR, and FT-IR confirmed the structure of IL. Another important finding related to this study was to find ~2.7 mmol/g of CO2 uptake at RT and an equilibrium pressure of 45 bar. As far our knowledge is concerned, imidazolium-based IL used in this study shows significant CO2 uptake as compared to values reported in the literature.
2. Experimental Section
2.1. Adsorbent
Ionic liquid (IL) 1-butyl-3-methylimidazolium methanesulfonate [C4mim][CH3SO3] was obtained from Sigma-Aldrich. Before characterization, IL dried at 313.15 K for one hour. Before sorption/desorption measurement, IL was activated by outgassed in situ ultrahigh vacuum conditions for 24 h at room temperature (RT).
2.2. Computational Method
The COSMO-RS software is an economical, quick, and adequate method to evaluate the thermodynamic, physical, and molecular characteristics of ILs [36]. This software works on unimolecular quantum study installed with authentic static thermodynamics that gives accurate data required for the calculation of molecular interactions in the ILs [37]. In this current report, COSMO-RS uses the BVP86/TZVP/DGA1 level of theory to screen charges on the surface of [C4mim][CH3SO3]. The structure of IL in the form of color coding 2D/3D, surface charge density, sigma profile, and sigma potential was calculated by using the COSMO-RS model. COSMO-RS uses an advanced thermodynamics path to elaborate the solubility process settled on the outcome of quantum chemical results. The parameters related to the qualitative description of IL are polarity, hydrogen bonding, lipophilicity, and hydrophilicity, which are estimated by using COSMO-RS three-dimensional (3D) sigma surface or screen and charge distribution. Furthermore, these kinds of properties make a model an important replacement for the calculation of activity coefficients.
2.3. Characterization
Fourier transform infrared (FTIR) spectroscopy pattern of [C4mim][CH3SO3] was obtained from Bruker ALPHA FT-IR spectrometer (Eco-ATR) using the KBr pellet process. The wavenumber range of the spectrum is 500-4000 cm-1. The 1H NMR (proton nuclear magnetic resonance) and 13C NMR (carbon 13 nuclear magnetic resonance) using Bruker AVANCE 400 NMR spectrometer were measured in a 5 mm precision glass NMR tube at RT in deuterated solvent DMSO-d6 or CDCl3 with tetramethylsilane (TMS) as the internal standard. In the splitting process, “s” is assigned as a singlet, “d” doublet, “dd” doublet of doublets, “t” triplet, and “m” multiplet. The coupling constant (J) in NMR is denoted as hertz (Hz). Both 1H NMR and 13C NMR patterns were analyzed by the Delta software package. The 1H NMR and 13C NMR analyses of the IL are as follows: 1H NMR (DMSO-d6): δ 0.90 (t, 3H, CH3), 1.25 (m, 2H, CH2), 1.76 (m, 2H, NCH2CH2), 2.33 (s, 3H, (SO3)CH3), 3.86 (s, 3H, NCH3), 4.17 (t, 2H, NCH2), 7.71 (s, 1H, ArH), 7.78 (s, 1H, ArH), 9.16 (s, 1H, ArH); 13C NMR (DMSO-d6): δ 13.73 (CH2CH3), 19.23 (CH2CH3), 31.82 (CH2CH2CH3), 36.18 (CH3S), 40.20 (NHCH3), 48.95 (NCH2CH3), 122.74 (NCH), 124.08 (NCH), 137.06 (NCHN). Fluorescence was recorded on Cary Eclipse fluorescence spectrometer (Agilent, USA) in a 1 cm cell assembled with a 150 W lamp and water bath. The emission spectrum was obtained in the range of 290 nm to 600 nm wavelength with 280 nm excitation wavelength using slit width 5 nm each. The fluorescence was performed in a methanol solvent. UV-visible absorption spectroscopy was performed in a T80 UV/visible spectrophotometer running from 300 to 650 nm at RT.
2.4. Carbon Dioxide Capture
High-pressure CO2 capture of [C4mim][CH3SO3] was executed at RT in a customized manometric instrument designed and built by the LMA group, now commercialized as iSorpHP by Quantachrome Instruments [38]. Before sorption/desorption measurements, IL was activated by being outgassed under ultrahigh vacuum conditions for 24 h at RT. Sorption/desorption isotherm was obtained at RT and 4.5 MPa. Gas expanded from a thermostat-dosing manifold of settled volume into the IL sample cell of the instrument. The instrument automatically measured the internal volume of the cell by using a helium-expansion calibration process. The volume of the gas discharged into the cell is measured by the pressure drop in the manifold. Postequilibrium is accomplished; the final pressure in the sample cell is recorded. The instrumental setup is established in an isothermal bath, and the inside temperature is administered and maintained by a proportional-integral-differential (PID) composed thermostat. The Redlich-Kwong equation is used to measure the accurate dead volume of the instrument. Approximately, 200 mg of sample was used for each CO2 sorption/desorption experiment. The CO2 was vented into the apparatus and maintained until equilibrium pressure was obtained. Post accomplishing equilibrium, CO2 gas isothermally discharged into a sample containing a steel reactor and a pressure drop of CO2 was recorded. The CO2 capture was measured from pressure drop using the ideal gas equation , where , , , , and have their usual meaning.
3. Results and Discussion
Figure 1 shows a chemical structure in the form of 2D/3D along with the COSMO-RS surface charge density of IL [C4mim][CH3SO3]. The sigma surface contains cation and anion with varying colors showing the strength of the hydrogen-bond donor and hydrogen-bond acceptor and nonpolarity. The color scale indicates the level of electropositivity and electronegativity nature of cation and anion. The blue area in the cation indicates a positive charge owing to its acidic hydrogen, and the red surface shows negative due to fluorine and oxygen, which are highly electronegative. Sigma () profile displayed in Figure 2 indicates a division of polar charges related to the molecular surface. The -profile shows charge densities, which are expressed to show a relative quantity of surface with polarity . The graph where () is plotted againstdisplays how far the molecular surface has traveled over the few intervals of polarity (). The information related to the effect of cation and anion in the IL was achieved from the sigma profile by utilizing the COSMO-RS model. The -profile of 1-butyl-3-methylimidazolium cation and methanesulfonate anion shown in Figure 2 extends from -0.03 (donor) to 0.03 (acceptor) with a maximum spike at about 22 and 15 (), respectively The region beyond +0.01 and –0.01 indicates intense polarity with a possibility to achieve hydrogen bonds (HBs). The area between ±0.01 is considered a nonpolar region. The highly polarized region beyond +0.01 belongs to imidazolium, hydrocarbon chain, and a methyl group. Similarly, a highly polarized region beyond -0.01 belongs to the methanesulfonate anion. Multiple peaks corresponding to anion may be due to the presence of highly negative oxygen and sulphur element. In addition, both cation and anion molecules can act as a donor and acceptor to make hydrogen bonds. The peak strength of the cation in the nonpolar region is high compared to the peak intensity of the anion [39]. Figure 3 shows the sigma potential profile of [C4mim][CH3SO3] IL. Total energy calculated from sigma potential is composed of two sections: (a) recovering free energy, which assists to reinstate the molecule to its originality and energy, is positive, and (b) unrestricted hydrogen bonding energy is negative in nature because of the hydrogen bond between hydrogen bond donor and hydrogen bond acceptor. The final negative charge indicates control of hydrogen bonding, while the final positive charge shows the influence of free energy reconstruction. The COSMO-RS histogram shows three main areas: (i) hydrogen bond donor beyond the region of –0.01 (b) hydrogen bond acceptor beyond the region of +0.01 , and (c) nonpolar region between Regions between +0.01 and –0.01 indicate nonpolar nature taking part mainly in CH- interaction [40].



FTIR in the form of % transmittance vs. wavenumber (cm-1) displayed in Figure 4 shows different types of mixed peaks. The broad characteristic band at 2960-3097 cm-1 corresponds to the symmetric and asymmetric vibrational state of the imidazole ring ν(C–H) and aliphatic ν(C–H) groups. The different kinds of vibrational bands found in the FTIR spectrum are summarized in Table 1. The absorption bands such as imidazole ν(ring stretching), imidazole H–C–C and H–C–N bending, SO3 symmetric stretch, in-plane imidazole ring bending, out-of-plane C–H bending of imidazole ring, C-S stretch, SO3 symmetric bend, and imidazole C2–N1–C5 bending were found at 1570.28, 1186.14, 1049.35, 850.52, 754.28, and 624.64 cm-1, respectively, confirming the presence of [C4mim][CH3SO3] IL. Furthermore, no red shift or blue shift was observed in the FTIR spectrum.

The residual proton in deuterated solvent DMSO-d6 found at 2.50-2.52 ppm was used as 1H NMR external reference. The positions and existence of ambient hydrogen atoms or protons in the [C4mim][CH3SO3] are shown in Figure 5 and Table 2. Assigned 1H NMR of IL calculated from the MNova software package is displayed in Figure 6. The presence of distinct chemical shifts in the 1H NMR spectrum confirms the structure of IL. The chemical shift found at 0.90 ppm indicates protons associated with carbon and hydrogen in the form of the CH3 group in the structure. The broad chemical shifts found at 1.23-1.29 ppm and 1.76 ppm show protons associated with carbon and hydrogen in the form of CH2 and NCH2CH2 groups. The chemical shift at 2.33 and 3.86 ppm shows protons associated with carbon, S, and nitrogen in the form of SO3CH3 and NCH3 groups. Similarly, the broad chemical shifts found at 7.71, 7.78, and 9.16 ppm show protons associated with carbon and hydrogen in the form of the ArH group.


The 13C NMR spectra of IL [C4mim][CH3SO3] were obtained by using Bruker AVANCE 400 NMR spectrometer having N2 as a carrier gas to fix the position and ambiance of distinct type of carbon atoms present. Figure 7 followed by Table 3 shows distinct chemical shifts for varying types of carbons present in the DMSO-d6 as a solvent. Assigned 13C NMR of IL [C4mim][CH3SO3] calculated from the MNova software package is displayed in Figure 8. The 13C chemical shifts found at 39.36-40.61 ppm correspond to solvent DMSO-d6. The 13C chemical shifts found at 13.73, 19.23, 31.82, 36.18, 40.20, 48.95, 122.74, 124.08, and 137.06 ppm show the presence of carbon atoms in the form of -CH2CH3, -CH2CH2CH3, -CH3S, -NHCH3, NCH2CH3, -NCH, and NCHN groups. The presence of different types of carbon atoms in the spectrum confirms the structure of IL.


Figure 9 displays the fluorescence measurement of [C4mim][CH3SO3] IL in polar solvent methanol. The emission peak of methanol related to interaction with IL was found at 310 nm. The high intensity peak of methanol indicates the interaction of methanol with IL via cation-π interactions and hydrogen bonding [41]. In addition, the high peak intensity of methanol suggests solubility of IL in methanol may be due to the presence of sulphur and oxygen groups in the anion region. The excitation spectra match in shape and position almost like the absorption spectrum of methanol solvent suggesting emission measured origin from the same compound. The spectrophotometric analysis in terms of UV-visible in the range of 300-650 nm of [C4mim][CH3SO3] IL is shown in Figure 10. The absorption spectra were performed in a 1 cm quartz cell. The absorbance behavior of IL in the UV-visible range in methanol solution was recorded at RT. The absorption spectra show a prominent broad absorption peak at 230 nm corresponding to transition of CH2 unit and group of C=N and S=O.


The CO2 capture analysis of [C4mim][CH3SO3] IL presented in Figure 11(a) shows adsorption (○) and desorption (●) isotherms. The isotherm is executed at high pressure using a static high-pressure apparatus at RT. The Benedict-Webb-Rubin (mBWR) equation of state (EOS) is an extension of virial EOS employed for CO2 isotherm [42]. The maximum CO2 sorption of the IL was found ~2.7 mmol/g at RT and 45 bar. Figure 11(b) shows the dependence of the maximum CO2 absorption capacity as a function of equilibrium pressure. Furthermore, we found that CO2 absorption capacity is found to have a monotonous dependence on the equilibrium pressure which indicates that there are no significant adsorbate–adsorbate interactions at low experimental pressures. It was found that CO2 capture is dependent on the bonding of CO2 with the anion as suggested by the spectroscopic and molecular simulation studies. The presence of sulphur-functionalized atoms in the IL was found to be supportive for bonding between CO2 and sulphur containing substituents present in the anion [43]. Furthermore, more sulphur atoms in anions followed by the presence of ammonium group with bulky alkyl chain groups in the cation may also reinforce overall CO2 capture [44]. At high pressure, the uptake of CO2 is more as compared to low pressure because our results show more than 80% uptake found at high pressure. It was found that the alky chain of the IL boosts CO2 uptake when the alkyl chain length was elevated from ethyl to hexyl. Blanchard et al. also found that CO2 solubility increases when the alkyl chain length was elevated from butyl to octyl [45]. The CO2 absorption capacity of imidazolium-based ionic liquids such as [BBT][TFSI], [HMIM][TFSI], [BMIM][GLY], [EMIM][TFSI], [BBT][BF4], [BMIM][BF4], and [EMIM][BF4] at RT was found to be 0.2523, 0.2335, 0.2150, 0.1808, 0.1553, 0.0972, and 0.0795 (Mol (CO2 abs)/Mol (IL)) (https://www.mdpi.com/2297-8739/10/3/192). In addition, it was observed that triazolium-based ILs show higher CO2 absorption capacity as compared to imidazolium-based ILs with the same anion. Another series of imidazolium-based ILs such as [EMIM][Tf2N], [EMMIM][Tf2N], [BMIM][PF6], [BMMIM][PF6], and [BMMIM][BF4] were utilized to perform CO2 absorption [J. AM. CHEM. SOC. 2004, 126, 5300-5308]. It was found that bis(trifluoromethylsulfonyl)-based imide anion shows high CO2 absorption capacity. Furthermore, there is a little difference in CO2 absorption between ILs having the tetrafluoroborate or hexafluorophosphate anion.

4. Conclusion
To reduce the effects of global warming and erratic climatic conditions generated by the emission of CO2 into the atmosphere, an efficient adsorbing material 1-butyl-3-methylimidazolium methanesulfonate [C4mim][CH3SO3] ionic liquid (IL) was chosen to study theoretically and experimentally. The IL was analyzed by COSMO-RS, a conductor-screening model to obtain various distinct parameters including sigma surface, sigma profile, and sigma potential. The striking feature of this model is that no experimental data, coefficients, or parameters are required to obtain these properties. Characterizations such as 1H NMR, 13C NMR, and FT-IR were used to confirm the structural characteristics of IL. The spectrochemical properties were calculated by using FTIR, 1H NMR, 13C NMR, UV-visible, and fluorescence. The maximum CO2 uptake performed at RT and an equilibrium pressure of 45 bar was found to be ~2.7 mmol/g. According to the hysteresis-based classifications, hysteresis follows type H3 classification indicating the presence of microporous and mesoporous kinds of materials.
Data Availability
Data are available on request.
Conflicts of Interest
The authors declare that they have no conflict of interest.
Acknowledgments
This research work was funded by Institutional Fund Projects under grant no. IFPNC-003-135-2020. Therefore, authors gratefully acknowledge technical and financial support from the Ministry of Education and King Abdulaziz University, Jeddah, Saudi Arabia.