Carbon Dioxide Solubility in Three Bis Tri (Fluromethylsulfonyl) Imide-Based Ionic Liquids

This study delves into the necessity of mitigating carbon dioxide (CO2) emissions, focusing on effective capture methods to combat global warming by investigating the solubility of CO2 in three ionic liquids (ILs), 1-Decyl-3-MethylimidazoliumBis (Trifluromethylsulfonyl Imide) [IL1], 1-Hexadecyl-3-Methyl imidazoliumBis (Trifluromethylsulfonyl Imide) [IL2] and Triethytetradecyl Ammonium Bis (Trifluromethylsulfonyl Imide) [IL3]. Solubility experiments were conducted at (30, 50 and 70) °C with pressures up to 1.5 MPa. The research shows [IL2] as the superior candidate for CO2 capture, with its longer alkyl chain, and is confirmed by its lower Henry’s Law constant. Utilizing the Peng Robinson equation of state, the study correlates well with the solubility measurements using three mixing rules. The study reveals promising results for IL1, IL2 and IL3 surpassing all other published ionic liquids including Selexol/Genesorb 1753, except for 1-Methyl-3-octylimidazolium bis(trifluoromethylsulfonyl)imide. Insights into the enthalpy and entropy of absorption underscore the significant impact of IL structure on CO2 solubility, emphasizing the potential of tailored ILs for advanced carbon capture strategies. In summary, this research highlights [IL2] as the optimal choice for CO2 capture, offering valuable contributions to the ongoing efforts in combating climate change.


Introduction
In the past two decades, the escalation of temperatures, largely because of the persistent elevation of atmospheric CO 2 levels, has presented noteworthy concerns for the world in general [1].This phenomenon has prompted widespread apprehension among researchers, scientists, and environmentalists regarding the broader implications of climate change.Consequently, there is a concerted effort within the research community to develop improved solvents or technologies capable of efficiently capturing CO 2 .A significant amount of CO 2 release originates directly from the burning of fossil fuels, particularly hydrocarbons like natural gas, oil and coal.Various technologies, including membrane separation/permeation, adsorption, absorption, and cryogenic distillation, are utilized for carbon dioxide (CO 2 ) capture from natural gas or flue gases.Among these methods, absorption, noted for its cost-effectiveness and energy efficiency, is widely adopted on a large scale [2].
Over the years, numerous mechanisms and processes have been explored for CO 2 absorption, contingent upon the partial vapor pressures of the gas mixture.According to Zhang (2021) [2], physical solvents and amines have garnered significant attention due to their effectiveness and relatively lower costs.Physical solvents, in particular, offer advantages over amines, notably requiring less energy for solvent regeneration and exhibiting non-corrosive properties that prolong equipment lifespan, thus reducing maintenance expenses.Ionic liquids can be used for natural gas sweetening and mixed with amines or used in flue gas treatment.
Ionic liquids (ILs) represent a category of physical solvents that utilize intermolecular forces or functional group incorporation for CO 2 capture.Their intrinsic properties, including low volatility, high thermal stability, and chemical resilience, offer distinct advantages.Minimal volatility facilitates nearly solvent-loss-free regeneration processes, thereby reducing operating costs.Additionally, high thermal and chemical stability minimizes degradation and corrosion risks, ensuring prolonged equipment durability.This technology addresses shortcomings observed in chemical solvents like Monoethanolamine (MEA), such as excessive energy consumption and solvent loss [3], potentially serving as a potential MEA substitute [4].
Ionic liquids are typically classified into two types: task-specific ILs (TSILs) and nonfunctionalized room temperature ILs (RTILs).The key difference lies in RTILs' ability to function as conventional physical absorbents, governed by Henry's Gas constant, whereas TSILs demonstrate both chemical and physical CO 2 solubility.This dual capability enhances TSILs' capacity for CO 2 capture, as elucidated by Vadillo et al. (2022) [5].
The objective of this investigation is to assess the performance of three ionic liquids (ILs), categorized as promising physical solvents, in capturing CO 2 across a range of temperatures (30 • C to 70 • C) and pressures relevant to industrial applications, up to 1.5 MPa utilizing an intelligent gravimetric microbalance (IGA-003).

