On the Diffusion of Ionic Liquids in ILs@ZIF-8 Composite Materials: A Density Functional Theory Study

Abstract

Recently, composite materials consisting of ionic liquids (ILs) and metal–organic frameworks (MOFs) have attracted a great deal of attention due to their fantastic properties. Many theoretical studies have been performed on their special structures and gas separation applications. Yet, the mechanism for the diffusion of ILs inside MOF channels still remains unclear. Here, the DFT calculations (e.g., rigid and relaxed potential energy surface, PES, scan) together with frontier orbital analysis, natural charge analysis, and energy decomposition analysis were performed to investigate the diffusion behavior of a typical IL, [C₄mim][PF₆], into the ZIF-8 SOD cage. The PES profiles indicate that it is quite difficult for the cation [C₄min]⁺ to diffuse into the cage of ZIF-8 through the pristine pores because of the large imidazole steric hindrance, which results in a large energy barrier of ca. 40 kcal·mol⁻¹ at the least. Interestingly, the PES reveals that a successful diffusion could be obtained by thermal contributions, which enlarge the pore size through swing effects at higher temperatures. For example, both [C₄mim]⁺ and [PF₆]⁻ could easily diffuse through the channel of the ZIF-8 SOD cage when the pore size was increased to 6.9 Å. Subsequently, electronic structure analyses reveal that the main interactions between [PF₆]⁻ or [C₄mim]⁺ and ZIF-8 are the steric repulsion interactions. Finally, the effects of the amounts of [C₄mim][PF₆] on the ZIF-8 structures were investigated, and the results show that two pairs of [C₄mim][PF₆] per SOD cage are the most stable in terms of the interaction between energies and structural changes. With these findings, we propose that the high-temperature technique could be employed during the synthesis of IL@MOF membranes, to enrich their family members and their industrial applications.

1. Introduction

Recently, composite materials consisting of ionic liquids (ILs) and metal–organic frameworks (MOFs) have attracted considerable attention because of their exceptional properties [1,2]; therefore, they have been used in many important fields, e.g., gas adsorption and separation, catalysis, and sensors [3,4,5,6,7,8,9,10,11,12,13]. In particular, the large number of adsorption sites of MOFs together with ILs play an important role in the gas mixture separation (e.g., CO₂/CH₄) [14,15,16,17]. For example, Ban et al. [18] loaded [C₄mim][TF₂N] into the cages of ZIF-8, and the results showed the adsorption selectivity for CO₂/CH₄ significantly increased from 7.5 to 41. Following this work, Kinik et al. [19] examined the [C₄mim][PF₆]/ZIF-8 system via the density functional theory (DFT) and Monte Carlo (MC) simulations. The authors found that the interactions between [C₄mim][PF₆] and ZIF-8 created certain new adsorption sites, which increases the ideal selectivity for CO₂/N₂ from 7.82 to 24.21. Burak et al. [20] obtained a similar result by incorporating [C₄mim][BF₄] into ZIF-8, of which the selectivity of 13.3 has been obtained for CO₂/N₂. Afterwards, a large number of systems have been reported with different combinations of ILs and MOFs [21].

In general, the reasons for ILs@MOFs composite materials with high gas separation or selectivity have been attributed to two main aspects. On the one hand, the rich hydrogen bond network of ILs [22,23] could effectively dissolve CO₂ molecules against other gases, e.g., CH₄ or N₂ [24]; hence, this improves the separation performance compared to pure MOF materials [25,26,27]. For example, Thomas et al. [28] performed DFT and Giant Canonical Monte Carlo (GCMC) simulations to investigate the selectivity of a series of IL@ZIF-8 systems. Their findings revealed that the hydrophobic fluorinated anions (e.g., [BF₄]⁻, [Tf₂N]⁻, and [PF₆]⁻) exhibited higher selectivity compared to hydrophilic non-fluorinated anions (e.g., [NO₂]⁻, [NO₃]⁻, and [SCN]⁻). Kavak et al. [29] investigated the CO₂ separation efficiency of five different IL-encapsulated MIL-53 systems, and they found that the CO₂/CH₄ separation selectivity of [C₄mim][PF₆]/MIL-53 was 2.8 times higher compared to the pure MIL-53. Zhang et al. [30] conducted an investigation on ILs with three different cation functional groups, and the authors revealed that ILs containing amino functional groups exhibited the highest CO₂ molecule affinity and CO₂/CH₄ adsorption selectivity.

On the other hand, it is believed that the introduction of ILs could also modify the structures of the MOFs, and subsequently increase the gas selectivity. Based on the experimental techniques (e.g., high-pressure X-ray diffraction, XRD) with molecular simulations, Fairen et al. [31] demonstrated that the imidazolate linker of the ZIF-8 changed because of the swing effects, and subsequently changed the gas selectivity [32]. Uzun and colleagues [33] coated 1-(2-hydroxyethyl)-3-methylimidazolium dicyanamide ([Hemim][DCA]) on the ZIF-8 surface to eliminate the nonselective voids, and the authors found that the CO₂ selectivity was increased by 45 times compared to the pure ZIF-8. Chang et al. [34] successfully constructed a zwitterionic MOF material (MOF UiO-66-SO₃⁻-NH₃⁺), and the adsorption capacity for CO₂ increased by 32–48%, because the positive and negative charges on the modified MOF material interacted with the ILs of [Emim][SCN]. Our previous theoretical studies [35,36] revealed that the ZIF-8 aperture configurations have a significant impact on the separation of CO₂/CH₄. Specifically speaking, the DFT calculations showed that the pristine ZIF-8 aperture (with a pore size of 3.4 Å) exhibits the best separation performance, which has the largest energy barrier difference between CO₂ and CH₄; additionally, the MD simulations revealed that ILs (e.g., [C₄mim][PF₆]) could effectively maintain pristine aperture configurations, and retained the high separation properties.

