Tian
Tian
a,
Michael T.
Wharmby
b,
José B.
Parra
c,
Conchi O.
Ania
c and
David
Fairen-Jimenez
*a
aDepartment of Chemical Engineering and Biotechnology, University of Cambridge, Pembroke St., Cambridge CB2 3RA, UK. E-mail: df334@cam.ac.uk
bDiamond Light Source Ltd., Diamond House, Harwell Science & Innovation Campus, Didcot, Oxon., OX11 0DE, UK
cInstituto Nacional del Carbón, INCAR-CSIC, Oviedo, Spain
First published on 26th February 2016
The flexibility and structure transition behaviour of ZIF-8 in a series of samples with different particle size has been studied using a combination of high-resolution N2 gas adsorption isotherms and, for the first time, a broad in situ PXRD and Rietveld analysis. During the stepped adsorption process, large particles showed a narrow adsorption/desorption pressure range with a shorter equilibrium time due to lower kinetic hindrance, deriving from higher amount of active sites. In situ PXRD showed that both the rotation of imidazole ring and a bend in the methyl group led to the gate opening of ZIF-8.
Among all ZIFs, ZIF-8 [Zn(mIM)2] (mIM = 2-methylimidazolate, C4H5N2−) is of particular interest due to its high thermal and chemical stability, and characteristic porosity.8,9 ZIF-8 presents a large BET area (SBETca. 1700 m2 g−1) and relative high crystal densities (0.95 g cm−3), and can be prepared in large and chemically robust monolithic morphologies with large volumetric adsorption capacities (SBET(vol) = 1660 m2 cm−3).10 ZIF-8 possesses a sodalite (SOD) topology containing relatively large pore cavities (ca. 11.6 Å diameter) interconnected by small windows (ca. 3.4 Å diameter).11 Due to its small window size, ZIF-8 was expected to separate molecules with different kinetic diameters. However, it was found that gas molecules with diameter larger than 3.4 Å, e.g. N2 (3.6 Å), could also be adsorbed.12
In a previous study, we combined the use of grand canonical Monte Carlo (GCMC) simulations and in situ powder X-ray diffraction (PXRD) to demonstrate that the structural transition of ZIF-8 was induced by gas adsorption, from an ambient to a high pressure structure, ZIF-8AP and ZIF-8HP, respectively.13 The structural change implies the reorientation of the mIM ligands and the increase in the size of both the 4- and 6-ring windows present in this material. In particular, the swing effect of the imidazole rings is responsible for the stepwise adsorption of N2 at 77 K at 0.02 p/p0. We have also shown the existence of the swing effect and phase transition in ZIF-8 during the adsorption of other alkanes at ca. 0.1 p/p0 and 125 K for methane; as well as 0.3 p/p0, 0.02 p/p0 and 6 × 10−3p/p0 at 273 K for ethane, propane and butane, respectively.14 Interestingly, in all these later cases the phase transition takes place without observing a stepped behaviour in the adsorption process (i.e. Type I adsorption isotherms). Using DFT calculations, we showed that the driving force for this transition was related to the insertion of additional molecules in the 4-ring windows, which in turn stabilizes the “high-loading” ZIF-8HP structure. Following this work, Ania et al. studied the structural transition of ZIF-8 using high resolution adsorption isotherms of different gases (CO, N2, O2, Ar) at different temperatures (77 and 90 K) and found that polarizability, size and shape of gas molecules affected the swing effect and phase transition behaviour.15 Contrary to our previous results,13 Ania et al.15 and Park et al.,1 observed hysteresis during the desorption of N2 and Ar at cryogenic temperatures as well as some other gases such as CO and O2, and reported the existence of two adsorption substeps around the transition pressure, attributing this behaviour to the reorganization of the adsorbed gas molecules and their interactions with the framework. In addition to these studies, the swing effect of ZIF-8 has been further studied both experimentally and computationally.16,17 For example, by using single crystal XRD and Raman spectroscopy, it has been observed that the structural transition of ZIF-8 was also induced by reducing the temperature beyond 150 K under N2.18,19 We and others have also shown recently the use of THz radiation and inelastic neutron scattering for the study of the flexibility modes of ZIF-8.20–22 Although the adsorption mechanism of multiple gases as well as the shape of the adsorption isotherms of ZIF-8 have been widely investigated, the origin of the substeps in the adsorption isotherms of various gases at cryogenic temperatures, and of the occurrence of a hysteresis loop during the desorption are still unclear.
During the preparation of this manuscript, Zhang et al. demonstrated that the particle size of ZIF-8 affected its structural transition during N2 adsorption at 77 K.23 By applying the osmotic framework adsorbed solution theory (OFAST),24 they observed that the pressure for the phase transition shifted to higher values when the particle size was reduced. This finding is especially relevant because ZIF-8 can be synthesized with particle sizes in the nanometre to micrometre range.25 In this regard, the flexibility of ZIF-8 is of particular importance as it affects the diffusivity of gas molecules through the porous network. Indeed, previous molecular dynamic studies have shown tremendous differences in the diffusivity of gas molecules between the rigid and the flexible structure of ZIF-8.14,26,27
In this work, we synthesised ZIF-8 with a wide range of particle sizes and studied the role of the particle size on the gas-induced structural transition on adsorption of N2 at 77 K, as well as the kinetics of the process. We also include, for the first time, a detailed study about the mechanisms of the phase transition during the adsorption process by using in situ PXRD experiments during N2 adsorption at cryogenic temperature on two different sized ZIF-8 samples.
ZIF-8-0.21, ZIF-8-0.61 and ZIF-8-1.6 were synthesised based on methods reported by Kida et al.29 Aqueous solutions of Zn(NO3)·6H2O (10 ml, 0.25 M) and 2-methylimidazole (90 ml, 2.77 M) were mixed and stirred for 24 h under ambient conditions to synthesise ZIF-8-0.21. ZIF-8-0.61 and ZIF-8-1.6 were synthesised by reducing the concentration of 2-methylimidazole to 1.66 M and 1.11 M, respectively, while keeping other conditions unchanged.
ZIF-8-10 and ZIF-8-98 were synthesised based on methods reported by Zhang et al.30 For ZIF-8-10, 2-methylimidazole (0.099 M) and sodium formate (0.198 M) were dissolved in 40 ml of methanol, which was then poured into a solution of Zn(NO3)·6H2O (40 ml, 0.049 M) in methanol. The mixture was placed in a Teflon-lined steel autoclave and heated at 363 K overnight. For ZIF-8-98, 2-methylimidazole (0.593 M) and sodium formate (0.281 M) were dissolved in 40 ml of methanol, and then poured into a solution of Zn(NO3)·6H2O (40 ml, 0.3 M) in methanol. The mixture was placed in a Teflon-lined steel autoclave and heated at 363 K overnight.
All as synthesised solids were collected by centrifugation, washed with ethanol (20 ml, 3 times) and dried at 373 K under vacuum.
High resolution N2 adsorption–desorption isotherms were measured at 77 K using a Micromeritics ASAP 2020 instrument in the relative pressure range from 10−6 to 1 for the adsorption branch and down to 10−3 for the desorption branch. The instrument was equipped with a molecular drag vacuum pump and three pressure transducers (0.1, 10, 1000 mmHg, uncertainty within 0.15% of reading) to enhance the sensitivity in the low-pressure range. All samples were evacuated overnight for 24 h at 423 K under dynamic vacuum prior to adsorption. Strict analysis conditions were programmed during the gas adsorption measurements to ensure equilibrium data in all cases. Consequently, the average elapsed time for the measurement of the isotherms was 90–120 h, with over 200 equilibrium points on average per isotherms. The saturation pressure of the gas was continuously measured throughout the analysis by means of a pressure transducer. For all isotherms, warm and cold freespace correction measurements were performed by using ultrahigh purity He gas (grade 5.0, 99.999% purity). Ultrahigh purity N2 (i.e. 99.9992%) was provided by Air Products.
Data for the in situ gas adsorption powder X-ray diffraction experiments were collected at beamline I11 at Diamond Light Source (Oxon., UK; λ = 0.825701 Å), using the I11 gas cell at 80 K.31–33 Both ZIF-8-0.14 and ZIF-8-98 samples were activated prior to the gas adsorption experiment by heating to 413 K under vacuum. Data were collected at 22 different pressures, 11 in the 0.0–0.1 p/p0 range and 11 in the 0.1–0.9 p/p0 range. Selected pressure points were fully refined by the Rietveld method using the TOPAS-Academic v5 suite.34 Subsequently, these structures were used as fixed points in a parametric Rietveld refinement,35 also performed using TOPAS-Academic v5. See ESI† for full details.
Materials | Particle size (μm) | BET areas (m2 g−1) |
---|---|---|
ZIF-8-0.14 | 0.141 ± 0.035 | 1740 |
ZIF-8-0.21 | 0.213 ± 0.041 | 1702 |
ZIF-8-0.61 | 0.608 ± 0.013 | 1739 |
ZIF-8-1.6 | 1.6 ± 0.4 | 1799 |
ZIF-8-10 | 10.0 ± 0.5 | 1809 |
ZIF-8-98 | 98 ± 34 | 1731 |
We further analysed the porosity for all samples by using high-resolution N2 adsorption/desorption isotherms at 77 K. The BET areas of ZIF-8 samples were calculated using the Rouquerol's consistency criteria.37Table 1 shows similar BET areas, ranging between 1700 and 1810 m2 g−1, for all the materials. Fig. 3 shows the adsorption isotherms of N2 on ZIF-8 samples; we used a semi-log plot to give more detail in the low pressure range. All samples presented the same characteristic stepwise adsorption isotherm widely reported for ZIF-81,13,15 and the swing effect of the mIM rings at a similar onset transition, of ca. 5 × 10−3p/p0. As particle size decreases, the samples showed a slight increase in the uptake during the plateau and close to saturation. This effect is attributed to the existence of interstitial spaces between ZIF-8 particles, something that is more important for smaller particle size.
Fig. 3 (inset) shows a magnification of the adsorption isotherms in the phase transition region for clarity, whereas Fig. 4 and S4† show the adsorption and desorption branches, where the differences between samples of different size are evident. On the one hand, samples with smaller particle size (i.e. ZIF-8-0.14 and ZIF-8-0.21) shifted the gate pressure for the adsorption substep to higher p/p0 values compared to samples with larger particle size (i.e. from 0.006 to 0.014 p/p0). On the other hand, samples with smaller particle size showed a wider hysteresis loop compared to larger particle samples. Interestingly, in the case of ZIF-8-98 (i.e. the sample with the largest particle size) the adsorption was completely reversible and no hysteresis was found. These phenomena agreed well with the previous reports from Ania et al., which indicate the presence of hysteresis15 – and from some of us which show no such hysteresis.13 As the hysteresis loop occurred at very low pressures, it cannot be explained by a capillary condensation mechanism characteristic of mesoporous materials.38 This is attributed to the desorption of N2 from an opened ZIF-8HP structure and hence gradual rearrangement of the remaining gas molecules during the transition to the initial closed ZIF-8AP structure.
Fig. 4 Low pressure hysteresis loops of ZIF-8 samples with different particle size. Closed symbols, adsorption; open symbols, desorption. |
The shift of the gate pressure to higher values and the wider range of the hysteresis loop for smaller particles suggests a higher energy barrier between ZIF-8AP and ZIF-8HP structures, and thus, a higher kinetic hindering. A similar observation was made by Watanabe et al.,39 using GCMC simulation, where they showed that the width of the hysteresis loop was inversely proportional to the energy barrier. This transition was also observed experimentally by Sakata et al.40 for the interpenetrated [Cu2(bdc)2(bpy)]n (bdc = 1,4-benzenedicarboxylate, bpy = 4,4′-bipyridine) system, which exhibited a cooperative guest-induced structural transformation from a non-porous closed phase to an open phase.
Fig. 5 Adsorption isotherms and equilibrium time for N2 at 77 K for (a) ZIF-8-98 and (b) ZIF-8-0.14. Black circles, N2 uptake; red empty squares, equilibrium time; grey shading, gate pressure region. |
There are three clear ranges taking into account the phase transition occurring at a gate pressure of ca. 5 × 10−3p/p0. First, before the gate pressure, all the ZIF-8 samples showed very different adsorption kinetics, with the equilibration process being much faster for samples with larger particle sizes (e.g. ZIF-8-98, ca. 25 min) than for those with smaller particle size (e.g. ZIF-8-0.14, ranging from 100 to 25 min). Second, the adsorption rate for all samples was significantly slowed near the gate pressure (up to ca. 80 min), before increasing at p/p0 > 0.01 (i.e. after the transition to ZIF-8HP), which is the third stage with all the samples behaving similarly, with faster kinetics than at lower pressures. The dependency of the adsorption equilibration time and particle size before the gate pressure implies differences in ZIF-8AP flexibility (i.e. the dynamic, free swing of the mIM rings to allow access of large molecules into the cavities) between different samples at these conditions. In turn, the increase of the equilibration time (i.e. slow kinetics) at the gate pressure would be explained by the permanent swing of the mIM rings, and the access and rearrangement of the gas molecules in the cavities, and therefore the phase transition of ZIF-8AP to ZIF-8HP observed within this pressure range.13 Above the transition pressure all samples present somewhat similar kinetics due to the opened structure. These observations are consistent with a model in which, below the gate pressure, the 3.4 Å windows in ZIF-8 change between open and closed configurations depending on environmental pressure gradient; whereas above the gate pressure, the windows remain open. In these conditions of static low and high pressure, a snapshot of the ZIF-8 structure would reveal (at least) two discrete, non-disordered phases ZIF-8AP and ZIF-8HP rather than a free swing of the mIM rings. This dynamic opening and closing process during the pressure swing depends on the particle size of ZIF-8.
Assuming that both the large and the small particles have the same structure, we are effecting the same change in both samples. From a thermodynamic point of view, the energy change should also be the same for both. As shown above, we found however important differences between samples, so the question would be: are there any differences in the composition or the structure of the different ZIF-8 crystals? The main difference is related to the ratio of external surface vs. bulk phase, which is much larger for smaller than for larger particles. This would imply that the external surface/bulk ratio of the 4-ring windows in large particles sizes is lower than in small particles. If the driving force for the phase transition is the adsorption of additional N2 molecules in the 4-ring window in the bulk14 – similar to the hand-glove model of enzymes – we will have a higher amount of active sites in larger particles, therefore reducing the activation energy for the phase transition.
Fig. S8† shows the first adsorption site (observed from 0.0371 bar and 0.816 molec per uc), associated to Site I as described previously for adsorption of methane in ZIF-8.14 At ca. 42 molec per uc (i.e. the first step at 0.0575–0.0685 bar), two additional N2 adsorption sites in the centre of the 6-ring window (4.4 molec per uc) and the centre of the cavity (2.1 molec per uc) were identified, associated to Site II and Site III, respectively.42 We termed this phase as the gate-closed large cell ZIF-8AP structure.
Interestingly, the PXRD patterns of ZIF-8-98 obtained at 0.0798 and 0.0868 bar showed a splitting of the diffraction peaks into three, indicating a mixture of different phases – probably due to the lack of equilibrium (Fig. S10†). These three phases were identified as: (i) the gate-closed ZIF-8AP structure (at p(N2) = 0 bar), (ii) the gate-closed large cell ZIF-8AP structure (observed at p(N2) = 0.0685 bar) and the (iii) gate-open ZIF-8HP structure (p(N2) > 0.0948 bar). Increasing adsorption pressure led to further N2 uptake and the completion of the phase transition to ZIF-8HP structure. This was accompanied by an increase in unit cell volume by ca. 0.9% (∼50 Å3); after the phase transition, the unit cell volume remains approximately constant. At ca. 52 molec per uc (i.e. the second adsorption substep at pressures above 0.100 bar), the rotation of the mIM rings increased to 16°, whereas the methyl group bent between 9 and 12°. A new adsorption site was identified in the 4-ring window with 3.9 molec per uc adsorbed (Fig. S8c†), and was associated to Site IV.14
Increasing the adsorption pressure to 0.3500 bar led only to a small increase in the number of adsorbed molecules (i.e. 35 molec per uc). One additional adsorption site was observed at this pressure: Site IV 0.5 molec per uc at the centre of the 4-ring window. In this range of pressure, both the rotation of the mIM rings and the bend of the methyl groups increased gradually up to 7.6 and 3.9°, respectively. This was similar to the values obtained for the large crystal ZIF-8-98 and the intermediate step during the adsorption process. Increasing the adsorbate pressure led to a jump in the N2 uptake, reaching 61 molec per uc at 1.3580 bar, as well as an increase in the rotation of the mIM rings and the bend of the methyl group up to 10 and 9°, respectively. This jump corresponds to the transition to the ZIF-8HP.
Overall, ZIF-8-0.14 showed a broader pressure range of transition compared with ZIF-8-98 (Fig. S6†). This is similar to the trend of phase transition in N2 isotherms at 77 K. The overall rotation of the imidazole for ZIF-98 and ZIF-8-0.14 is 17° and 10°, respectively as shown in Fig. S11 and S12.†
Footnote |
† Electronic supplementary information (ESI) available: Particle size distribution, gas adsorption, BET areas calculation, and in situ X-ray diffraction. See DOI: 10.1039/c6dt00565a |
This journal is © The Royal Society of Chemistry 2016 |