CO 2 Solubility Validation Test
Similar to the density validation, tests were conducted to confirm the repeatability of the IGA-003 and its adherence to the relevant experimental Standard Operating Procedures (SOPs).During this validation phase, the IL, [BMIM][BF4], was employed to evaluate CO 2 absorption.A comparative study was performed, aligning the obtained values with those documented by Shiflett and Yokozeki (2005) [6], who employed the same ionic liquid and referenced density values from a publication (Figure 1).advantages over amines, notably requiring less energy for solvent regeneration and exhibiting non-corrosive properties that prolong equipment lifespan, thus reducing maintenance expenses.Ionic liquids can be used for natural gas sweetening and mixed with amines or used in flue gas treatment.Ionic liquids (ILs) represent a category of physical solvents that utilize intermolecular forces or functional group incorporation for CO2 capture.Their intrinsic properties, including low volatility, high thermal stability, and chemical resilience, offer distinct advantages.Minimal volatility facilitates nearly solvent-loss-free regeneration processes, thereby reducing operating costs.Additionally, high thermal and chemical stability minimizes degradation and corrosion risks, ensuring prolonged equipment durability.This technology addresses shortcomings observed in chemical solvents like Monoethanolamine (MEA), such as excessive energy consumption and solvent loss [3], potentially serving as a potential MEA substitute [4].
Ionic liquids are typically classified into two types: task-specific ILs (TSILs) and nonfunctionalized room temperature ILs (RTILs).The key difference lies in RTILs' ability to function as conventional physical absorbents, governed by Henry's Gas constant, whereas TSILs demonstrate both chemical and physical CO2 solubility.This dual capability enhances TSILs' capacity for CO2 capture, as elucidated by Vadillo et al. (2022) [5].
The objective of this investigation is to assess the performance of three ionic liquids (ILs), categorized as promising physical solvents, in capturing CO2 across a range of temperatures (30 °C to 70 °C) and pressures relevant to industrial applications, up to 1.5 MPa utilizing an intelligent gravimetric microbalance (IGA-003).

CO2 Solubility Validation Test
Similar to the density validation, tests were conducted to confirm the repeatability of the IGA-003 and its adherence to the relevant experimental Standard Operating Procedures (SOPs).During this validation phase, the IL, [BMIM][BF4], was employed to evaluate CO2 absorption.A comparative study was performed, aligning the obtained values with those documented by Shiflett and Yokozeki (2005) [6], who employed the same ionic liquid and referenced density values from a publication (Figure 1).

ILs CO2 Solubility
The CO2 solubility of IL1, IL2 and IL3 was measured at various temperatures and pressures up to 1.5 MPa.The obtained data for the two ILs were plotted in Figure 2 and recorded in

ILs CO 2 Solubility
The CO 2 solubility of IL1, IL2 and IL3 was measured at various temperatures and pressures up to 1.5 MPa.The obtained data for the two ILs were plotted in Figure 2 and recorded in Table 1 below.Standard uncertainty u(x) = 0.006, standard uncertainty u(T) = 0.1 K, standard uncertainty u(P) = 0.0008 MPa.

Simulation Interaction Parameters
The experimental data were correlated using three mixing rules: (a) PR + vdW1, (b) PR + vdW2, and (c) PR + WS + NRTL.Table 2 below provides an overview of the calculated interaction binary parameters and their corresponding average absolute deviations (AADs %) for each mixing rule.The average absolute deviations (AADs %) for the vdW1, vdW2, and WS-NRTL rules applied to CO 2 absorption in these ILs were determined to be 6.52%, 0.67%, and 0.62%, respectively, for IL1, 6.21%, 1.16%, and 1.10% for IL2, and 5.14%, 0.99%, and 0.94% for IL3.Notably, the Wong-Sandler (WS-NRTL) mixing rule exhibited the lowest average absolute deviation in both cases, indicating its superiority as the preferred choice among these options.