Yet, there are still open questions following the reported mechanistic studies; that is to say, how do the ILs diffuse into the cage of MOFs (e.g., ZIF-8)? And, how many ILs could be accommodated per cage? In this vein, we performed DFT calculations to investigate the diffusion behavior of ILs into MOF cages by taking [C₄mim][PF₆] and ZIF-8 as an example. Firstly, the potential energy surface (PES) scan was performed to study the passing of [C₄mim]⁺ and [PF₆]⁻ through different ZIF-8 apertures. Secondly, several electronic analysis methods, including frontier orbital energies, natural charge analysis, and energy decomposition analysis, were used to investigate the interactions between [C₄mim]⁺/[PF₆]⁻ and ZIF-8. Lastly, a series of IL-encapsulated ZIF-8 with different loading amounts were investigated to obtain stabilities of ionic pairs in the SOD cage of ZIF-8.

2. Results

First of all, we performed the DFT-based PES scan calculations. In general, small energy barriers would indicate a smooth diffusion of ILs inside the pore structures of MOFs. The PES scan of the [C₄mim]⁺ and [PF₆]⁻ passing through the pristine ZIF-8 aperture (with a pore size of 3.4 Å) is shown in Figure 1a. It is shown that both [C₄mim]⁺ and [PF₆]⁻ almost cannot pass through the pore because of their large energy barriers. For [C₄mim]⁺, an energy barrier of 135. 5 kcal·mol⁻¹ was obtained, and the corresponding structure is the one when the COM of imidazole almost perpendicularly overlaps with the COM of the pristine ZIF-8 aperture. For [PF₆]⁻, three cases have been considered (see Materials and Methods (Section 3)). The lowest energy barrier, which is 69.17 kcal·mol⁻¹, was found for case III, in which three F atoms pointed to the plane of the ZIF-8 apertures. This finding could be explained by their molecular diameters, in which case I has a diameter of 3.29 Å, case II has a diameter of 3.20 Å, and case III has a diameter of 2.94 Å (Figure S1). Meanwhile, as shown in Figure S2, the results of the energy decomposition analysis show that the repulsive interactions between the [PF₆]⁻ and ZIF-8 are the smallest when [PF₆]⁻ adopted case III. After obtaining the configurations with the highest energy barriers, we then refined the PES with the relaxed scan (see Figure 1b), and we found that both [C₄mim]⁺ and [PF₆]⁻ still have relatively high energy barriers, being 39.87 kcal·mol⁻¹ and 16.61 kcal·mol⁻¹, respectively, which still prohibited [C₄mim]⁺ and [PF₆]⁻ from diffusing into the SOD cage of ZIF-8 at finite temperatures, e.g., room temperature. In detail, only one energy barrier was found when [PF₆]⁻ moved from the left to the right side of the pore, and the highest energy point corresponds to the structure of which the P atom is located at the plane of the pristine ZIF-8, with three F atoms located at each side (see Figure S3a). The average distance between the H atoms of ZIF-8 and the F atoms of the [PF₆]⁻ is computed to be 2.52 Å. While for the [C₄mim]⁺, two energy barriers were found when [C₄mim]⁺ moved from the left to the right of the pristine ZIF-8 aperture. The first one, 14.60 kcal·mol⁻¹, corresponds to the passing of the alkyl chain, in which the average distance between the alkyl chain and the ZIF-8 aperture is 2.61 Å. The second one, 39.87 kcal·mol⁻¹, corresponds to the passing of the imidazole ring, in which the average distance between the imidazole ring and the ZIF-8 aperture is only 1.66 Å (Figure S3b); hence, such small distances result in large steric hindrances or strong repulsion interactions.

Afterwards, the frontier orbitals, including the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) were plotted in Figure 2. In general, the results indicate that there are no significant orbital interactions between the [C₄mim]⁺/[PF₆]⁻ and the pristine ZIF-8 aperture within the frontier orbital scheme. For example, both the LUMOs and HOMOs of [C₄mim]⁺-ZIF-8 and [PF₆]⁻-ZIF-8 are almost the same compared to those of the free pristine ZIF-8 aperture; in other words, the frontier orbitals mainly consist of the contributions from the ZIF-8 aperture. The HOMO–LUMO gaps of [C₄mim]⁺-ZIF-8 and [PF₆]⁻-ZIF-8 are also identical to those of the free pristine ZIF-8 aperture, e.g., 0.64 eV versus 0.65 eV. Moreover, the natural population analysis (NPA) was performed, and the natural charges as summarized in Table S1. The values show that a certain number of charges (ca. 0.08 e) have been transferred between [C₄mim]⁺ and [PF₆]⁻ and the pristine ZIF-8 apertures.

