Hari P.
Paudel
*,
Wei
Shi
,
David
Hopkinson
,
Janice A.
Steckel
and
Yuhua
Duan
*
National Energy Technology Laboratory, United States Department of Energy, Pittsburgh, Pennsylvania 15236, USA. E-mail: Hari.Paudel@netl.doe.gov; yuhua.duan@netl.doe.gov
First published on 13th January 2021
Understanding of zeolitic imidazolate framework-8 (ZIF-8) interaction with different gas molecules is crucial when ZIF-8 is used in gas separation. Computational studies based on density functional theory (DFT) can be used to investigate gas interactions and diffusion mechanisms that can be directly correlated with experimental observations. Here we present our studies based on DFT calculations on CO2 and CH4 gas adsorption and diffusion in the bulk of ZIF-8. We evaluate the structural and electronic properties of bulk ZIF-8, and determine the most stable adsorption sites and the corresponding diffusion barriers for CO2 and CH4 molecules. Our calculations incorporate long-range dispersion interactions to describe the weak interactions between adsorbate molecules and the framework. We analyze the adsorption and diffusion properties in relation with the material's volume expansion. We find that the CO2 and CH4 adsorption energies at the most stable adsorption sites are 5.01 and 4.47 kcal mol−1, respectively. The diffusivity of CO2 is found to be about two times that of CH4. Our calculated diffusion coefficients were found to have the same order of magnitude with the experimental results. Furthermore, our calculations indicate that CO2 and CH4 diffusivities in fixed ZIF-8 (all ZIF-8 framework atoms are fixed) are 5–9 times lower than the corresponding diffusion values in flexible ZIF-8 (all the framework atoms are allowed to move).
There have been significant efforts devoted to characterizing zeolite imidazolate frameworks (ZIFs) in order to efficiently and economically use them in membrane technology to separate gas components from mixtures. ZIFs have selective adsorption properties that enable them to discriminate on a particular type of gas molecule to be adsorbed more than others.13–16 When compared to traditional chemical sorbent materials (mostly amines), ZIFs have (a) high CO2 loading capacity, (b) lower energy consumption, (c) high thermal and chemical stability, and (d) lower corrosion rates.17–20 Due to these properties and the presence of permanent porosity, ZIFs outperform typical metal organic frameworks (MOFs) for gas separation. One of the most studied systems within the family of zeolitic imidazolium frameworks is ZIF-8. This system has a sodalite topology and a Zn(mim)2 stoichiometric representation, where mim represents 2-methylimidazolate.21 Within the same topology, ZIFs can be synthesized with different metal centers such as Zn, Co, Cu and Fe that are tetrahedrally coordinated and linked by imidazolate ligands.22,23 Importantly, the strength of the bonding between the metal ions and the ligands is thought to provide a high chemical stability.21,24 The structure of ZIF-8 consists of a central nanopore with a 11.6 Å diameter which is surrounded by six and four window pores of 3.4 and 0.8 Å diameters, respectively.25,26 The size of 3.4 Å is comparable to the kinetic diameters of CO2 (3.3 Å) and CH4 (3.8 Å) molecules, making ZIF-8 an ideal material for gas mixture separations. A total of 17 polyforms of Zn(mim)2 have been synthesized to date but many of them do not meet the optimal requirements needed for gas separation. In particular, they either lack the optimal pore size, collapse as they interact with guest molecules or are mechanically weak.27
In addition to gas separation applications, ZIF-8 has a number of other applications such as removal of oil from a water surface,28 on-demand drug delivery,29 heterogenous catalysts,30 chemical vapors and gas sensing and detection.31 ZIF-8 can be integrated with sensitive transducers to achieve high sensitivity in the microwave wavelength regime for CO2 and CH4 detection.32 A review on ZIF-8 with areas of applications is given in ref. 33.
By using high resolution transmission electron microscopy (HRTEM), Zhu et al. revealed experimentally the surface and interfacial structures of ZIF-8.14 In their experiment, an armchair type surface with a (110) index was found to be stable with Zn2+ ions capped by 2-mim as viewed along the [110] direction. This result was shown to be supported by theoretical calculations obtained by using a first principles density functional theory-tight binding (DFT-TB) approach.14 While Zhu et al. work provides an important step towards modelling the surface and interface of ZIF-8, an understanding of the adsorption and diffusion properties at the atomistic scale within the bulk, from the bulk to the surface, and MOF/polymer interfaces is still lacking. Semino et al. work on ZIF-8 interfaced with polymer of intrinsic microporosity-1 (PIM-1) provides a comprehensive theoretical model for the surface and interface while simulating a composite structure.34 By integrating experiment and theory, Hobday et al. revealed the most important adsorption sites of CH4, O2, N2 and Ar gas molecules within the framework of ZIF-8.20 While the later work provides gas adsorption properties at different pressures, the mechanism of diffusion and its changes upon loading of different gas molecules are still unexplained. Molecular properties can be altered greatly by applying an external electric field. Knebel et al. measured the gas permeances for H2, CO2 and CH4 under an applied electric field with a strength of 500 V mm−1.35 The authors performed a comparative study of single-gas permeances with and without the application of an in situ-applied field perpendicular to the MOF layer. Their results indicate that ZIF-8 undergoes a crystallographic structural transformation into polyforms under the applied field, altering the value of permeances. There are several experimental and theoretical studies performed to understand the relaxation in size of pores called “gate opening” to allow gas molecules with size larger than the equilibrium pore diameter of the framework to permeate smoothly through the openings. On the theoretical side, most of the quantitative reports on gate opening and diffusions of gas molecules rely upon molecular dynamics and Monte-Carlo simulations,16,36–40 and some upon DFT and DFTB especially for adsorption energy calculations.20,41,42 Garberoglio et al. have used DFTB to calculate diffusion of H2, CO2 and CH4 molecules in MOF materials with several hundred atoms in a unit cell.43 Recently, Verploegh et al. have introduced an accurate force field which is based on DFT parameterization for ZIFs.44 This work inherently considers the framework flexibility unlike in the cases where generic force fields are empirically tuned and fitted to sets of structural data obtained from the experiment.45–47 Fischer and Bell used DFT with periodic boundary conditions including Grimme's D2 dispersion to identify H2 and CO2 adsorption sites in ZIF-8 and other isostructural systems.48,49 The average interaction energies between H2 and the framework, and between CO2 and the framework were found to be −10.4 and −24.6 kJ mol−1. By combining ab initio MD with classical MD calculations, improved results could be achieved for framework flexibility while CO2 and CH4 diffuse in the pore.50 There are other similar studies on ZIF structures by using different flavors of DFT exchange–correlation functionals.51–53
While there are studies on experiments and theory-simulations that unravel adsorption and diffusion kinetics for gas molecules such as CO2 and CH4 in ZIF-8 materials at the atomistic-level, a deeper understanding of the diffusion of these gas molecules in relation to adsorption and framework expansion using fundamental concepts could help in improving the membrane design and optimal performance. Calculations of surface and interface properties of ZIF-8 relevant to gas mixture separation using theoretical modelling based entirely on quantum mechanics can be challenging due to the size of the framework unit cell. Nevertheless, it is possible to obtain a deeper insight into adsorption and diffusion kinetics such as diffusion pathways, energy profiles and energy hierarchy of adsorption sites for different gas molecules in the bulk structure from quantum mechanical calculations. Diffusion barrier profiles with and without atomic position relaxation allow us to quantitatively explain the orientation dependent interaction of the guest molecule with the framework and its structural changes at pore sites.
There are force field parameters available for ZIF-8 in the open literature and one can naturally think of using the MD/MC method over DFT for such a system.54 In ref. 55, DFT based force field parameters were developed for MOF materials including ZIF-8. Therefore, using the DFT method it is desirable to understand the fundamentals of adsorption and diffusion properties of ZIF-8. This opens a possibility of extending the DFT method to many other systems for which the development of a classical force field is typically a lengthy and daunting work. Therefore, in systems for which force field parameters generally lack, DFT will be an alternative method with relatively low computational cost because developing reasonable classical force field parameters is a time-consuming task. In this work, we show that the DFT method can predict the gas diffusivity and gas saturated loading in ZIF-8.
In this study, we perform first principles DFT calculation with Grimme D3 dispersion corrections to study the adsorption and diffusion properties of CO2 and CH4 in ZIF-8. Using the nudged elastic band (NEB) approach we calculate the energy profile for gas molecules while they diffuse in the bulk structure.56–58 We identify the hierarchy of adsorption sites for gas molecules interacting with the bulk framework atoms. The bond angle of CO2 remains unaltered in ZIF-8, which confirms that the interaction between the framework and guest molecule is mainly of the van der Waals type. The calculated maximum loading amount identified in this work can be compared with the experimental value42 obtained at the saturation limit at high pressure limit. The diffusion energy profiles we presented in this work are for single molecules of CO2 or CH4 under dynamic equilibrium in the DFT framework. CO2 and CH4 are linearly and spherically symmetrical molecules, whose center of mass coincide with the geometrical center of the molecule. There are several possible diffusion pathways for the guest molecule within the framework. We provided only the pathways with the minimum diffusion barriers. The response of the framework to the guest molecule is constrained as well as relaxed to see the impact of the change in framework–guest molecule interaction on the diffusion barriers. We observe that there is about 0.88 kcal mol−1 change in the barrier energy for CO2 and 1.25 kcal mol−1 for CH4 while relaxing both the framework and molecule and the molecule only. This is evidence of framework flexibility in porous materials like ZIF-8, where there are different geometrical pore sizes and not all of them allow for gas diffusion. The size description of pores is used in this work to identify the migration pathways while mapping out the images along a particular path. In fact, the probability of diffusing the molecule through a geometrical pore size of 0.8 Å is close to zero. We calculate the diffusivities corresponding to the obtained energy barriers. The order of magnitude of our calculated diffusion coefficients agrees with the experimental results measured at 298 K.42 We find that membranes made from ZIF-8 have good adsorption and diffusion selective properties for CO2 and CH4 molecules, making these membranes useful in gas separation processes.
EDFT–D = EKS–DFT + EDisp | (1) |
![]() | (2) |
The ZIF-8 structure obtained from the experiment was cleaned to remove disorder and relaxed without any symmetry constraints to allow full relaxation of both the lattice parameters and atomic positions. In order to calculate the diffusion energy profile, the nudged elastic band (NEB) approach was used to map the migration pathways of CO2 and CH4 molecules between two local minima.58 The NEB approach has been shown to predict energy barriers with high accuracy in many solid systems.64,65 The initial and final states were confirmed to be minima based on the lack of imaginary frequency modes.
![]() | (3) |
![]() | ||
Fig. 1 Optimized geometry of CO2 adsorption sites in ZIF 8. Fragments a–f which are labeled I, II, III, IV, V and VI, respectively, show six different representative sites for CO2 adsorption. |
We calculate the most favorable binding site of the CO2 molecule inside the framework unit cell (see the ESI†). The calculated value of −5.01 kcal mol−1 for the adsorption energy is comparable with the experimentally reported value for the isosteric heat of adsorption for the CO2 molecule at low coverage limit. In Fig. 1, we also show six different binding sites for CO2. For example, site I lies close to the metal center whereas site III is near the pore of size 0.8 Å. In site I, due to a strong repulsion from the framework's linker atoms we observe that the CO2 molecule is bent by an angle of 3.5° from its original linear structure, indicating that a slight charge transfer occurs between the framework atoms and the CO2 molecule. Detailed information on the adsorption of CO2 is provided in the ESI.†
To understand the change in adsorption energy because of the increase in coverage of the framework–molecule volume, we calculate the adsorption energy per CO2 for an increasing number of adsorbed CO2 molecules. The adsorption energy is calculated both by allowing a change in the volume of the unit cell and by holding the volume fixed. As the loading of molecules increases, the adsorption energy per CO2 molecule fluctuates slightly until 23 molecules are loaded in the pore. After reaching this loading, the adsorption energy per molecule rapidly decreases (in absolute value) as shown in Fig. 2a. Perez-Pellitero et al. reported an experimental loading saturation for CO2 of 8.5 mmol g−1 (∼23 CO2 per unit cell) at 303 K at a pressure of around 30 bar.66 The peak loading amount before the decrease in the DFT-calculated adsorption energy per molecule in Fig. 2a can be considered as an upper limit of the loading amount seen experimentally. After the dotted line in Fig. 2a, the adsorption energy falls in average with a scale of 0.38 kcal per mole per molecule added. Overall, the volume of the unit cell is found to increase noticeably only after the loading of 16 molecules as shown in Fig. 2b where we show the percentage volume change and adsorption energy as a function of the number of CO2 molecules in the cage. It is interesting to note from Fig. 2b that there is a significant change in the volume of the unit cell before the start of the rapid fall in adsorption energy. Our calculations show that the expansion in lattice parameters is uniform in all directions and the change in unit cell angles is negligibly small. In addition, we find that the change of volume by about more than 2% effects the loading amount noticeably as the adsorption energy per molecule rapidly falls after that.
As the loading of CH4 molecules increases in ZIF-8, the change in adsorption energy per molecule is found to be within in 0.1 kcal mol−1 until about 17 molecules (∼6.2 mmol g−1) are adsorbed by the framework for both cases of fixed and relaxed volume as shown in Fig. 4a. After this loading, the adsorption energy weakens with an average of 0.18 kcal mol−1 for each molecule added. This attenuation is more than two times steeper than the one observed in the case of CO2. Perez-Pellitero et al. measured the adsorption of CH4 at 303 K, and observed about 4 mmol g−1, but saturation was not achieved in their measurements.66 The application of sufficiently high pressure in certain directions in the experiment results in an anisotropic volume expansion, dilation along the applied force and contraction perpendicular to it. This results in a higher number of gas molecules coming closer to the framework's wall in the perpendicular direction. This results in a different loading amount from the situation of uniform expansion in all directions. Knebel et al.35 have demonstrated changes in lattice parameters upon electric field poling with a strength of 0.02 eV Å−1 followed by a change in phase transformation of ZIF-8 into different polyforms which subsequently reduced the separation selectivity for molecules like H2, CO2 and CH4. There is a difference of less than 0.1 kcal per mole in the adsorption energy at low loading when the volume of the unit cell is relaxed compared to when the unit cell volume is fixed. As we allow changes in the atomic positions and lattice parameters, for the adsorption of the first 5 CH4 molecules, the volume is found to decrease by 0.15% as shown in Fig. 4b. After the adsorption of 17 molecules, the average rate of increase of volume expansion is 13.6 Å3 per molecule for CH4 whereas for CO2 it was 15.5 Å3 per molecule. Therefore, at higher loading, the change in volume and adsorption energy as a function of the number of molecules is found to be larger in the case of CO2 as compared to the case of CH4.
As the number of molecules inside the cage increases, the interactions among the gas molecules and between the gas molecules and the framework increase significantly, and so the expansion in volume occurs.67 The orientation of CO2 molecules becomes important at higher loading. The adsorption energy changes faster with the number of molecules in the framework for CH4 than for CO2. This sets up the lower loading limit for CH4 than for CO2. As the loading increases, the available space within the cage is smaller for CH4 due to its larger size than for CO2 for an equal amount of loading. In such a case, there is a repulsion between two molecules, which is stronger for CH4 that effectively reduces the loading limit for CH4 than for CO2.
Under the NEB approach, the migration path of the molecule is relaxed in every step and is mapped out vectorially along the direction of the final image. The different vibrational modes of the framework and molecule may give inconsistent results due to the quasi-static responses of their motion under the NEB assumption. To avoid that while calculating the barrier energy, in addition to the entire material relaxation, we fix the MOF and relax the molecule only and vice versa. Fixing the framework and relaxing the molecule results in overall a rather coherent vibrational mode than the ones resulting from the quasi-static responses.
In Fig. 5 we show the initial, transition and final states of diffusing molecules of CO2 (Fig. 5a–c) and CH4 (Fig. 5d–f). For distinct visualization of diffusing molecules, the position of each image within the unit cell is identified with reference to the (111) and (110) Miller planes. The initial image of CO2 lies in the (110) Miller plane near the 6-MR window. The (110) plane bisects the 6-MR ring along the [100] direction as shown in Fig. 5a. The minimum energy diffusion pathway for CO2 lies along the [110] direction. Similarly, the initial state is located slightly off the [110] plane near the 6-MR window as shown in Fig. 5d. The minimum energy diffusion pathway for CH4 is shown to intersect the (110) and (111) planes, which can be seen from Fig. 5d and e. From the initial state, in which the CO2 (CH4) molecule is located within the large pore, the minimum energy path takes the molecule through the 6-MR window and into the next large pore, wherein lies the final state. Unlike the CH4 case, the energetic pathway for CO2 is sensitive to the orientation of the molecule. The barrier value is minimized when the CO2 axis is perpendicular to the plane of the 6-MR window.
To understand the diffusion of CO2 and CH4 from the initial state in one pore, and finally to the final state in the neighboring pore in ZIF-8, a number of steps between the initial and final states are mapped out under the NEB approach as shown in Fig. 5. These pathways are studied in two cases: (1) with the atomic positions of all atoms allowed to relax and (2) with the MOF atomic positions frozen and the CO2 or CH4 atomic positions relaxed. The lattice parameters are kept constant in both cases. In (1), the calculation yields the quasi-vibrational modes of molecules and frameworks together with a combined center of mass oscillation as all ions are relaxed. In (2), the calculations provide the vibrational modes of CO2 against the framework, which also reflects the attempt frequency of CO2 while diffusing through the pores. Fig. 6 shows the minimum energy pathways obtained for CO2 and CH4 for both cases (1) and (2).
The diffusion barrier for CO2 is found to be 3.43 kcal mol−1 in the flexible ZIF-8 (case (1)). The diffusion energy barrier is increased by 25% (∼1 kcal mol−1) in the fixed ZIF-8 (case (2)) as shown in Fig. 6a. Although other gas diffusion paths may exist, this CO2 diffusion energy barrier suggests that CO2 diffusion at 298 K in the flexible ZIF-8 is about 5.0 times larger than in the fixed ZIF-8. This CO2 diffusion difference in the flexible and fixed ZIF-8 structures obtained from the quantum chemistry calculations is similar to the results obtained from the molecular dynamics simulations by using a classical flexible force-field for ZIF-8. Similarly, the diffusion barrier for CH4 diffusion is found to be 3.38 kcal mol−1 in case (1) and is increased by 37% in case (2). Fixing framework atoms has a higher impact on the barrier in the case of CH4. This is expected as the kinetic diameter of the CH4 molecule is larger than that of the window. CH4 experiences a strong repulsion while passing through the window. The peak in the pathway in Fig. 6 for both CO2 and CH4 arises due to gas molecules interacting relatively strongly while they lie near the 6-MR window as shown in Fig. 5(b and e).
The diffusion coefficient can be calculated using the Einstein–Smoluchowski relation:69
![]() | (4) |
![]() | (5) |
Molecule | Methods | D diff (m2 s−1) | Sources |
---|---|---|---|
CO2 | Flexible framework | 8.77 × 10−10 | This work |
Fixed framework | 1.7 × 10−10 | ||
Infrared microscopy | 1.5 × 10−10 | Ref. 42 | |
DFT optimized | 9.31 × 10−11 | Ref. 50 | |
Flexible framework (MD) | 7.93 × 10−10 | Ref. 50 | |
Flexible framework (MC/MD) | 6.2 × 10−9 | Ref. 16 | |
Infrared microscopy | 1.56 × 10−10 | Ref. 70 | |
CH4 | Flexible framework | 3.52 × 10−10 | This work |
Fixed framework | 4.6 × 10−11 | ||
Infrared microscopy | 1.0 × 10−10 | Ref. 42 | |
Flexible framework (MD) | 4.47 × 10−12 | Ref. 50 | |
Flexible framework (MC/MD) | 6.4 × 10−11 | Ref. 16 |
We extend the above calculation to capture the diffusion dynamics from one pore to the next nearest neighbor pore of the same size. The molecule diffuses in through one of the pores in a particular unit cell, migrates across the subsurface, and diffuses out through the neighbor pore of the same unit cell as shown by the initial and final states in Fig. 7 for CO2. The diffusion barrier is symmetrical with two lobes occurring across the two pores.
From the results presented here, we find that loading of individual gas molecules of CO2 until around 2% change in volume does not significantly affect its adsorption energy. At higher loading limit, after loading of about 17 molecules, the volume started changing dramatically with a simultaneous decrease in adsorption energy in the case of CH4 whereas a rapid decrease in adsorption energy is found to occur after loading of 23 molecules in the case of CO2. These results indicate that the maximum loading limit (as in experiment) is more sensitive to adsorption energy for CH4 loading than for CO2 where up to 2% of the framework's volume expansion change in adsorption energy is negligibly small. In the case of a large number of molecules diffusing in the framework, the crowding effect may alter the diffusion coefficient due to interactions among the gas molecules in addition to the molecule–framework interaction. Nevertheless, our results advance the understanding of diffusion mechanisms of gas molecules and can be potentially useful in designing membranes.
By default, DFT does not consider any temperature effect. There exist many closely spaced local energy minima within the ZIF-8 frameworks for the adsorbate molecule. The interaction between them is weak and is mostly of the van der Waal type. It could be challenging to locate global minima in such cases using the DFT method as it does not employ statistics for phase space sampling. Phase space sampling using MD simulation is helpful to locate a global minimum, where statistics is rigorously applied. However, the DFT method could be an alternative method to achieve a good estimate in the system like ZIF-8. This opens up an avenue for the possibility of using the DFT method for diffusion property calculations in other MOF materials for which force field parameters are not readily available and calculation of those parameters is a daunting task. The difference in many closely spaced energy sites is most likely within the DFT energy bar of 0.025 kcal mol−1. In addition, if we compare the computational costs and the significance of expected accuracy in the results, phase space sampling using MD by developing reasonable force field parameters in many instances would not improve the overall quality of the work appreciably. Importantly, ZIF-8 has a highly symmetrical conventional unit cell.
At 0 K, the Gibbs free energy G is equal to the enthalpy H of the system. As the temperature is increased, the entropy of the system changes and the change in entropy contributes to the Gibb's free energy ΔG = ΔH – TΔS. The error in adsorption and diffusion barriers due to a rise of temperature is substantially mitigated due to the cancellation of such error in total energy while talking a difference. In solid systems like fcc Al, the temperature dependence of enthalpy and entropy is shown to nearly cancel out each other under harmonic and quasi-harmonic approximations.71 Nevertheless, our results obtained for volume dependence of adsorption energy and diffusion coefficients provide good estimates and help in understanding of major diffusion channels in ZIF-8 materials.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0re00416b |
This journal is © The Royal Society of Chemistry 2021 |