Henry's Law Constant, Enthalpies and Entropies
The Henry's Law constants (H) were evaluated through a comprehensive analysis involving the plotting of fugacity against mole fraction and the subsequent fitting of a secondorder trend line.This analysis enabled the derivation of the slope of the second-order equation, which serves as an essential parameter in determining Henry's Law constant [7].Similarly, the estimation of entropy involved plotting the natural logarithm of Henry's Law constant against the natural logarithm of temperature (T) for each temperature of the ionic liquid.The resulting negative enthalpy variation signifies an exothermic absorption process.The operation of the absorption of a gas in a liquid results in a decrease in the volume of the molecules from the gas phase to the liquid phase (equivalent to compression), leading to a reduction in the entropy.The negative entropy variation with the increase in temperature, as illustrated in Table 3 below, results from the decrease in the number of absorbed CO 2 moles and a reduction in the solubility.
H (MPa) Figure 4 showcases a notable CO 2 absorption capacity exhibited by the studied ionic liquids in comparison to Selexol/Genesorb 1753 [8], a high-performance solvent utilized in gas treatment facilities.Given their superior CO 2 absorption performance, the ILs investigated in this research can be regarded as promising solvents.
Molecules 2024, 29, x FOR PEER REVIEW 7 of 13 in gas treatment facilities.Given their superior CO2 absorption performance, the ILs investigated in this research can be regarded as promising solvents.

Materials
Table 4 lists the detailed ILs used in this research along with their purity, chemical structures and nomenclature.

Density Measurement
Before and after each experiment, at the corresponding temperature, the Anton Paar Density and Speed of Sound instrument (DSA 5000, Anton Paar, Graz, Austria) was tested

Materials
Table 4 lists the detailed ILs used in this research along with their purity, chemical structures and nomenclature.in gas treatment facilities.Given their superior CO2 absorption performance, the ILs investigated in this research can be regarded as promising solvents.

Materials
Table 4 lists the detailed ILs used in this research along with their purity, chemical structures and nomenclature.

Density Measurement
Before and after each experiment, at the corresponding temperature, the Anton Paar Density and Speed of Sound instrument (DSA 5000, Anton Paar, Graz, Austria) was tested Molecules 2024, 29, x FOR PEER REVIEW 7 of 13 in gas treatment facilities.Given their superior CO2 absorption performance, the ILs investigated in this research can be regarded as promising solvents.

Materials
Table 4 lists the detailed ILs used in this research along with their purity, chemical structures and nomenclature.

Density Measurement
Before and after each experiment, at the corresponding temperature, the Anton Paar Density and Speed of Sound instrument (DSA 5000, Anton Paar, Graz, Austria) was tested in gas treatment facilities.Given their superior CO2 absorption performance, the ILs investigated in this research can be regarded as promising solvents.

Materials
Table 4 lists the detailed ILs used in this research along with their purity, chemical structures and nomenclature.

Density Measurement
Before and after each experiment, at the corresponding temperature, the Anton Paar Density and Speed of Sound instrument (DSA 5000, Anton Paar, Graz, Austria) was tested

Materials
Table 4 lists the detailed ILs used in this research along with their purity, chemical structures and nomenclature.

Density Measurement
Before and after each experiment, at the corresponding temperature, the Anton Paar Density and Speed of Sound instrument (DSA 5000, Anton Paar, Graz, Austria) was tested

Density Measurement
Before and after each experiment, at the corresponding temperature, the Anton Paar Density and Speed of Sound instrument (DSA 5000, Anton Paar, Graz, Austria) was tested with air and double-distilled water at atmospheric pressure to ensure the error was less than 0.00005 g•cm −3 when compared to the reference values for water and air stored in the instrument's database.To further ensure the accuracy of the density measurements for the ILs utilized in this study, the density (m/v) of N-methyldiethanolamine (MDEA) with a concentration purity of ≥99% was initially determined using the Anton Paar density meter.These measured data were then compared with those published by Karunarathne et al. (2020) [7] for MDEA with a similar purity level.The comparison revealed an average deviation (AAD) of 0.014%, indicating excellent repeatability/accuracy in the density measurements.Subsequently, the densities of three specific ILs were determined at atmospheric pressure conditions (P = 98.3 kPa) and temperatures ranging from 20 • C to 70 • C using the density meter.The obtained data are provided in Figure 5 and Table 5 below.
mospheric pressure conditions (P = 98.3 kPa) and temperatures ranging from 20 °C to 70 °C using the density meter.The obtained data are provided in Figure 5 and Table 5 below.Standard uncertainty u (ρ) = 0.0002 g/cm 3 ; standard uncertainty u (T) = 0.01 K; standard uncertainty u (P) = 0.8 kPa.