To have a deep understanding of the interactions between [C₄mim]⁺/[PF₆]⁻ and the pristine ZIF-8 aperture, the GKS-EDA energy decomposition analysis has been performed based on the [PF₆]⁻-ZIF-8(3.4 Å) and [C₄mim]⁺-ZIF-8(3.4 Å) structures and the results are plotted in Figure 3. Generally, the interaction energies are decomposed according to the following equations:

${\Delta E}^{total}={\Delta E}^{ex}+{\Delta E}^{rep}+{\Delta E}^{corr}+{\Delta E}^{ele}+{\Delta E}^{pol}+{\Delta E}^{disp},$

where ${\Delta E}^{ex}$, ${\Delta E}^{rep}$, ${\Delta E}^{corr}$, ${\Delta E}^{ele}$, ${\Delta E}^{pol}$, and ${\Delta E}^{disp}$ represent exchange, repulsion, correlation, electrostatic, polarization, and dispersion components, respectively. As shown in Figure 3, the interaction energies are dominated by three terms, including the repulsion, exchange, and electrostatic energies, of which the repulsions are the main reason for the high energy barriers shown in Figure 1. The rest of the components are relatively small (e.g., less than 28 kcal·mol⁻¹ in absolute values, and with a total ratio being less than 6.6%). Here, we take [C₄mim]⁺-ZIF-8(3.4 Å) as an example and discuss those three terms individually: (1) The repulsion component. This term is responsible for steric repulsion, and it consists of the destabilizing interactions between occupied orbitals of the fragments. Usually, this term is repulsive (with a positive computed value). Here, a value of 304.27 kcal·mol⁻¹ was obtained, which can be attributed to the large steric repulsion between the imidazole ring of [C₄mim]⁺ and the pristine ZIF-8 aperture. As discussed in the PES scan section, the pristine ZIF-8 aperture has a pore size of only 3.4 Å; hence, this leaves the imidazole ring few spaces to pass through. (2) The exchange component. This term essentially is related to electrons with the same spin, exchanging their positions in degenerate orbitals to increase the stability of electronic structure states. Here, we found a value of −162.2 kcal·mol⁻¹, in the case of [C₄mim]⁺-ZIF-8. (3) The electrostatic component. This term is the energy between the unperturbed charge distributions of the prepared fragments, which is usually attractive. As shown in Figure 3, the computed electrostatic interaction is −50 kcal·mol⁻¹, which is much smaller compared to the repulsion and exchange interactions. This finding is somehow consistent with the NPA charge analysis, in which the pristine ZIF-8 is slightly charged (ca. 0.08 e) via the charge transfer from [C₄mim]⁺.

Our previous MD simulations showed that the pristine ZIF-8 aperture could be distorted because of the thermal oscillations [36]. When the simulations were performed at a temperature of 300 K, an average of 6.8 Å movements was obtained for the atoms defined by the size and shape of the pores. In other words, we could conclude that the ZIF-8 structure is very flexible against the temperatures, and the pore size might be modified by varying the synthesis temperatures. Following that study, here, we computed the PES scan of [C₄mim]⁺ and [PF₆]⁻ passing through the ZIF-8 aperture with the pore size of 4.1 Å, 4.2 Å, 4.7 Å, and 6.9 Å, respectively, (the pore structures are depicted in Materials and Methods (Section 3)), and the results are summarized in Figure 4. Generally, the results show that the [PF₆]⁻ is able to freely pass through these pores with energy barriers under zero. However, [C₄mim]⁺ can only pass through the pore with a diameter of 6.9 Å. Together with our previous room-temperature MD simulations, we conclude that a higher temperature might be needed to obtain a larger swing effect and to encapsulate [C₄mim][PF₆] into the ZIF-8 SOD cage. These findings are rather qualitatively consistent with a reported experiment study, in which the authors employed a two-step adsorption/infiltration method to incorporate [C₄mim][PF₆] into ZIF-8 [37]. In that work, [C₄mim][PF₆] were firstly adsorbed on the outer surface of ZIF-8, then through heat treatment at 105 °C, [C₄mim][PF₆] molecules infiltrated into the SOD cage. Afterwards, the temperature was cooled down to room temperature, and then the [C₄mim][PF₆]-encapsulated ZIF-8 composite materials were obtained.

Lastly, to study how many pairs of [C₄mim][PF₆] could be stable inside the ZIF-8 SOD cage, and to examine the impact of their contents on the structures of the ZIF-8 apertures, the structural optimizations were performed by adding one, two, and three pairs of [PF₆][C₄mim] in the pristine ZIF-8 SOD cage. As shown in Figure 5a, it is shown that the interaction energies between [C₄mim][PF₆] and ZIF-8 were strongest when two pairs of [C₄mim][PF₆] were added, with a computed value of −44.81 kcal·mol⁻¹, whereas the interactions energies of one and three pairs were computed to be −35.69 kcal·mol⁻¹ and −40.92 kcal·mol⁻¹. As such, we conclude that most of the ZIF-8 SOD cages could accommodate two pairs of [C₄mim][PF₆], and some of them are even able to accommodate three pairs. Overall, a value slightly larger than 2 would be obtained for the pairs of [C₄mim][PF₆] per SOD cage. These findings are comparable to the experimental results obtained by Ban et al. [18], in which an average of 1.4 pairs of [C₄mim][TF₂N] was identified per ZIF-8 SOD cage. The smaller value of 1.4 versus 2.0 could be attributed to the fact that the anion [TF₂N]⁻ is relatively larger than [PF₆]⁻ in space size. Compared to the experimental values [38], we found that the volume of the ZIF-8 SOD cage was increased when [C₄mim][PF₆] was encapsulated, of which the smallest value was obtained (0.9%) for the case with two pairs. Moreover, we also observed certain confinement effects from the results shown in Figure 5c, e.g., the structure of one pair of [C₄mim][PF₆] is very similar compared to the free pair (e.g., gas phase calculations), while the two and three pairs show structures closer to the condensed phase [39].

