Louis R.
Redfern
and
Omar K.
Farha
*
International Institute of Nanotechnology, Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, USA. E-mail: o-farha@northwestern.edu
First published on 17th October 2019
As the field of metal–organic frameworks (MOFs) continues to grow, the physical stability and mechanical properties of these porous materials has become a topic of great interest. While strategies for synthesizing MOFs with desirable chemical functionalities or pore sizes have been established over the past twenty years, design principles to modulate the response of MOFs to mechanical stress are still underdeveloped. The inherent porosity of these frameworks results in many interesting and sometimes unexpected phenomena upon exposure to elevated pressures and other physical stimuli. Beyond its fundamental importance, an understanding of mechanical properties (e.g. bulk modulus, shear modulus, Young's modulus, linear compressibility, and Poisson's ratio) plays an essential role in the post-synthetic processing of MOFs, which has implications in the successful transition of these materials from academic interest to industrial relevance. This perspective provides a concise overview of the efforts to understand the mechanical properties of MOFs through experimental and computational methods. Additionally, current limitations and possible future directions for the field are also discussed briefly.
The compression of porous materials has long been a fascinating topic of study, as the “empty” space in the structures yields surprising and unexpected behaviour at high pressures. Early efforts in this area investigated zeolites under high-pressure conditions.6–10 While zeolites and MOFs are notably distinct in many ways, prior analysis of the former helps to provide insight into the response of the latter to mechanical stress. Still, the inclusion of structural organic linkers in MOFs precludes the direct transfer of conclusions regarding zeolites to these hybrid materials. Likewise, although the fundamental advances in understanding the mechanical properties of MOFs is certain to influence the study of new porous materials (e.g. covalent organic frameworks), each class demands proper investigation.
The rich structural diversity of MOFs provides a near limitless number of interesting species to examine, but it comes with a cost in that the sheer number of distinct frameworks cannot be feasibly measured for every property. Given this limitation, it is important to consider structure–property relationships that can enable a more rapid and heuristic evaluation of materials design. Such relationships can only be drawn from careful, thorough studies that systematically vary a single structural characteristic of a MOF. Unfortunately, experimental data are often complicated by subtle differences between MOFs resulting from batch-to-batch variance (e.g. defect density, porosity, and guest loading). Moreover, many experiments rely on advanced techniques such as diamond anvil cell (DAC) sample environments and synchrotron radiation which puts significant constraints on the availability of such experimental data. Fortunately, computational simulation is a powerful tool for studying the mechanical properties of MOFs and has been used extensively in stand-alone and joint investigations.11,12 These two complementary approaches are essential for the successful analysis of the mechanical properties of MOFs.
The examples highlighted in this review are selected to summarize many of the major advances and efforts toward understanding the mechanical properties of MOFs over the past decade. The linker, node, and structure of each MOF discussed herein is displayed in Table 1. Several reviews on this topic have been written that delve into the fine details of individual investigations;13–16 however, this perspective is intended to provide a broad context of the field while discussing in depth several reports that stand out in the field. Furthermore, this work is not meant to be a tutorial for conducting high-pressure experiments, as this is covered extensively elsewhere.17 The studies discussed herein are organized by into broad classes of MOFs, emphasizing the often varied results obtained from subtle differences in experimental design, even when examining closely related or identical materials.
MOF | Organic linker | Inorganic node | Structure | Ref. |
---|---|---|---|---|
ZIF-3 |
![]() |
![]() |
![]() |
25 |
ZIF-4 |
![]() |
![]() |
![]() |
22–24 and 27–29 |
ZIF-8 |
![]() |
![]() |
![]() |
20–22, 24, 26, 27 and 35 |
ZIF-62 |
![]() |
![]() |
![]() |
29 |
HKUST-1 (Cu-BTC) |
![]() |
![]() |
![]() |
37–42 |
MOF-5 (IRMOF-1) |
![]() |
![]() |
![]() |
43–47 |
MIL-53(M) |
![]() |
![]() |
![]() |
48–54 |
UiO-66 |
![]() |
![]() |
![]() |
16, 40, 57–62 and 67 |
UiO-67 |
![]() |
![]() |
![]() |
63 and 67 |
UiO-abdc |
![]() |
![]() |
![]() |
63 |
Zr6O4(OH)4(edb)6 (1) |
![]() |
![]() |
![]() |
64 |
Zr6O4(OH)4(sdc)6 (2) |
![]() |
![]() |
![]() |
64 |
MIL-140 |
![]() |
![]() |
![]() |
65 and 66 |
NU-901 |
![]() |
![]() |
![]() |
67 |
Co2(4,4′-bpy)3(NO3)4 (3) |
![]() |
![]() |
![]() |
70 |
Zn2(L)2(dabco) (4) |
![]() |
![]() |
![]() |
71 |
MOF-74 |
![]() |
![]() |
![]() |
72 and 73 |
Cu24 isophthalate (L1/L2) (5) |
![]() |
![]() |
![]() |
74 |
DUT-60 |
![]() |
![]() |
![]() |
75 |
Sc2(BDC)3 |
![]() |
![]() |
![]() |
77 |
MOF-520 |
![]() |
![]() |
![]() |
80 |
Since the first reports of DAC diffraction experiments of ZIF-8 [Zn(2-methylimidazolate)2] in 2009,20,21 many groups have examined this diverse class of MOFs to understand their behaviour at high pressure and temperature. Cheetham and coworkers found that across seven distinct materials, the elastic modulus and hardness tend to decrease as solvent accessible volume increases,22 and reported the high-pressure phase of ZIF-4 [Zn(imidazolate)2] following a single-crystal X-ray diffraction (SCXRD) study.23 In addition to experimental efforts, computational analysis has provided invaluable insight into the mechanism and implications of high pressure on MOFs. Coudert and coworkers probed the compression of ZIF-4 and ZIF-8, finding that amorphization of these MOFs could be explained by invoking a “shear-mode softening” mechanism.24 Ryder and Tan sought to isolate the role of topology in the mechanical properties of a series of MOFs with identical chemical components. Their findings demonstrate that the spatial orientation of nodes and linkers can have a significant impact on the stability of MOFs, as evidenced by the exceptionally low Young's and shear moduli of ZIF-3 [Zn(imidazolate)2].25 Modifications to the electronic properties26 and steric bulk27 of the organic linkers have been shown to enhance or diminish the various mechanical properties of the porous frameworks. Recent research interest in the concept of liquid- and glass-phase MOFs has brought about redoubled efforts to map out the high-pressure–high-temperature phase space of ZIF-4 (ref. 28) and ZIF-62 [Zn(benzimidazolate)0.25(imidazolate)1.75].29
Among the first experimental investigations of ZIFs under high pressures, Moggach, Bennett, and Cheetham conducted a SCXRD study on ZIF-8 in a DAC.20 High-pressure diffraction data were collected with a mixture of methanol and ethanol as a pressure transmitting fluid. Upon raising the pressure to 0.18 GPa, the crystal structure swelled unexpectedly rather than compressing. This counterintuitive phenomenon of unit cell expansion under pressure is now found throughout the literature and is attributed to “hyperfilling” of the MOF pores with the hydrostatic fluid, enlarging the framework at moderate pressures. At 1.47 GPa, the ZIF-8 sample underwent a single-crystal to single-crystal phase transition in which the imidazolate ligands twist and increase the accessible pore volume (Fig. 1). Shortly after this report, Chapman, Halder, and Chupas published their own study utilizing powder X-ray diffraction (PXRD) of ZIF-8 in a DAC using Fluorinert™ FC-75 as a nonpenetrating pressure transmitting fluid. Under these conditions, the bulk modulus (K0) was estimated to be 6.5(4) GPa, with clear amorphization of the material at pressures beyond 0.34 GPa.21 In addition to diffraction studies, Chapman et al. monitored the porosity of pressure-treated ZIF-8 using nitrogen adsorption, revealing a precipitous drop in surface area and pore volume following exposure to pressure. The stark contrast between the behaviour of ZIF-8 in the studies of Chupas and Cheetham is due to the choice of hydrostatic fluid—FC-75 is too large to fill the pores of ZIF-8 and support the framework, leading to rapid amorphization.
![]() | ||
Fig. 1 Accessible pore volume of ZIF-8 at (a) 0 GPa and (b) 1.47 GPa. The twist of the imidazolate ligands opens new channels of porosity within the structure. Adapted with permission from ref. 20. Copyright 2009, Wiley-VCH. |
While experimental studies are essential for understanding the pressure response of MOFs, the insight that is gained from computational simulations is crucial for understanding the mechanism of compression and deformation. In 2013, Coudert and coworkers shed light on the mechanism of ZIF-8 and ZIF-4 amorphization through a molecular dynamics study.24 This investigation monitored the elastic constants that describe the mechanical properties of ZIF-8 as a function of pressure. Notably, the elastic constant corresponding to the shear modulus (C44) rapidly drops as the pressure is increased from 0.0 GPa to 0.35 GPa (Fig. 2), indicating that ZIF-8 becomes highly susceptible to shear forces as hydrostatic pressure increases. This behaviour is known as shear-mode softening and leads to the instability of the framework at pressures above 0.4 GPa. This study also indicates that ZIF-4 displays similar shear-mode softening, suggesting that this mechanism of amorphization may be generalizable to other ZIFs and MOFs.
![]() | ||
Fig. 2 Plot of ZIF-8 unit cell parameter and elastic constants vs. pressure. C11, C12, and C44, correspond to the moduli for axial compression, dilation upon compression, and shear, respectively. The precipitous decline in C44 indicates a drop in the shear modulus at elevated pressure, a phenomenon known as “shear-mode softening”. Adapted with permission from ref. 24. Copyright 2013, American Chemical Society. |
As the field of MOFs continues to advance, attention has been drawn to non-crystalline, glass30–32 and liquid33,34 phases of MOFs. In this area, Bennett and coworkers have explored the high-pressure and high-temperature phase diagram of ZIF-62 and ZIF-4 with both experimental and computational methods.29 Remarkably, this study reveals that the melting point of both ZIFs is lowered significantly at elevated hydrostatic pressures, opening the door to possible synthetic strategies to achieve liquid and glass MOFs that normally decompose before melting at ambient pressure. Further investigations into the high-pressure and high-temperature behaviour of ZIF-4 from Bennett and coworkers reveals four distinct high-pressure–high-temperature crystalline phases.28 The crystal structure of two of the phases were determined by PXRD refinement techniques, while structural refinement beyond space group and unit cell assignment was intractable for the other two phases. The phase diagram of ZIF-4 derived from this study is strikingly complex (Fig. 3), emphasizing the need for thorough studies to outline and understand the polymorphism of these frameworks. These exciting findings invite researchers to probe the high-pressure and high-temperature space for other classes of MOFs.
![]() | ||
Fig. 3 Pressure–temperature phase diagram of ZIF-4. The complexity of the diagram highlights the need for thorough studies to elucidate the numerous phases of MOFs. Reproduced with permission from ref. 28. Copyright 2019, American Chemical Society. |
The examples highlighted here emphasize several pervasive characteristics of high-pressure studies of MOFs: (1) seemingly minor experimental details (e.g. choice of hydrostatic fluid) can have a drastic impact on results, (2) the combination of experiment and simulation is vital to build a deeper understanding of the processes at play at high pressure, and (3) the field of high-pressure–high-temperature behaviour of MOFs contains boundless space to explore in future efforts. As MOF technologies continue to advance and ZIFs become incorporated into composite materials, such as encapsulating enzymes and polymer fibres, the mechanical properties of these advanced materials will become a topic of interest. Moreover, the behaviour of ZIFs under dynamic pressures (e.g. shockwave or impact) has already gained attention35 and will continue to be an interesting area to examine for both fundamental understanding and practical applications.
In the first high-pressure diffraction study of HKUST-1, Chapman, Halder, and Chupas probed the compression of the MOF using PXRD in a DAC in the presence of different pressure-transmitting fluids.38 This important study first demonstrated the now well-known phenomenon of pressure-induced hyperfilling of MOF pores with sufficiently small fluids. For small-molecule fluids, such as isopropanol and methanol/ethanol/water mixtures, HKUST-1 exhibits two distinct regimes of compression: at lower pressures the apparent bulk modulus is exceptionally high (K0 = 117.6(1) GPa), while at higher pressure compression occurs at a more reasonably expected rate, with bulk moduli between 25.9(5) and 41.9(4) GPa. Conversely, when a non-penetrating pressure transmitting fluid (in this case, Fluorinert™ FC-70) is used, the observed bulk modulus is 29.5(7) GPa throughout the hydrostatic limit of the fluid. The plot of unit cell volume vs. pressure (Fig. 4) exemplifies the difference between pressure transmitting fluids. This work demonstrates the need for appropriate hydrostatic fluid selection to accurately determine the bulk modulus of the intrinsic framework.
![]() | ||
Fig. 4 Changes in the unit cell volume of HKUST-1 as a function of pressure in the presence of different pressure transmitting fluids. The different initial slopes demonstrate the importance of non-penetrating fluids for the accurate determination of bulk moduli in MOFs. Adapted with permission from ref. 38. Copyright 2008, American Chemical Society. |
While Moggach and coworkers again observed pore hyperfilling in the presence of N,N-diethylformamide during SCXRD experiments in a DAC,46 Hu and Zhang found that MOF-5 undergoes rapid amorphization at pressures as low as 3.5 MPa in the absence of pressure-transmitting fluid.47 In this study, a number of ex situ experiments were conducted after exposure of MOF-5 to pressures ranging from 0.0 to 10.3 MPa. The loss of crystallinity observed by PXRD was corroborated by a steady decrease in the surface area as the pressure treatment increased (Fig. 5), indicating pressure-induced pore collapse. Furthermore, Raman spectroscopy revealed significant changes in the relative intensity of peaks after treatment with 10.3 MPa, further supporting structural changes under pressure. These results highlight the importance of pressure-transmitting fluid in studies of MOF compression, as the behaviour of these frameworks is highly dependent on the surrounding environment and the presence of fluid in the pores.
![]() | ||
Fig. 5 Surface area of MOF-5 as a function of pressure. Elevated pressures were obtained using an anvil press in the absence of pressure transmitting fluid. Surface area values were determined using N2 adsorption measurements at 77 K. Adapted with permission from ref. 47. Copyright 2010 by American Physical Society. |
Wine-rack type MOFs have been predicted to exhibit interesting anisotropy in their elastic properties, and, in some cases, negative linear compressibility (i.e. expansion of a structure along one axis upon application of pressure) is expected to arise due to this anisotropy.53 This phenomenon was demonstrated experimentally by Serra-Crespo et al. in an X-ray diffraction study examining MIL-53(Al) and MIL-53-NH2(Al) in a DAC.54 While the unit cell volume decreases as expected upon increasing pressure, the individual unit cell parameters reveal interesting behaviour in both MOFs. Lattice parameters corresponding to the a and c crystallographic axes continuously decrease as a function of pressure, while the b axis, corresponding to the long axis of the diamond-shaped pores (Fig. 6), increases at low pressures (up to ∼3 GPa). This example of negative linear compressibility is an important step in understanding the mechanism of compression in MOFs and highlights the importance of exploring these phenomena further.
![]() | ||
Fig. 6 (a) Unit cell volume of MIL-53(Al), (b) unit cell parameters of MIL-53(Al), and (c) unit cell parameters of NH2-MIL-53(Al) as a function of pressure. The parameter corresponding to the b-axis exhibits a notable increase at low pressures—an example of negative linear compression. Adapted from ref. 54 – published by the Royal Society of Chemistry. |
Experimental evidence of the impact of defects on the compression of UiO-66 was reported by Dissegna et al., corroborating Coudert's findings that increasing defect density generally yields a lower bulk modulus.62 By conducting these experiments in a water-filled cell rather than a traditional DAC, the pressure control at low pressures is exceptionally precise, allowing for measurements in 25 MPa intervals from 0–4 GPa. As the defect density increases from ∼3% to ∼26% of linkers missing, the bulk modulus decreases as expected; however, at higher defect densities (∼28% missing linkers), the bulk modulus of UiO-66 appears to increase (Fig. 7). This result comes as a surprise, inviting further experimental and computational study to understand how the presence of more extensive defects impacts the mechanical properties of MOFs.
![]() | ||
Fig. 7 Unit cell volume of UiO-66 as a function of pressure. Each sample exhibits different defect density. A more shallow slope indicates a higher bulk modulus, while steeper slopes correspond to lower bulk modulus values. Adapted with permission from ref. 62. Copyright 2019, American Chemical Society. |
In the Farha group's first venture into the field of MOFs under pressure, we conducted a systematic study into the compression of two topological families of Zr-MOFs: UiO-type MOFs (fcu topology) and the NU-900 series of wine-rack type MOFs (scu topology).67 The aim of this study was to experimentally probe the effects of linker length and porosity across different topologies, allowing for more broadly generalizable conclusions to be drawn from the data. In general, we found that while linker length and nearest-node distance correlate well with bulk modulus for each series, only void fraction serves as a good predictor for both families of MOFs studied. Interestingly, two samples exhibited significantly lower bulk moduli than expected based on the compressibility of similar materials and previous computational results. Careful analysis of ambient condition single crystal structure models and unit cell volume changes reveal that the linkers of these two samples are “pre-distorted” and thereby do not resist compression as well as undistorted analogues (Fig. 8). It is worth noting that experimental surveys of the tremendous structural diversity of Zr-MOFs using DACs is currently hampered by the lack of appropriate non-penetrating pressure transmitting fluids for the numerous mesoporous frameworks. In order to properly determine the bulk modulus of these highly porous materials, the development of new pressure transmitting fluids is essential.
![]() | ||
Fig. 8 (a) MOFs with undistorted structures are more resistant to compression than (b) those with distortions in their structures at ambient pressure. Adapted with permission from ref. 67. Copyright 2019, American Chemical Society. |
While Zn(CN)2 normally forms a dense, non-porous structure, Lapidus et al. demonstrated that at elevated pressures (∼1 GPa) the crystal structure can undergo phase transitions to yield an open framework with lower density.69 Interestingly, these phase transitions depend greatly on the pressure transmitting fluid, emphasizing the role of the fluid in the mechanism of the transition. In a study of a 4,4′-bipyridine MOF with Co nodes (3), Zhou et al. observed a phase transition followed by negative linear compression.70 These studies are an example of the power of pressure to alter MOF structures and imbue new and desirable properties on the material.
Fundamental understanding of the mechanical properties of MOFs has led many to ask how these properties can be modified and controlled to suit the needs of different applications. Great effort has thus been put toward developing post-synthetic modifications that modulate the mechanical stability of MOFs. Henke, Li, and Cheetham demonstrate the anisotropy and guest-dependence of the elastic modulus and hardness of MOFs (4) pillared with DABCO (diazabicyclo[2.2.2]octane), indicating that simply modifying solvent choice can alter these properties.71 Several efforts have sought to improve the mechanical properties of MOF-74 by grafting the framework onto graphene or including N,N′-dimethylethylenediamine guest molecules in the pores.72,73 While the gas sorption properties remained fairly similar, these post-synthetic processes have an impact on the elastic modulus, hardness, Young's modulus, and shear modulus.
In the ever-expanding field of porous materials, new and fascinating materials are continuously created and characterized. Recently, Lal et al. synthetized an interesting MOF-like material (5) by cross-linking the dangling ligands on metal–organic polyhedra.74 While this structure is far from a traditional, crystalline MOF, nanoindentation studies reveal that its hardness correlates with the crosslinking density, akin to trends seen in other MOFs with additional supporting linkers between metal nodes. As new materials are developed, the foundational knowledge gained through years of studying mechanical properties of MOFs and other porous materials will help to guide and direct the design of these structures.
Among the most interesting phenomena observed in the articles discussed above are pressure-induced phase transitions. While ball milling has been a long-established method for MOF synthesis,76 the application of pressure to achieve new phases of these scaffolds presents largely unexplored synthetic opportunities. Recently, Moggach and coworkers observed structural changes to a Sc-based MOF (Sc2BDC3) that alters its pore shape to accommodate molecules (e.g. 2-methylbutane) that cannot fit inside the MOF under ambient conditions (Fig. 9).77 This, along with the introduction of porosity into Zn(CN),69 are shining examples of the power of pressure to achieve new structures with desirable properties, rather than acting only as a destructive force that results in amorphization. Given the propensity for pore hyperfilling in MOFs at elevated pressure, we envision the potential for controlled phase transitions from MOFs comprised of 2-D sheets to 3-D structures and vice versa. Understanding the role of pressure in the formation of new MOF phases will enable the synthesis of novel materials that cannot be obtained with traditional solvothermal synthesis.
![]() | ||
Fig. 9 Pore shape modification of Sc2BDC3 upon exposure to 2-methylbutane at elevated pressure. Rotation of the phenyl rings reveals a larger cavity that can accommodate molecules that cannot fit in the pores at ambient conditions. Reprinted with permission from ref. 77. Copyright 2015, American Chemical Society. |
While the flexibility of the organic linker in MOFs has been shown to correlate with compression of the material, recent synthetic efforts have revealed numerous examples of frameworks with inherent flexibility under mild stimuli, such as changing solvent polarity.78,79 Investigating the mechanical properties of these highly flexible frameworks is crucial to understand the differences between these pliable scaffolds and more rigid MOFs. We envision two strategies for designing materials that can withstand high pressures: reinforce the structure to resist compression, or intentionally include a “bend, don't break” motif that allows the MOF to deform under pressure and rebound to the original structure upon release of the pressure. In an example of structural reinforcement, Yaghi and coworkers recently demonstrated the viability of “molecular retrofitting” to improve the stability of a framework that readily undergoes amorphization under pressure.80 By installing 4,4′-biphenyldicarboxylate into MOF-520 [Al8(μ-OH)8(HCOO)4(1,3,5-benzenetribenzoate)4], crystallinity was maintained up to 5.5 GPa, while the pristine MOF amorphized at pressures over 2.8 GPa (Fig. 10). On the other hand, the capacity for flexibility to impart structural resilience and maintain MOF functionality requires further study.
![]() | ||
Fig. 10 MOF-520 (left) demonstrates improved stability under pressure after incorporation of 4,4′-biphenyldicarboxylic acid “girder” linkers (right). Adapted with permission from ref. 80. Copyright 2017, American Chemical Society. |
This journal is © The Royal Society of Chemistry 2019 |