Solubility Analysis
The IGA-003 analyzer, manufactured by Hiden Isochema Ltd. (Warrington, UK), incorporates cutting-edge components to facilitate gas sorption experiments effectively.A Polyscience water bath, providing precise temperature control within the reactor chamber, up to 343 K in this study, ensures optimal conditions for accurate gas sorption measurements.This system integrates a water jacket around the reactor chamber to provide equal heat/temperature all around.Within the core vessel of the reactor chamber, the IL is placed, creating a controlled environment for the experimental procedure.The  Standard uncertainty u (ρ) = 0.0002 g/cm 3 ; standard uncertainty u (T) = 0.01 K; standard uncertainty u (P) = 0.8 kPa.

Solubility Analysis
The IGA-003 analyzer, manufactured by Hiden Isochema Ltd. (Warrington, UK), incorporates cutting-edge components to facilitate gas sorption experiments effectively.A Polyscience water bath, providing precise temperature control within the reactor chamber, up to 343 K in this study, ensures optimal conditions for accurate gas sorption measurements.This system integrates a water jacket around the reactor chamber to provide equal heat/temperature all around.Within the core vessel of the reactor chamber, the IL is placed, creating a controlled environment for the experimental procedure.The experiment gas-carbon dioxide (CO 2 )-is introduced into the controlled system to study their sorption characteristics.
Pressure regulation maintains precise pressure conditions for gas injection into the reactor chamber.The Mass Flow Control (MFC) system plays an important role in controlling the rate of experiment gases getting into the reaction chamber, maintaining uniform gas flow for accurate absorption study.The IGA-003 microbalance takes note of the change in ionic liquid mass during gas absorption, providing essential data on ILs' sorption behavior.In addition to the microbalance (IGA-003), a weight counterbalance enclosure compensates the weight of the reaction chamber and other components like the weight of the string, buoyancy, etc., allowing the IGA to measure specific weight changes in the sample material only during gas sorption precisely and accurately.The experimental procedure is automated and overseen by the data acquisition (DAQ) and control system, ideally operated through an onsite computer.This integrated system manages variables such as weight changes, pressure, gas flow, and temperature, maintaining the repeatability, efficiency, and accuracy of experiments performed with the IGA-003.

Thermodynamic Modelling
We utilized the Peng-Robinson Equation of State (PR EoS), as indicated in Equation ( 1), a well-established model renowned for its accuracy in describing the phase behavior of fluid mixtures, particularly at elevated pressures and temperatures.Through this model, we were able to estimate the absorption of CO 2 over a range of pressures up to 1.5 MPa and temperatures reaching 70 • C. The selection of the PR EoS was deliberate, as it leverages critical components of both liquids and gases involved in the system, aligning closely with our research objectives.
The coefficients of this model were predicted using various mixing rules, as stated below [7]: • van der Waals one single binary interaction parameter; • van der Waals two binary interaction parameters; • NRTL model combined with Wong-Sandler mixing rules (WS-NRTL).

van der Waals Mixing Rules
Two distinct mixing rules devised by van der Waals, namely van der Waals two (vdW2) and van der Waals one (vdW1), were employed to estimate the mixture variables represented by a m and b m [9].The vdW1 mixing rule involves estimating a single interaction parameter (lij), whereas vdW2 entails estimating two (2) binary interaction parameters (lij and kij).The parameter a m for vdW2 and vdW1 mixing rules was determined using Equation (2), while Equation (5) facilitated the estimation of a ij based on the temperaturedependent binary interaction parameter, kij.The co-volume factor bm for both vdW2 and vdW1 was estimated through Equations ( 3) and (4), respectively.Furthermore, the interaction parameter bij was determined using Equation (6). where NRTL model combined with Wong-Sandler mixing rule (WS-NRTL).
The modeling approach in this study integrates Wong-Sandler mixing rules, incorporating Equations ( 7) and ( 8) to determine the liquid-gas mixture's attractive force parameter 'a' and co-volume parameter 'b'.Utilizing the Non-Random Two Liquid (NRTL) model, the activity coefficient and excess Gibbs energy are computed, as specified in Equations ( 13)-( 19).Binary interaction parameters (τ_ji, τ_ij, τ_ij_kij) are instrumental in estimating the mixture parameters, where τ_ji and τ_ij represent the NRTL parameters, and g_ij and g_jj denote the interaction energies between molecules 'i' and 'j', as elucidated by Nath and Henni (2020) [9].Additionally, the value of α in Equations ( 14) and ( 15) is arbitrarily set at 0.3 in this study.