3. Materials and Methods

The crystal structure of the ZIF-8 was obtained from the Cambridge Crystallographic Data Centre (CCDC) Crystal Library [24], and the simplified SOD cage presentation is depicted in Figure 6a. We choose the five typical pore apertures of the ZIF-8 structure as the result of the rotation of the zinc–imidazole–zinc and methyl functional groups, as described by our previous work, [35], including “Pristine (3.4 Å)”, “Closed” (4.1 Å), “Semi-open (4.2 Å)”, “Closed (4.7 Å)”, and “Open” (6.9 Å), respectively, which are shown in Figure 6b–f.

As demonstrated in Figure 7a, the [PF₆]⁻ and [C₄mim]⁺ were firstly placed 10 Å away from the aperture of the ZIF-8; then, gradually (e.g., with a step size of 0.1 Å), penetrate passed through the center of the aperture (the center of mass, COM). In principle, [PF₆]⁻ and [C₄mim]⁺ might have many orientations when passing through the five pore structures as shown in Figure 6. Here, we considered only the orientations that have possible low energy barriers. To summarize Figure 7b, three different orientations of [PF₆]⁻ have been studied, denoted as case I, II, and III, with one, two, and three fluorine atoms (F) towards the aperture, respectively. For the [C₄mim]⁺, one case was chosen towards the aperture, of which the imidazole ring is almost vertical to the aperture.

With the above-described models, rigid PES scans with single-point (SP) energy calculations were performed, where all atoms in the aperture of ZIF-8 and [PF₆]⁻/[C₄mim]⁺ were kept frozen. Afterward, we performed the relaxed PES scan within the regions near the highest energy barrier positions from the rigid PES scans (e.g., ca. −2 Å to 2 Å). In these cases, the atoms in the aperture of ZIF-8 were kept frozen, while [PF₆]⁻/[C₄mim]⁺ were fully optimized. Both PES scans were carried out with the Gaussian16 software package [40]. In the next discussion of both rigid and relaxed PES scans, we mainly focus on the largest energy barriers, which are the energy differences between the highest point on the PES and the free IL plus ZIF-8 structures. Subsequently, the orbital energies and charge transfer were computed with natural population analysis (NPA) [41] using the Multiwfn v.3.7 [42] software package and visualized using the VESTA v.3.5.8 [43] software package. To quantitatively study the interactions between [PF₆]⁻/[C₄mim]⁺ and ZIF-8, the energy decomposition analysis was performed based on the generalized Kohn–Sham energy density analysis (GKS-EDA) method [44], which is implemented in the GAMESS program package [45]. All above-mentioned DFT calculations were performed at the B3LYP level of theory, with a basis set of 6-311+G* [46], together with the damping scheme of Becke and Johnson, which denoted as DFT-D3 (BJ) [47,48].

Lastly, to study how many pairs are stable inside the SOD cage of ZIF-8, the periodic structures of ILs@ZIF-8 were examined by the BLYP functional with the dispersion correction scheme (DFT-D3). The geometry optimizations were performed to account for both changes in atomic positions and lattice dimensions, which were updated with the efficient L-BFGS algorithm. Within the DFT calculations, the DZVP-MOLOPT-SR-GTH basis set was adopted for Zn atoms, while the other atoms used the TZVP-MOLOPT -GTH basis set. The energy cutoff was set to 400 Ry, and the energy convergence for the self-consistent field (SCF) calculation was set to 1 × 10⁻⁵ Hartree. All those periodic calculations were performed by employing the CP2K’s Quickstep module [49].

4. Conclusions

In summary, the DFT calculations and several electronic structure analyses have been performed to study the diffusion behavior of the ILs, [C₄mim][PF₆], through the different aperture configurations of ZIF-8, and their stabilities inside the ZIF-8 cage. The results indicate that the original aperture configuration (with a 3.4 Å pore size) eventually prohibits the diffusion of [C₄mim][PF₆] into the SOD cage of ZIF-8, as a minimum energy barrier of 39.87 kcal·mol⁻¹ had been identified, which is mainly due to steric hindrance from the imidazole ring of the [C₄mim]⁺. The energy decomposition analysis based on the GKS-EDA method revealed that the large repulsion component is the main reason for the high energy barriers when [C₄mim]⁺ and [PF₆]⁻ pass through the pristine ZIF-8 aperture. Nevertheless, we found that the diffusion properties could be largely enhanced by modifying the ZIF-8 apertures, e.g., the pore size, and it is proved that this could be achieved via thermal contributions (e.g., employing a higher temperature during the synthesis procedure). Moreover, we found that two pairs of [C₄mim][PF₆] per SOD cage is the stable state, according to interaction energies and volume changes. Certain confinement effects were obtained as well, in which the structure of one pair of [C₄mim][PF₆] is similar compared to the free pair, while the two and three pairs show structures closer to the condensed phase.

