Open Access Article
Zhehao
Huang
*,
Tom
Willhammar
* and
Xiaodong
Zou
*
Department of Materials and Environmental Chemistry, Stockholm University, Stockholm, SE-106 91, Sweden. E-mail: zhehao.huang@mmk.su.se; tom.willhammar@mmk.su.se; xzou@mmk.su.se
First published on 15th December 2020
Porous crystalline materials such as zeolites, metal–organic frameworks (MOFs) and covalent organic frameworks (COFs) have attracted great interest due to their well-defined pore structures in molecular dimensions. Knowing the atomic structures of porous materials is crucial for understanding their properties and exploring their applications. Many porous materials are synthesized as polycrystalline powders, which are too small for structure determination by X-ray diffraction. Three-dimensional electron diffraction (3DED) has been developed for studying such materials. In this Minireview, we summarize the recent developments of 3DED methods and demonstrate how 3DED revolutionized structural analysis of zeolites, MOFs, and COFs. Zeolites and MOFs whose structures remained unknown for decades could be solved. New approaches for design and targeted synthesis of novel zeolites could be developed. Moreover, we discuss the advances of structural analysis by 3DED in revealing the unique structural features and properties, such as heteroatom distributions, mixed-metal frameworks, structural flexibility, guest–host interactions, and structure transformation.
Single crystal X-ray diffraction (SCXRD) is currently the most predominant method for crystal structure determination at the atomic level, but requires large crystals (approx. 5 × 5 × 5 μm3 for advanced synchrotron sources) with sufficient quality. Many zeolites, MOFs, and COFs are obtained as nano- or submicron-sized crystals. It is often difficult and sometimes even impossible to grow large crystals, for example, 1D and 2D crystals based on reaction kinetics and thermodynamics. Meanwhile, the sample quantities can be very limited, which prevents testing of enough synthesis conditions for growing large crystals. Although powder X-ray diffraction (PXRD) can be used for structure determination of polycrystalline materials, peak overlap often limits unambiguous indexing and correct intensity estimation of reflections. Structure determination from PXRD becomes more problematic when a crystal has large unit cell dimensions and/or the sample contains multiple phases. During the past two decades, three-dimensional electron diffraction (3DED) techniques,10,11 also known as microcrystal electron diffraction (MicroED),12,13 have been developed to overcome the barriers for structural analysis of nano- and submicron-sized crystals, which drastically accelerated the development in the fields of zeolites, MOFs, and COFs.
In this Minireview, we summarize recent advances in the development of 3DED methods, and their applications for structural analysis of zeolites, MOFs, and COFs, with emphasis on novel materials (Scheme 1). We describe how the detailed atomic structures obtained by 3DED can reveal the key structural features, uncover structure–property relationships, and provide an increased understanding of the properties of such materials. We believe that a timely critical overview on 3DED is of great importance to advance our understanding towards the development of porous crystalline materials, where large crystals are no longer mandatory for single crystal structural analysis.
![]() | ||
| Scheme 1 Novel structures and topologies of zeolites, MOFs, and COFs can be determined by using 3DED techniques from nano- and sub-microcrystals, revealing their underlying properties. | ||
The revolution in electron crystallography started in late 2000s, when Kolb et al. and Hovmöller et al. independently developed software to control three-dimensional electron diffraction (3DED) data collection (Fig. 1). This is achieved by stepwise rotation of a crystal along an arbitrary axis and collecting ED patterns at each angle. They applied different approaches to achieve fine sampling of reciprocal lattices and integrated intensities. Kolb et al. developed automated diffraction tomography (ADT)10 to control the goniometer tilt and applied beam precession using dedicated hardware.23 Hovmöller et al. developed rotation electron diffraction (RED) software to control both the goniometer and the electron beam, combining coarse crystal rotation with fine beam tilt11,24 (Table 1). They have also developed software to process the ED patterns from different angles and reconstruct 3D reciprocal space (Movie S1†). This new development of 3DED no longer requires alignment of crystals at a zone-axis, which both simplifies and speeds up the data collection. More importantly, dynamical effects are reduced when crystals are oriented off zone axes. As a result, intensities obtained from 3DED are less dynamical and in most cases can be treated as kinematical intensities, so that similar routines and software for structure determination by SCXRD can be applied on 3DED data (Fig. 1).25
![]() | ||
| Fig. 1 Comparison of PXRD and 3DED data obtained from polycrystalline samples. PXRD patterns are collected from millions of crystallites but contain a limited number of peaks due to overlapping. Meanwhile, 3DED data are collected from different orientations of a single crystal, and provide thousands of diffraction peaks.31 3DED data contain a missing wedge due to the rotation angle limited by the TEM while PXRD patterns are complete. Reproduced from ref. 31 with permission from the Elsevier, copyright 2014. | ||
| Protocol name | Data collection strategy | Data collection speed | EM mode | Beam tilt | Beam procession | Publication year | Reference |
|---|---|---|---|---|---|---|---|
| ADT/PEDT | Stepwise | Slow | STEM | No | Yes | 2007 | 10 and 26 |
| RED | Stepwise | Slow | TEM | Yes | No | 2010 | 11 and 24 |
| EDT | Stepwise | Slow | TEM | Yes | No | 2013 | 27 |
| MicroED | Stepwise | Slow | TEM | No | No | 2013 | 12 |
| Rotation electron diffraction | Continuous | Fast | TEM | No | No | 2013 | 28 |
| MicroED | Continuous | Fast | TEM | No | No | 2014 | 13 |
| Fast EDT | Continuous | Fast | TEM | No | No | 2015 | 29 |
| cRED | Continuous | Fast | TEM | No | No | 2018 | 30 |
Porous materials often suffer from electron beam damage, which diminishes the crystallinity during data acquisition. This drawback has been overcome with the development of continuous rotation data collection methods,13,28–30 another milestone in the development of 3DED. In contrast to stepwise methods, the crystal rotates continuously during the entire data collection, and ED patterns are simultaneously recorded as a movie. Because the crystals can move during rotation, Cichocka et al. developed a crystal-tracking solution to ensure that the same crystal area is used during the data acquisition.30 By adjusting the goniometer rotation speed and using a fast detector, a full 3DED dataset can be acquired very fast, within one minute (Movie S2†). The short data collection time combined with a low electron dose rate (<0.1 e s−1 Å−2) can effectively reduce the total electron dose applied on the analytes, and thus minimizes beam damage. Notably, because intensities are integrated on each ED pattern by using the continuous rotation methods, more accurate intensities can be obtained compared to step-wise rotation. Several groups developed their own protocols for continuous rotation electron diffraction data collection, including microcrystal electron diffraction (MicroED),13 fast electron diffraction tomography (Fast EDT29), rotation electron diffraction28 and continuous rotation electron diffraction (cRED)30 (Table 1). By applying 3DED, thousands of diffraction peaks can be obtained from a single nanocrystal, while only a limited number of peaks are observed in a PXRD pattern (Fig. 1). The high number of datasets ensures an unambiguous determination of unit cells and space groups. Moreover, it allows reliable structure solution and refinement despite the presence of dynamical effects.
To process 3DED data, several programs, such as RED data processing,24 ADT 3D,32 and PETS33 have been developed. Software developed and widely used for X-ray crystallography, such as XDS34 and DIALS,35 can be easily adapted for processing continuous rotation electron diffraction data. After extracting the intensities, structure solution and refinement are performed by applying similar approaches to those for SCXRD, using direct methods (SHELX,36 SIR37) or charge flipping (JANA38). These programs treat intensities from 3DED as kinematical. Atomic scattering factors for electrons are used instead of those for X-rays. Various structures of porous crystalline materials, which could not be solved by X-ray diffraction, have now been determined by 3DED, and the number keeps growing.
The strong interactions between electrons and matter allow obtaining high-resolution 3DED data from nano- and submicron-sized crystals. On the other hand, the strong interactions also lead to multiple dynamical scattering events that cause the intensities of electron diffraction to deviate from kinematical intensities. Although dynamical effects can be reduced by 3DED techniques where data are collected at arbitrary off-zone orientations, structure refinements against 3DED data still lead to relatively high R1 values, which indicate the difference between the intensities calclated from the model and the observed intensities. This has been a key concern in the field of electron crystallography. Nevertheless, recently studies on zeolites39 and MOFs40 have proved that the structural models obtained from 3DED are similar to those obtained by SCXRD, with average deviations of atomic positions by less than 0.10 Å. Moreover, structure determination by 3DED also shows high reproducibility. The deviations of atomic positions obtained from 3DED data of different individual crystals are less than 0.10 Å.39,41 The high accuracy is attributed to the high ratio (usually >10) of the number of reflections to the number of refined parameters. Despite the high accuracy, however, because the structure refinement software used is developed for SCXRD, which assumes the intensities to be kinematical. Thus, the resulting high R1 values mainly show that the observed 3DED intensities are dynamical and deviate from kinematical intensities calculated from the model. To take dynamical effects into account and further improve the structure refinement, Palatinus et al. developed dynamical refinement42 that compares observed PEDT intensities with dynamical intensities. They show that more accurate structural models including hydrogen positions could be obtained, and R1 values were greatly reduced.
An important advantage of 3DED is that each dataset can be collected in a few minutes. Thus, it is possible to collect datasets from many different individual crystals. The different datasets can be merged to overcome the missing-wedge problem and achieve high data completeness. Based on this, a new approach, serial rotation electron diffraction (SerialRED), has recently been developed to achieve fully automated crystal screening and cRED data collection.43 It is capable of screening up to 500 crystals per hour. A data processing pipeline using hierarchical cluster analysis (HCA) was developed to automatically process each dataset, group the phases present in the sample and find the best matching datasets for subsequent data merging and structural analysis (Fig. 2). This method enables high-throughput phase analysis and structure determination, which is especially powerful for studies of porous crystalline materials.
![]() | ||
| Fig. 2 Dendrogram of the hierarchical cluster analysis (HCA) based on unit cells. Two clusters are identified, corresponding to mordenite and ZSM-5, shown as yellow and black clusters respectively (marked by circles).43 Reproduced from ref. 43 with permission from the International Union of Crystallography, copyright 2019. | ||
| Name | Code | Channel dimension and sizes | 3DED method | Publication year | Name | Code | Channel dimension and sizes | 3DED method | Publication year |
|---|---|---|---|---|---|---|---|---|---|
| a Disordered structures. | |||||||||
| ITQ-43 | 3D, 28 × 12 × 12 | ADT | 2011 | SCM-14 | SOR | 3D, 12 × 10 × 10 | RED | 2017 | |
| ITQ-51 | IFO | 1D, 16 | RED | 2013 | ZSM-43 | MRT | 2D, 8 × 8 | RED | 2017 |
| EMM-23a | -EWTa | 3D, 24 × 10 × 10 | RED | 2014 | IM-18a | UOEa | 1D, 10 | RED | 2018 |
| ERS-18 | EEI | 2D, 8 × 8 | RED | 2014 | PST-13 | POR | 3D 8× | cRED | 2018 |
| IM-17 | UOV | 3D, 12 × 10 × 8 | ADT | 2014 | PST-14 | ||||
| PKU-16 | POS | 3D, 12 × 11 × 11 | RED | 2014 | SYSU-3 | -SYT | 3D, 24 × 8 × 8 | cRED | 2018 |
| PST-6 | PSI | 1D, 10 | RED | 2014 | ECNU-21 | EWO | 1D, 10 | cRED | 2019 |
| CIT-7 | CSV | 2D, 10 × 8 | RED | 2015 | EMM-37 | ETV | 3D, 10 × 8 × 8 | cRED | 2019 |
| ITQ-53 | -IFT | 3D, 14 × 14 × 14 | RED | 2015 | SCM-15 | 3D, 12 × 12 × 10 | cRED | 2019 | |
| ITQ-54 | -IFU | 3D, 20 × 14 × 12 | RED | 2015 | SSZ-27 | 1D, 8 | cRED | 2019 | |
| PST-20 | 3D, 8 × 8 × 8 | RED | 2015 | ECNU-23 | 2D, 12 × 8 | cRED | 2020 | ||
| ZSM-25 | MWF | 3D, 8 × 8 × 8 | RED | 2015 | IDM-1a | 3D, 17 × 17 × 10 | cRED | 2020 | |
| CIT-13a | CTHa | 2D, 14 × 10 | RED | 2016 | PST-5 | 2D, 10 × 8 | cRED | 2020 | |
| EMM-26 | EWS | 2D, 10 × 10 | RED | 2016 | PST-24a | 2D, 10 × 10 | cRED | 2020 | |
| 3D, 10 × 10 × 10 | |||||||||
| ITQ-58 | 2D, 8 × 8 | Fast-EDT | 2016 | RUB-5a | 2D, 8 × 8 | ADT | 2020 | ||
51 and IDM-1
52) were determined by 3DED. ITQ-43 with extra-large 28-ring channels was the first example of ab initio structure determination of zeolites by 3DED.45 The structure was solved by direct methods using SIR2008, from ADT data collected at 100 K. All 20 Ge/Si and 13 out of 24 oxygen positions were directly located. The missing oxygen atoms were either found from difference Fourier maps or placed according to the geometry. It is common that extra-large pore zeolites have interrupted frameworks with terminal T-OH groups. The terminal groups can be further disordered, such as in EMM-23,51 and IDM-1.52 The average structures of EMM-23 and IDM-1 were also determined by 3DED. In these cases, solid-state NMR played a complementary role in determining the disordered terminals in the framework. Disorder can be utilized to alter the pore sizes. For instance, the channel size of EMM-23 can be varied from 21- to 24-ring channels by tuning the occupancy of one terminal T-site from one to zero.
Using a hybrid pixel detector with fast data readout, fast cRED data collection (0.5–3 min per dataset) can be achieved. Combined with the crystal tracking system, it enables collecting multiple datasets over a large angular range to reach a high data completeness. The aluminosilicate EMM-37 also crystallizes in a triclinic system.54 Using cRED, the merged dataset reached a completeness of 79.8% with a resolution of 0.70 Å, allowing refinement of anisotropic ADP of each individual atom without applying any restraints. cRED revealed that EMM-37 has a 3D 10- and 8-ring channel system. Notably, the pore widows became narrowed after calcination. This property is crucial for EMM-37 to selectively separate ethane over propane.
![]() | ||
| Fig. 3 The structural model of EMM-26. (a) The cavity in EMM-26 where the organic template resides. (b) The structural model viewed along the [100] direction. The Si-rich sites are shown in blue and B-rich sites in pink and magenta. Bridging O atoms have been omitted for clarity.55 Reproduced from ref. 55 with permission from the Royal Chemical Society, copyright 2016. | ||
![]() | ||
| Fig. 4 Structural transformation from germanosilicate CIT-13 into high silica zeolites ECNU-21 and ECNU-23.57 The structures of all three zeolites were discovered by cRED. Reproduced from ref. 57 with permission from the John Wiley & Sons, Inc., copyright 2020. | ||
3DED has also been used for ab initio structure determination of (silico)aluminophosphates such as ITQ-51,50 PST-5,60 PST-6,61 PST-13
62 and PST-14.62 PST-5 and PST-13 were synthesized using diethylamine (DEA) as the organic structure directing agent (OSDA). Both structures contain penta-coordinated framework Al atoms connected by OH groups, which could be removed by calcination to form fully tetra-coordinated frameworks PST-6 and PST-14, respectively (Fig. 5). Interestingly, while the framework topology of PST-14 remained the same as that of PST-13, the framework topology of PST-5 changed completely upon dehydration. It is revealed by cRED that the structural changes between PST-5 and PST-6 involve a 3D–3D topotactic transformation between the dcc/d4r and nsc. Based on the structural transformations observed in PST-5/PST-6, AlPO-21/AlPO-25 and AlPO-C/AlPO-D pairs, a new approach to generate novel structures in aluminophosphate systems was proposed. A series of new synthetically feasible zeolite structures were predicted.60 This approach could enable targeted synthesis of novel phosphate-based materials, which adds to the existing ADOR approach applicable on silicate-based zeolites.
![]() | ||
| Fig. 5 Structural transformation of the aluminophosphates PST-13 to PST-14 (a) and PST-5 to PST-6 (b) by dehydration. While the transformation from PST-13 to PST-14 only involves the loss of bridging –OH groups and no change in the framework topology,62 that from PST-5 to PST-6 involves transformation of double-crankshaft chain (dcc) and double four ring (d4r) units to narsarsukite-type chain (nsc) units. This leads to topologically different framework structures. The pore system of the material changes from the 2D channel system of PST-5 to a 1D system.60 Reproduced from ref. 60 and 62 with permission from the John Wiley & Sons, Inc., and Springer Nature, respectively, copyright 2018 and 2020. | ||
63 and ZSM-25,64 and a 2D natural zeolite cowlesite65 were solved from 3DED data after remaining elusive for more than 30 years. The structures of ZSM-43 and cowlesite were solved from high-resolution (∼1.0 Å) data by direct methods. However, the RED data from ZSM-25 had a low resolution (2.5 Å) due to beam damage. The data were sufficient for the determination of the unit cell and space group, but not enough for ab initio structure solution. Fortunately, two structures in the IZA Database of Zeolite Structures,44 namely rho (RHO) and paulingite (PAU) were found to be closely related to ZSM-25. All three have the same space group and their unit cell parameters differ by 20 and 10 Å, respectively. More importantly they show a similar intensity distribution of reflections in reciprocal space, indicating that they belong to the same RHO family. The structure of ZSM-25 was thus solved using the strong reflections, with their intensities from RED, and the corresponding crystallographic structure factor phases from PAU. The structural details of ZSM-25 revealed by 3DED opened up an isoreticular route for targeted synthesis of new zeolites. By adopting this new strategy, five new embedded isoreticualr zeoliltes in the RHO family were synthesized (Fig. 6).64,66,67
![]() | ||
| Fig. 6 Framework representations of cross sections (ca. 12 Å thick) of RHO-G1 to RHO-G8 in the RHO family of embedded isoreticular zeolites.66 Reproduced from ref. 66 with permission from the John Wiley & Sons, Inc., copyright 2016. | ||
68 and SSZ-61.69 PST-24 is one example, whose average structure was determined ab initio by cRED.70 The reconstructed 3D reciprocal lattice of PST-24 is featured by alternating rows of sharp spots and diffuse lines (Fig. 7a), indicating the presence of disorder. Nevertheless, from the cRED data, the average unit cell could be determined to be a = 24.14(4) Å, b = 5.21(2) Å, c = 21.76(3) Å, α = 90°, β = 111.4(2)°, γ = 90°. The possible space groups were deduced as C2, Cm, and C2/m. Based on this structural information, intensities were extracted from the sharp Bragg reflections, from which an average structure was obtained directly. The revealed structure is composed of cas-zigzag chains and d5r units. However, the d5r units in the average structure are so close to each other that only every second d5r units in the column along b-axis can exist simultaneously. This resulted in disorders in PST-24 with either open or closed 10-ring channels (Fig. 7b and c), and a channel system varying locally from 2D to 3D (Fig. 7e). HRTEM imaging further shows that PST-24 crystals are 1D ordered and 2D disordered. A similar 3DED approach was used to determine the average structure of the disordered silicogermanate IM-18 (*UOE), with HRTEM imaging revealing the locations and domain sizes of the multidimensional stacking disorders.71 Furthermore, 3DED was also used to quantify stacking disorder in zeolite Beta.72 It is important to note that compared to HRTEM imaging, 3DED data acquisition is simpler and requires much lower electron doses. In addition, 3DED can be performed on a standard TEM and there is no need to align the crystals. On the other hand, HRTEM imaging gives additional information about local structural details, and is complementary to 3DED for studying disordered materials.
![]() | ||
| Fig. 7 (a) Reconstructed reciprocal lattice from a crystal of PST-24. Viewing along (b) [001] and (c) [101] directions, different arrangements of the d5r units (in blue) result in open or closed 10-ring channels. (d) The pockets, which are formed instead of the additional channels, are observable, and one of them is marked in red. (e) One hypothetical model of the channel system of PST-24. The 10-ring channels along the [101] direction are randomly distributed (shown in light blue).70 Reproduced from ref. 70 with permission from the John Wiley & Sons, Inc., copyright 2020. | ||
![]() | ||
| Fig. 8 (a) Structural model of bismuth subgallate, showing a one-dimensional structure with the framework stabilized by intermolecular H-bonding. The 3DED data collected by (b) stepwise RED and (c) cRED with sample cooling, showing a significant improvement on data resolution.75 Reproduced from ref. 75 with permission from the Royal Chemical Society, copyright 2017. | ||
Notably, 3DED is particularly powerful for determining unexpected novel structures. For example, cRED was used to determine the structure of PCN-226, which was crystallized in nano-size.84 From cRED data with 1.15 Å resolution, PCN-226 was revealed to contain rare Zr-oxo chains connected by tetrakis(4-carboxyphenyl)porphyrin (TCPP) ligands (Fig. 9a). The unique chain-based structure revealed by cRED showed an unexpected new topology, which was utilized as a favorable system for electrochemical applications. COK-47 is another novel nanocrystalline MOF determined from ab initio using cRED data.82 Its structure is built from Ti-oxo sheets linked together by 4,4′-biphenyldicarboxylate (bpdc2−) anions (Fig. 9b).
![]() | ||
| Fig. 9 Structural models and corresponding cRED data of (a) PCN-226, a Zr-oxo chain based structure,84 and (b) COK-47, a Ti-oxo sheet based structure.82 All of them are obtained as nanocrystals, and the structures are determined using cRED data. Reproduced from ref. 82 and 84 with permission from the American Chemical Society, and John Wiley & Sons, Inc., respectively, copyright 2020 and 2019. | ||
As aforementioned in zeolite materials, 3DED can be used to determine preferred locations of different metals based on intensities of electrostatic potentials and coordination environments. This is particularly important for mixed-metal MOFs, where the positions of different metal cations play important roles in their physical and chemical properties. For example, the Ti(IV) and Zr(IV) cations as well as the linkers in PCN-415 could be located directly by ab initio structure determination using cRED data, which show that Ti(IV) and Zr(IV) cations form bimetallic [Ti8Zr2O12(COO)16] clusters.78 The structural details of PCN-415 provided important insights on the mechanism of its photocatalytic properties, and could explain how electrons transfer from the valence bands at the NH2- functionalized benzene-1,4-dicarboxylic acid (BCD) to the conduction bands of the Ti(IV) cations.
Cooperative structural changes are one of the unique properties of MOFs. The structural transformability can be triggered by the host–guest interactions as well as by external stimuli.86,87 However, the structural changes may lead to loss of crystallinity that prevents the structural analysis by SCXRD. 3DED was applied to gain an insight into the structural changes of MOFs at atomic levels, and to understand their properties and functionalities. The first example of novel flexible MOFs studied by 3DED was PCN-128,74 which is built from 8-connected Zr6 clusters and a chromatic linker 4′,4′′′,4′′′′′,4′′′′′′′-(ethene-1,1,2,2-tetrayl)tetrakis-(([1,1′-biphenyl]-4-carboxylic acid)) (H4ETTC). PCN-128 exhibits a piezofluorochromic behavior, with reversible color changes from white (PCN-128W) to yellow (PCN-128Y) upon the change of the exterior or interior pressures. RED revealed that the color change was associated with the change of the bond angles within the rectangular ETTC ligand, making PCN-128W a microscissor lift. The flexible structure led to a bathochromic shift in luminescence, which is a promising property for applications in photocatalysis and sensing (Fig. 10a). Another flexible MOF SU-100, with a composition [Bi(BPT)]·2MeOH (BPT = biphenyl-3,4,5-tricarboxylate), was solved by cRED.81 SU-100 exhibits a solvent-dependent flexibility and undergoes reversible volume changes. Different from most flexible MOFs, such as PCN-128 where the flexibility occurs through conformational changes of organic linkers, the breathing phenomenon in SU-100 was associated with its inorganic building unit, a Bi2O12 dimer, through the changes of the O–Bi–O angles (Fig. 10b).
![]() | ||
Fig. 10 Illustration of structural changes in flexible MOFs (a) PCN-128 74 and (b) SU-100,81 whose structures are revealed by 3DED. Reproduced from ref. 74 and 81 with permission from the American Chemical Society, copyright 2015 and 2019. | ||
Beyond the framework structures, 3DED has also shown to be effective in studying host–guest interactions in MOF nanocrystals. The structure of Co-CAU-36, which is linked through Ni-4-tetraphosphonophenylporphyrins (Ni-H4TPPP4−), was ab initio determined by merging eight high-resolution cRED datasets (0.83–1.00 Å).41 Importantly, to reduce the beam damage and prevent the removal of guest molecules by vacuum, the sample was cooled at 96 K before being transferred to the TEM. The positions of all non-H atoms in the framework as well as most atoms of the guest solvent 1,4-diazabicyclo[2.2.2]octane (DABCO) and H2O molecules in the pores could be located directly from the structure solution using SIR-2014. The remaining atoms were found from subsequent structure refinements. Furthermore, strong hydrogen bonds were identified between the solvent molecules and the framework. It is worth noting that the strong hydrogen bonds are essential to keep the guest molecules stable and ordered in the pores. Importantly, the structural model could be obtained independently from each of the eight cRED datasets. The atomic positions only deviate on average by 0.03(2) Å, for the framework atoms, 0.10(6) Å for the DABCO molecules and by 0.16(11) Å for the water molecules. This indicates a high reproducibility of structure determination by 3DED.
The importance of 3DED for structure determination of low crystalline MOFs is demonstrated on UU-100,92 a Zr-MOF synthesized using a widely studied molecular hydrogen evolution catalyst (HEC) cobaloxime as the organic linker and Zr6 cluster as the inorganic building unit. The PXRD pattern shows broad peaks, indicating a low crystallinity of UU-100 (Fig. 11a). The cRED data contain both sharp spots and streaks, indicating the presence of structural disorders (Fig. 11b). Although the unit cell parameters and space group could be determined directly from the cRED data, the data resolution was too low (ca. 2 Å along the c*-axis and 4 Å along a*- and b*-axes) for ab initio structure solution by direct methods. Instead, the Patterson method was applied to locate the positions of strongly scattering Zr6 clusters. Consequently, the structural model of UU-100 was obtained based on the location of the Zr6 clusters, the unit cell dimensions, the space group, and the chemical linkage information (Fig. 11c). The PXRD pattern was used to further confirm the unit cell and space group (Fig. 11a), while simulated ED patterns were compared to the experimental patterns to confirm the structural model. The structure revealed the spatial proximity of the redox-active sites, which promotes high electron hopping rates, and is crucial for applying UU-100 as electrocatalysts.
![]() | ||
| Fig. 11 (a) Pawley fitting profiles of the PXRD pattern of UU-100 (λ = 1.5418 Å). (b) 3D reciprocal lattice of UU-100 reconstructed from 3DED data. (c) Structure of the cobaloxime linker, and the structural model of UU-100 showing open reticular pore structures.92 Reproduced from ref. 92 with permission from the American Chemical Society, copyright 2019. | ||
Obtaining structural information is particularly important for mesoporous MOFs, whose large structures often result in low crystallinity and small crystal sizes. Two highly stable mesoporous MOFs, PCN-777 and PCN-333 are among such examples.90,91 Both PCN-777 and PCN-333 were designed and synthesized based on reticular chemistry. The same 4,4′,4′′-s-triazine-2,4,6-triyl-tribenzoate (TATB) was used as the organic linker, while different metal clusters, Zr6 clusters for PCN-777 and Al3 clusters for PCN-333, were used as inorganic building units. Their unit cells and space groups were determined from the RED data, as a = 55.57 Å and F
3m for PCN-777, and a = 126.42 Å and Fd
m for PCN-333. The structural models were obtained based on the possible topologies that are compatible with the building units. The revealed large cavities, 3.8 nm for PCN-777 and 1.1, 3.4 and 5.5 nm for PCN-333, were further utilized to accommodate large molecules and enzymes for catalysis. With an improved data collection method such as cRED, it is now possible to obtain high-resolution data and perform ab initio structure determination of mesoporous MOFs, as demonstrated recently on mesoporous Ni-MOF BUT-3394 and Zr-MOF CAU-45.95
The first single crystal COF structure determined by 3DED was COF-320 (Fig. 12a).98 RED data were collected at low temperature (LT, 89 K) to reduce the beam damage and improve the data quality. The structure of COF-320 with 9-fold interpenetration was directly determined, with tetra-(4-anilyl)methane (TAM) and 4,4′-biphenyldialdehyde (BPDA) connected in a dia framework. The resolution of RED data was lower when data were acquired at room temperature (RT). Nevertheless, based on the key structural information extracted from 3DED, it was possible to model the structure. The RT structure of COF-320 has the same connectivity and degree of interpenetration as the LT structure, but with different pore sizes.
![]() | ||
| Fig. 12 (a) COF-320, with different structures at 298 K and 89 K, which were both determined by RED.98 (b) A series of 3D-TPB-COFs, whose structures were determined by cRED.99 Reproduced from ref. 98 and 99 with permission from the American Chemical Society, and the John Wiley & Sons, Inc., respectively copyright 2013 and 2019. | ||
Continuous rotation electron diffraction data collection methods, such as cRED, can provide higher data resolution than RED. This has been demonstrated on three isostructural 3D-TPB-COFs built by tetra(p-aminophenyl)methane (TAPM) and 1,2,4,5-tetraphenylbenzene (TPB) with -H, -Me, or -F substituents, respectively. cRED data of all three COFs have a resolution better than 1.0 Å, allowing an accurate ab initio structure determination where all non-H atoms can be found directly.99 The structure solutions show that the 3D-TPB-COFs have the same topology with a five-fold interpenetrated pts net (Fig. 12b).
Despite the meritorious examples, development and application of 3DED techniques are still in the dawn. More efforts are needed to further develop 3DED techniques and explore their applications in porous materials. While the demands on TEM expertise have been considerably reduced with the establishment of practical protocols for 3DED data collection, knowledge on both TEM and crystallography is still required to collect and process 3DED data. Full automation of the data collection and processing procedures would be beneficial for non-experts to use 3DED for structural analysis. SerialRED43 has been developed as a prototype for such purposes, yet many improvements are desirable. Although 3DED can provide accurate structural information regarding the frameworks, as well as guest molecules, dynamical effects caused by multiple scattering are still a serious concern. To uncover more structural details, such as locations of OSDAs, doped heteroatoms, and hydrogen atoms,42 as well as absolute structures,102 multiple scattering needs to be properly addressed. Considering the importance of dynamic behaviors in porous crystalline materials, further development on in situ 3DED is crucial to gain the fundamental understanding on this aspect. Another remaining challenge is to prevent the loss of crystallinity caused by vacuum and electron beam damage.
In conclusion, structural studies of porous crystalline materials at atomic levels are indispensable to fundamentally understand their physical and chemical properties, which in turn is important for the development of new materials. As TEMs are widely available in laboratories around the world compared to synchrotron facilities, we foresee that the impacts of 3DED for structural analysis will continue to increase, which will certainly accelerate research in the fields of porous materials.
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: 10.1039/d0sc05731b |
| This journal is © The Royal Society of Chemistry 2021 |