Critical Properties Calculations
To utilize the model effectively, it is imperative to possess thermo-critical properties for both the gas and the solvent under investigation.In our study, we employed the modified Lydersen-Joback-Reid group contribution method [9] to ascertain the critical temperature (Tc), the acentric factor (ω) and the critical pressure (Pc) of the ILs.The Lydersen-Joback-Reid method is a chemical group contribution method revisited and updated by Valderrama and Rojas [10] for the prediction of the critical properties and the acentric factor of the ionic liquids.The full description and an Excel file are reported in reference [10].The following equations are used: where A = 0.5703, B = 1.0121,C = 0.2573, D = 6.75 and Pb = 1.01325 bar.
In the equations, M is in g/mol, Tb and Tc are in K and Pc in bar, and n is the number of different groups.A table in reference [9] provides the values of Mi, ∆Tc, ∆Pc, ∆T b for a variety of chemical groups for ILs with and without rings.Table 6 shows a summary of the critical properties of the ionic liquids studied.

Binary Interaction Parameter Optimization
The thermodynamic models employed in this research were developed by Dr. Nath and Dr. Kazi as reported in reference [9] using the MATLAB R2024a software, utilizing a bubble point algorithm.As illustrated in Equation (24), the optimization of binary interaction parameters was carried out using the Nelder-Mead simplex method, utilizing the 'fminsearch' function integrated within MATLAB to minimize the error in the objective function.Given the likelihood of experimental data inaccuracies at lower pressures, the optimization of binary interaction parameters was conducted within a pressure range of 0.1 to 1.5 MPa.

Err
3.6.Henry's Law Constant, Enthalpy of Absorption and Entropy of Solvation Henry's Law constant (H) for the two ILs was evaluated by analyzing the slope of the second-order polynomial of a plot with the mole fractions (x) against the fugacity (f) of CO 2 at all temperatures.This constant, denoted as H i , represents the ratio of the solute (i) fugacity to its mole fraction in the solvent (j) at infinite dilution, occurring at a specific temperature within a specific pressure scope.Within this research, f V i and f L i signify the solute's fugacity in the vapor and liquid phases, respectively, while y i and x i imply the mole fractions of the gas in the vapor and liquid phases, respectively, as outlined by Huseynov (2014) [11].After the Henry's Law constant (H) is determined from the absorption data points, the subsequent procedure involves evaluating the entropy of solvation (∆S ∞ ) at infinite dilution and the enthalpy of absorption (∆H ∞ ) at infinite dilution using the provided equations.

Discussion and Conclusions
The observed trends in Henry's Law constants for CO 2 absorption in IL1, IL2 and IL3 reveal a consistent increase with temperature.This thermal dependence suggests the enhanced solubility of CO 2 gas in ionic liquids at lower temperatures, aligning with typical gas solvation behavior.
The alkyl chain lengths and anion fluorination play crucial roles in determining the absorption behavior of CO 2 in the studied ionic liquids (IL1, IL2 and IL3).IL1, composed of 1-decyl-3-methylimidazolium bis(trifluoromethylsulfonyl) imide, possesses a shorter alkyl chain compared to both IL2 and IL3.This difference in alkyl chain length significantly impacts the solvation properties of the ionic liquids, influencing their interaction with CO 2 molecules.
The observed trends in Henry's Law constants reveal interesting insights into the effect of alkyl chain length on gas solubility.Generally, relatively longer alkyl chains, as in IL2, tend to promote stronger interactions with gas molecules, resulting in lower Henry's Law constants compared to shorter alkyl chain counterparts, such as IL1 and IL3.This phenomenon can be attributed to the increased surface area and van der Waals interactions provided by longer alkyl chains, facilitating greater solubility of CO 2 in IL2.
Furthermore, the fluorination of the anion in IL1, IL2 and IL3 introduces additional effects on gas solvation behavior.The presence of fluorine atoms in the anion enhances the polarity and can induce specific interactions with gas molecules, contributing to the overall solvation process.This effect is evident in the observed Henry's Law constants, where IL1, IL2 and IL3 exhibit impressive performance in terms of CO 2 absorption compared to non-fluorinated counterparts reported in the published literature.
Regarding the obtained thermodynamic parameters, the negative enthalpy values (∆H ∞ ) obtained for CO 2 absorption in IL1, IL2 and IL3 indicate exothermic solvation processes.The enthalpy value calculated was −12.80 kJ/mol in IL1, −10.33 kJ/mol in IL2 and −13.79 kJ/mol in IL3, further confirming the exothermic nature of the absorption processes.Additionally, the negative entropy values imply a decrease in disorderliness within the system upon solvation, reflecting the ordered arrangement of solvent molecules around the solute during absorption.
In conclusion, the alkyl chain lengths and anion fluorination effect significantly influence the absorption behavior of CO 2 in IL1, IL2 and IL3.The observed trends in Henry's Law constants highlight the importance of molecular structure and interactions in gas solvation processes.The exothermic nature of the absorption processes and the decrease in entropy upon solvation further underscore the thermodynamic aspects of gas solvation in ionic liquids.
Comparative analysis with literature data highlights the favorable performance of IL1, IL2 and IL3 in terms of Henry's Law constants for CO 2 absorption.Despite their impressive performance, further investigation into other ionic liquid formulations, such as 1-octyl-3-methyl imidazolium bis(trifluoromethyl) imide, may provide valuable insights into optimizing gas absorption processes.
Overall, the detailed investigation provides valuable insights into the thermodynamics of gas solvation in IL1, IL2 and IL3, offering potential avenues for further research and optimization in gas separation technologies.The observed trends underscore the promise of fluorinated ionic liquids for CO 2 absorption applications, while also suggesting the need for additional studies to elucidate the underlying molecular mechanisms and optimize the performance of these systems.