Author Contributions

Investigation, L.L. (Longlong Liu) and K.J.; Data Curation, L.L. (Longlong Liu) and K.J.; Writing—Original Draft, L.L. (Longlong Liu) and K.J.; Writing—Review and Editing, Q.C. and L.L. (Lei Liu); Supervision, Q.C. and L.L. (Lei Liu). All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Materials. Further inquiries can be directed to the corresponding authors.

Conflicts of Interest

The authors declare that they have no known competing financial interests.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures

Figure 1. The potential energy surface curves for the anion [PF6]− and cation [C4mim]+ passing through ZIF-8 structures: (a) the rigid PES scan and (b) the relaxed PES scan.

Figure 2. Frontier orbitals and their energies, with an isosurface value of 0.05 e/Å3.

Figure 3. Energy decomposition analysis calculated with the GKS-EDA method at the B3LYP-D3(BJ)/6-311+G* level of theory for the [PF6]−/[C4mim]+-ZIF-8(3.4) systems.

Figure 4. Potential energy surface curves for [C4mim]+ (a) and [PF6]− (b) passing through different ZIF-8 apertures.

Figure 5. Interaction energies (a), volume changes (b), and geometric structures (c) of the ZIF-8 SOD cage containing one, two, and three pairs of [PF6][C4mim], respectively.

Figure 6. (a) The SOD cage presentation of the ZIF-8, and its typical (pristine) pore structure is highlighted in red color; (b–f) five different pore structures of ZIF-8 with their diameters (ref. [33]). Color legend: C, black; N, blue; Zn, light grey; H, white.

Figure 7. (a) Schematic representation of [PF6]−/[C4mim]+ passing through ZIF-8 apertures. (b) The geometric structures of [PF6]−/[C4mim]+, and their orientations when passing through the ZIF-8 aperture.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/molecules29081697/s1: Figure S1: Three different ways and their diameters of the [PF₆] pathing through the ZIF-8 pore; Figure S2: Interaction energy components calculated using the GKS-EDA method at the B3LYP-D3(BJ)/6-311+G* level for the different [PF₆]⁻ ways towards the ZIF-8 aperture; Figure S3: The structures for the points with the highest interaction energies in the cases of [PF₆]⁻/[C₄mim]⁺ and ZIF-8(3.4), respectively, and the average distance refers to the closest distances between the fluorine atom on the [PF₆]⁻, or the hydrogen atoms on the [C₄mim]⁺ and the hydrogen atom on the aperture; and Table S1: NPA charge analysis of [PF₆]⁻/[C₄mim]⁺ and ZIF-8(3.4) at theoretical levels of B3LYP-D3(BJ)/6-311+G*.(values are in e).

References