Figure 4
Figure4showcases a notable CO2 absorption capacity exhibited by the studied ionic liquids in comparison to Selexol/Genesorb 1753[8], a high-performance solvent utilized

Figure 4 .
Figure 4. Comparison of Henry's Law constants for CO2 in ILs obtained in this work and Selexol as reported by Aravind et al. (2012) [8].

Figure 4 .
Figure 4. Comparison of Henry's Law constants for CO 2 in ILs obtained in this work and Selexol as reported by Aravind et al. (2012) [8].

Figure 4 .
Figure 4. Comparison of Henry's Law constants for CO2 in ILs obtained in this work and Selexol as reported by Aravind et al. (2012) [8].

Figure 4 .
Figure 4. Comparison of Henry's Law constants for CO2 in ILs obtained in this work and Selexol as reported by Aravind et al. (2012) [8].

Figure 4 .
Figure 4. Comparison of Henry's Law constants for CO2 in ILs obtained in this work and Selexol as reported by Aravind et al. (2012) [8].
et al 2012) IL3Molecules 2024, 29, x FOR PEER REVIEW 7 of 13 in gas treatment facilities.Given their superior CO2 absorption performance, the ILs investigated in this research can be regarded as promising solvents.

Figure 4 .
Figure 4. Comparison of Henry's Law constants for CO2 in ILs obtained in this work and Selexol as reported by Aravind et al. (2012) [8].

Figure 5 .
Figure 5. Densities of the ILs used in this work.

Figure 5 .
Figure 5. Densities of the ILs used in this work.

Table 2 .
Optimized binary interaction parameters for vdW1, vdW2 and WS-NRTL with their corresponding average absolute deviation percentages (AAD %).

Table 3 .
Henry's Law constant, enthalpy and entropy of solvation between CO 2 and ILs.Law constant (H) implies the IL's capability to absorb a larger quantity of CO 2 .In this study, [IL2] demonstrated the highest CO 2 absorption capacity, while [IL1] exhibited the least.When juxtaposing the estimated Henry's Law constants from this research with those of other ILs documented by previous researchers, [OMIM] [TF2N] was the only published IL with the highest solubility among those reported in this study, as seen in Figure3below.

Table 3 .
Henry's Law constant, enthalpy and entropy of solvation between CO2 and ILs.Law constant (H) implies the IL's capability to absorb a larger quantity of CO2.In this study, [IL2] demonstrated the highest CO2 absorption capacity, while [IL1] exhibited the least.When juxtaposing the estimated Henry's Law constants from this research with those of other ILs documented by previous researchers, [OMIM] [TF2N] was the only published IL with the highest solubility among those reported in this study, as seen in

Table 4 .
Detailed list of all ionic liquids.

Table 4 .
Detailed list of all ionic liquids.

Table 4 .
Detailed list of all ionic liquids.

Table 4 .
Detailed list of all ionic liquids.

Table 4 .
Detailed list of all ionic liquids.

Table 4 .
Detailed list of all ionic liquids.

Table 5 .
Density of the ILs used in this work.

Table 5 .
Density of the ILs used in this work.