1. Li, X.; Chen, K.; Guo, R.; Wei, Z. Ionic Liquids Functionalized MOFs for Adsorption. Chem. Rev.; 2023; 123, pp. 10432-10467. [DOI: https://dx.doi.org/10.1021/acs.chemrev.3c00248] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/37450853]

2. Durak, O.; Zeeshan, M.; Habib, N.; Gulbalkan, H.C.; Alsuhile, A.A.A.M.; Caglayan, H.P.; Kurtoğlu-Öztulum, S.F.; Zhao, Y.; Haslak, Z.P.; Uzun, A. et al. Composites of porous materials with ionic liquids: Synthesis, characterization, applications, and beyond. Microporous Mesoporous Mater.; 2022; 332, 111703. [DOI: https://dx.doi.org/10.1016/j.micromeso.2022.111703]

3. Wu, K.; Miao, X.; Zhao, H.; Liu, S.; Fei, T.; Zhang, T. Selective Encapsulation of Ionic Liquids in UiO-66-NH₂ Nanopores for Enhanced Humidity Sensing. ACS Appl. Nano Mater.; 2023; 6, pp. 9050-9058. [DOI: https://dx.doi.org/10.1021/acsanm.3c01727]

4. Qian, Y.; Zhang, F.; Pang, H. A Review of MOFs and Their Composites-Based Photocatalysts: Synthesis and Applications. Adv. Funct. Mater.; 2021; 31, 2104231. [DOI: https://dx.doi.org/10.1002/adfm.202104231]

5. Friess, K.; Izák, P.; Kárászová, M.; Pasichnyk, M.; Lanč, M.; Nikolaeva, D.; Luis, P.; Jansen, J.C. A Review on Ionic Liquid Gas Separation Membranes. Membranes; 2021; 11, 97. [DOI: https://dx.doi.org/10.3390/membranes11020097] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33573138]

6. Fernandez, E.; Saiz, P.G.; Peřinka, N.; Wuttke, S.; Fernández de Luis, R. Printed Capacitive Sensors Based on Ionic Liquid/Metal-Organic Framework Composites for Volatile Organic Compounds Detection. Adv. Funct. Mater.; 2021; 31, 2010703. [DOI: https://dx.doi.org/10.1002/adfm.202010703]

7. Tuffnell, J.M.; Morzy, J.K.; Kelly, N.D.; Tan, R.; Song, Q.; Ducati, C.; Bennett, T.D.; Dutton, S.E. Comparison of the ionic conductivity properties of microporous and mesoporous MOFs infiltrated with a Na-ion containing IL mixture. Dalton Trans.; 2020; 49, pp. 15914-15924. [DOI: https://dx.doi.org/10.1039/D0DT02576C]

8. Zeeshan, M.; Nozari, V.; Keskin, S.; Uzun, A. Structural Factors Determining Thermal Stability Limits of Ionic Liquid/MOF Composites: Imidazolium Ionic Liquids Combined with CuBTC and ZIF-8. Ind. Eng. Chem. Res.; 2019; 58, pp. 14124-14138. [DOI: https://dx.doi.org/10.1021/acs.iecr.9b02415]

9. Yoshida, Y.; Fujie, K.; Lim, D.W.; Ikeda, R.; Kitagawa, H. Superionic Conduction over a Wide Temperature Range in a Metal–Organic Framework Impregnated with Ionic Liquids. Angew. Chem. Int. Ed.; 2019; 58, pp. 10909-10913. [DOI: https://dx.doi.org/10.1002/anie.201903980] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31140203]

10. Li, H.; Tuo, L.; Yang, K.; Jeong, H.-K.; Dai, Y.; He, G.; Zhao, W. Simultaneous enhancement of mechanical properties and CO₂ selectivity of ZIF-8 mixed matrix membranes: Interfacial toughening effect of ionic liquid. J. Membr. Sci.; 2016; 511, pp. 130-142. [DOI: https://dx.doi.org/10.1016/j.memsci.2016.03.050]

11. Khan, M.S.; Li, Y.; Yang, L.; Yan, Z.C.; Li, D.-S.; Qiu, J.; Xu, X.; Yang, H.Y. Improving capacitive deionization performance through tailored iodine-loaded ZIF-8 composites. Desalination; 2024; 579, 117486. [DOI: https://dx.doi.org/10.1016/j.desal.2024.117486]

12. Khan, M.S.; Leong, Z.Y.; Li, D.-S.; Qiu, J.; Xu, X.; Yang, H.Y. A mini review on metal–organic framework-based electrode materials for capacitive deionization. Nanoscale; 2023; 15, pp. 15929-15949. [DOI: https://dx.doi.org/10.1039/D3NR03993E] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/37772477]

13. Khan, M.S.; Li, Y.; Li, D.-S.; Qiu, J.; Xu, X.; Yang, H.Y. A review of metal–organic framework (MOF) materials as an effective photocatalyst for degradation of organic pollutants. Nanoscale Adv.; 2023; 5, pp. 6318-6348. [DOI: https://dx.doi.org/10.1039/D3NA00627A] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/38045530]

14. Chen, W.; Zhang, Z.; Yang, C.; Liu, J.; Shen, H.; Yang, K.; Wang, Z. PIM-based mixed-matrix membranes containing MOF-801/ionic liquid nanocomposites for enhanced CO₂ separation performance. J. Membr. Sci.; 2021; 636, 119581. [DOI: https://dx.doi.org/10.1016/j.memsci.2021.119581]

15. Zeeshan, M.; Keskin, S.; Uzun, A. Enhancing CO₂/CH₄ and CO₂/N₂ separation performances of ZIF-8 by post-synthesis modification with [BMIM][SCN]. Polyhedron; 2018; 155, pp. 485-492. [DOI: https://dx.doi.org/10.1016/j.poly.2018.08.073]

16. Krokidas, P.; Moncho, S.; Brothers, E.N.; Castier, M.; Economou, I.G. Tailoring the gas separation efficiency of metal organic framework ZIF-8 through metal substitution: A computational study. Phys. Chem. Chem. Phys.; 2018; 20, pp. 4879-4892. [DOI: https://dx.doi.org/10.1039/C7CP08456K] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29384175]

17. Gulbalkan, H.C.; Haslak, Z.P.; Altintas, C.; Uzun, A.; Keskin, S. Assessing CH₄/N₂ separation potential of MOFs, COFs, IL/MOF, MOF/Polymer, and COF/Polymer composites. Chem. Eng. J.; 2022; 428, 131239. [DOI: https://dx.doi.org/10.1016/j.cej.2021.131239]

18. Ban, Y.; Li, Z.; Li, Y.; Peng, Y.; Jin, H.; Jiao, W.; Guo, A.; Wang, P.; Yang, Q.; Zhong, C. et al. Confinement of Ionic Liquids in Nanocages: Tailoring the Molecular Sieving Properties of ZIF-8 for Membrane-Based CO₂ Capture. Angew. Chem. Int. Ed. Engl.; 2015; 54, pp. 15483-15487. [DOI: https://dx.doi.org/10.1002/anie.201505508] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26515558]

19. Kinik, F.P.; Altintas, C.; Balci, V.; Koyuturk, B.; Uzun, A.; Keskin, S. [BMIM][PF6] Incorporation Doubles CO₂ Selectivity of ZIF-8: Elucidation of Interactions and Their Consequences on Performance. ACS Appl. Mater. Interfaces; 2016; 8, pp. 30992-31005. [DOI: https://dx.doi.org/10.1021/acsami.6b11087] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27783899]

20. Koyuturk, B.; Altintas, C.; Kinik, F.P.; Keskin, S.; Uzun, A. Improving Gas Separation Performance of ZIF-8 by [BMIM][BF4] Incorporation: Interactions and Their Consequences on Performance. J. Phys. Chem. C; 2017; 121, pp. 10370-10381. [DOI: https://dx.doi.org/10.1021/acs.jpcc.7b00848]

21. Ali, S.A.; Khan, A.U.; Mulk, W.U.; Khan, H.; Nasir Shah, S.; Zahid, A.; Habib, K.; Shah, M.U.H.; Othman, M.H.D.; Rahman, S. An Ongoing Futuristic Career of Metal–Organic Frameworks and Ionic Liquids, A Magical Gateway to Capture CO₂; A Critical Review. Energy Fuels; 2023; 37, pp. 15394-15428. [DOI: https://dx.doi.org/10.1021/acs.energyfuels.3c02377]

22. Dong, K.; Liu, X.; Dong, H.; Zhang, X.; Zhang, S. Multiscale Studies on Ionic Liquids. Chem. Rev.; 2017; 117, pp. 6636-6695. [DOI: https://dx.doi.org/10.1021/acs.chemrev.6b00776] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28488441]

23. Chen, J.; Dong, K.; Liu, L.; Zhang, X.; Zhang, S. Anti-electrostatic hydrogen bonding between anions of ionic liquids: A density functional theory study. Phys. Chem. Chem. Phys.; 2021; 23, pp. 7426-7433. [DOI: https://dx.doi.org/10.1039/D0CP06718K] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33876102]

24. Kwon, H.T.; Jeong, H.-K.; Lee, A.S.; An, H.S.; Lee, J.S. Heteroepitaxially Grown Zeolitic Imidazolate Framework Membranes with Unprecedented Propylene/Propane Separation Performances. J. Am. Chem. Soc.; 2015; 137, pp. 12304-12311. [DOI: https://dx.doi.org/10.1021/jacs.5b06730] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26364888]

25. Thomas, A.; Prakash, M. Tuning the CO₂ adsorption by the selection of suitable ionic liquids at ZIF-8 confinement: A DFT study. Appl. Surf. Sci.; 2019; 491, pp. 633-639. [DOI: https://dx.doi.org/10.1016/j.apsusc.2019.06.130]

26. Mohamed, A.M.O.; Moncho, S.; Krokidas, P.; Kakosimos, K.; Brothers, E.N.; Economou, I.G. Computational investigation of the performance of ZIF-8 with encapsulated ionic liquids towards CO₂ capture. Mol. Phys.; 2019; 117, pp. 3791-3805. [DOI: https://dx.doi.org/10.1080/00268976.2019.1666170]

27. Gao, W.; Zheng, W.; Sun, W.; Zhao, L. Understanding the Effective Capture of H₂S/CO₂ from Natural Gas Using Ionic Liquid@MOF Composites. J. Phys. Chem. C; 2022; 126, pp. 19872-19882. [DOI: https://dx.doi.org/10.1021/acs.jpcc.2c05667]

28. Thomas, A.; Ahamed, R.; Prakash, M. Selection of a suitable ZIF-8/ionic liquid (IL) based composite for selective CO₂ capture: The role of anions at the interface. RSC Adv.; 2020; 10, pp. 39160-39170. [DOI: https://dx.doi.org/10.1039/D0RA07927H] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35518407]

29. Kavak, S.; Polat, H.M.; Kulak, H.; Keskin, S.; Uzun, A. MIL-53(Al) as a Versatile Platform for Ionic-Liquid/MOF Composites to Enhance CO₂ Selectivity over CH₄ and N₂. Chem.–Asian J.; 2019; 14, pp. 3655-3667. [DOI: https://dx.doi.org/10.1002/asia.201900634] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31339661]

30. Zhang, Z.; Jia, X.; Sun, Y.; Guo, X.; Huang, H.; Zhong, C. Pore engineering of ZIF-8 with ionic liquids for membrane-based CO₂ separation: Bearing functional group effect. Green Chem. Eng.; 2021; 2, pp. 104-110. [DOI: https://dx.doi.org/10.1016/j.gce.2020.10.007]

31. Fairen-Jimenez, D.; Moggach, S.A.; Wharmby, M.T.; Wright, P.A.; Parsons, S.; Düren, T. Opening the Gate: Framework Flexibility in ZIF-8 Explored by Experiments and Simulations. J. Am. Chem. Soc.; 2011; 133, pp. 8900-8902. [DOI: https://dx.doi.org/10.1021/ja202154j] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21553843]

32. Zhang, C.; Lively, R.P.; Zhang, K.; Johnson, J.R.; Karvan, O.; Koros, W.J. Unexpected Molecular Sieving Properties of Zeolitic Imidazolate Framework-8. J. Phys. Chem. Lett.; 2012; 3, pp. 2130-2134. [DOI: https://dx.doi.org/10.1021/jz300855a] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26295759]

33. Zeeshan, M.; Nozari, V.; Yagci, M.B.; Isık, T.; Unal, U.; Ortalan, V.; Keskin, S.; Uzun, A. Core–Shell Type Ionic Liquid/Metal Organic Framework Composite: An Exceptionally High CO₂/CH₄ Selectivity. J. Am. Chem. Soc.; 2018; 140, pp. 10113-10116. [DOI: https://dx.doi.org/10.1021/jacs.8b05802]

34. Chang, Y.; Wang, L.; Jiang, Z.; Zhang, R.; Zhu, H.; Zhang, D.; Zhu, J.; Kong, X.; Huang, H. Zwitterionic metal–organic framework with highly dispersed ionic liquid for enhancing CO₂ capture. Sep. Purif. Technol.; 2023; 326, 124828. [DOI: https://dx.doi.org/10.1016/j.seppur.2023.124828]

35. Yu, T.; Cai, Q.; Lian, G.; Bai, Y.; Zhang, X.; Zhang, X.; Liu, L.; Zhang, S. Mechanisms behind high CO₂/CH₄ selectivity using ZIF-8 metal organic frameworks with encapsulated ionic liquids: A computational study. Chem. Eng. J.; 2021; 419, 129638. [DOI: https://dx.doi.org/10.1016/j.cej.2021.129638]

36. Yu, T.; Cai, Q.; Lian, G.; Liu, L. Molecular dynamics studies on separation of CO₂/CH₄ by the ionic liquids encapsulated ZIF-8. J. Membr. Sci.; 2022; 644, 120117. [DOI: https://dx.doi.org/10.1016/j.memsci.2021.120117]

37. Guo, Z.; Zheng, W.; Yan, X.; Dai, Y.; Ruan, X.; Yang, X.; Li, X.; Zhang, N.; He, G. Ionic liquid tuning nanocage size of MOFs through a two-step adsorption/infiltration strategy for enhanced gas screening of mixed-matrix membranes. J. Membr. Sci.; 2020; 605, 118101. [DOI: https://dx.doi.org/10.1016/j.memsci.2020.118101]

38. Park, K.S.; Ni, Z.; Côté, A.P.; Choi, J.Y.; Huang, R.; Uribe-Romo, F.J.; Chae, H.K.; O’Keeffe, M.; Yaghi, O.M. Exceptional chemical and thermal stability of zeolitic imidazolate frameworks. Proc. Natl. Acad. Sci. USA; 2006; 103, pp. 10186-10191. [DOI: https://dx.doi.org/10.1073/pnas.0602439103] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16798880]

39. Bühl, M.; Wipff, G. Insights into Uranyl Chemistry from Molecular Dynamics Simulations. ChemPhysChem; 2011; 12, pp. 3095-3105. [DOI: https://dx.doi.org/10.1002/cphc.201100458] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21960471]

40. Frisch, G.A.; Trucks, M.J.; Schlegel, G.W.; Scuseria, H.B.; Robb, G.E.; Cheeseman, M.A.; Scalmani, J.R.; Barone, G.; Petersson, V. Gaussian 16, Rev. C.01; Gaussian, Inc.: Wallingford, CT, USA, 2016.

41. Weinhold, F.; Landis, C.; Glendening, E.D. What is NBO analysis and how is it useful?. Int. Rev. Phys. Chem.; 2016; 35, pp. 399-440. [DOI: https://dx.doi.org/10.1080/0144235X.2016.1192262]

42. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem.; 2012; 33, pp. 580-592. [DOI: https://dx.doi.org/10.1002/jcc.22885] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22162017]

43. Momma, K.; Izumi, F. VESTA: A three-dimensional visualization system for electronic and structural analysis. J. Appl. Crystallogr.; 2008; 41, pp. 653-658. [DOI: https://dx.doi.org/10.1107/S0021889808012016]

44. Su, P.; Jiang, Z.; Chen, Z.; Wu, W. Energy decomposition scheme based on the generalized Kohn-Sham scheme. J. Phys. Chem. A; 2014; 118, pp. 2531-2542. [DOI: https://dx.doi.org/10.1021/jp500405s] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24611964]

45. Schmidt, M.W.; Baldridge, K.K.; Boatz, J.A.; Elbert, S.T.; Gordon, M.S.; Jensen, J.H.; Koseki, S.; Matsunaga, N.; Nguyen, K.A.; Su, S. et al. General atomic and molecular electronic structure system. J. Comput. Chem.; 1993; 14, pp. 1347-1363. [DOI: https://dx.doi.org/10.1002/jcc.540141112]

46. Stephens, P.J.; Devlin, F.J.; Chabalowski, C.F.; Frisch, M.J. Ab initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields. J. Phys. Chem.; 1994; 98, pp. 11623-11627. [DOI: https://dx.doi.org/10.1021/j100096a001]

47. Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem.; 2011; 32, pp. 1456-1465. [DOI: https://dx.doi.org/10.1002/jcc.21759] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21370243]

48. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys.; 2010; 132, 154104. [DOI: https://dx.doi.org/10.1063/1.3382344] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20423165]

49. Kühne, T.D.; Iannuzzi, M.; Del Ben, M.; Rybkin, V.V.; Seewald, P.; Stein, F.; Laino, T.; Khaliullin, R.Z.; Schütt, O.; Schiffmann, F. et al. CP2K: An electronic structure and molecular dynamics software package—Quickstep: Efficient and accurate electronic structure calculations. J. Chem. Phys.; 2020; 152, 194103. [DOI: https://dx.doi.org/10.1063/5.0007045] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33687235]