Thermodynamic Insight in the High-Pressure Behavior of UiO-66: Effect of Linker Defects and Linker Expansion

Plasticity of Metal-Organic Framework Crystals: Thermally Activated Collapse of Nanopores

Unrecoverable deformation or plasticity can be generated in crystalline metal-organic frameworks (MOFs) by compressive loading with different rates in various applications. Herein, plastic behaviors of MOF HKUST-1 crystals are investigated by a series of in situ strain-rate-dependent compression tests. The yield strength is found to significantly increase with increasing strain rate, following a logarithmic dependence. Our reactive molecular dynamics simulations illustrate that the yielding of crystalline HKUST-1 is induced by the irreversible collapse of its nanopores, which can be accelerated by thermal activation at finite temperatures. Based on this mechanism together with the reaction rate theory, we derive an analytical expression relating the yield strength of MOFs and strain rate, which fits experimental findings well. Overall, this work can expand our current understanding of MOF plasticity, which is of importance for the mechanical shaping and various applications of MOF crystals.

An automated protocol to construct flexibility parameters for classical forcefields: applications to metal-organic frameworks

In this work, forcefield flexibility parameters were constructed and validated for more than 100 metal-organic frameworks (MOFs). We used atom typing to identify bond types, angle types, and dihedral types associated with bond stretches, angle bends, dihedral torsions, and other flexibility interactions. Our work used Manz's angle-bending and dihedral-torsion model potentials. For a crystal structure containing N atoms in its unit cell, the number of independent flexibility interactions is 3(N atoms - 1). Because the number of bonds, angles, and dihedrals is normally much larger than 3(N atoms - 1), these internal coordinates are redundant. To reduce (but not eliminate) this redundancy, our protocol prunes dihedral types in a way that preserves symmetry equivalency. Next, each dihedral type is classified as non-rotatable, hindered, rotatable, or linear. We introduce a smart selection method that identifies which particular torsion modes are important for each rotatable dihedral type. Then, we computed the force constants for all flexibility interactions together via LASSO regression (i.e., regularized linear least-squares fitting) of the training dataset. LASSO automatically identifies and removes unimportant forcefield interactions. For each MOF, the reference dataset was quantum-mechanically-computed in VASP via DFT with dispersion and included: (i) finite-displacement calculations along every independent atom translation mode, (ii) geometries randomly sampled via ab initio molecular dynamics (AIMD), (iii) the optimized ground-state geometry using experimental lattice parameters, and (iv) rigid torsion scans for each rotatable dihedral type. After training, the flexibility model was validated across geometries that were not part of the training dataset. For each MOF, we computed the goodness of fit (R-squared value) and the root-mean-squared error (RMSE) separately for the training and validation datasets. We compared flexibility models with and without bond-bond cross terms. Even without cross terms, the model yielded R-squared values of 0.910 (avg across all MOFs) ± 0.018 (st. dev.) for atom-in-material forces in the validation datasets. Our SAVESTEPS protocol should find widespread applications to parameterize flexible forcefields for material datasets. We performed molecular dynamics simulations using these flexibility parameters to compute heat capacities and thermal expansion coefficients for two MOFs.

Impact of Vertex Functionalization on Flexibility of Porous Organic Cages

Efficient carbon capture requires engineered porous systems that selectively capture CO2 and have low energy regeneration pathways. Porous liquids (PLs), solvent-based systems containing permanent porosity through the incorporation of a porous host, increase the CO2 adsorption capacity. A proposed mechanism of PL regeneration is the application of isostatic pressure in which the dissolved nanoporous host is compressed to alter the stability of gases in the internal pore. This regeneration mechanism relies on the flexibility of the porous host, which can be evaluated through molecular simulations. Here, the flexibility of porous organic cages (POCs) as representative porous hosts was evaluated, during which pore windows decreased by 10-40% at 6 GPa. POCs with sterically smaller functional groups, such as the 1,2-ethane in the CC1 POC resulted in greater imine cage flexibility relative to those with sterically larger functional groups, such as the cyclohexane in the CC3 POC that protected the imine cage from the application of pressure. Structural changes in the POC also caused CO2 adsorption to be thermodynamically unfavorable beginning at ∼2.2 GPa in the CC1 POC, ∼1.1 GPa in the CC3 POC, and ∼1.0 GPa in the CC13 POC, indicating that the CO2 would be expelled from the POC at or above these pressures. Energy barriers for CO2 desorption from inside the POC varied based on the geometry of the pore window and all the POCs had at least one pore window with a sufficiently low energy barrier to allow for CO2 desorption under ambient temperatures. The results identified that flexibility of the CC1, CC3, or CC13 POCs under compression can result in the expulsion of captured gas molecules.

Node Distortions in UiO-66 Inform Negative Thermal Expansion Mechanisms: Kinetic Effects, Frustration, and Lattice Hysteresis

In metal-organic frameworks (MOFs) the interplay between the dynamics of individual components and how these are constrained by the extended lattice can yield unusual emergent phenomena. For the archetypal Zr-MOF, UiO-66, we explore the cooperative dynamics of a Zr-node transformation that gives rise to negative thermal expansion (NTE). Using in situ synchrotron X-ray scattering, with powder diffraction and pair distribution function (PDF) analyses, we identify lattice hysteresis and a thermal ramp-rate-dependence of the thermal expansion. Specifically, kinetic trapping of distorted node states formed at high temperature, leads to broad variability in the apparent thermal expansion which ranges from large positive to large negative thermal expansion with coefficients of thermal expansion (CTE) from +45 to -80 × 10-6K-1. Time-resolved relaxation studies at selected temperatures suggest that when equilibrated UiO-66 is intrinsically NTE, with a CTE of -35 × 10-6K-1. Kinetic trapping of the node-distorted state following high temperature activation has broad implications for characterization and applications of these Zr-MOFs; the nonequilibrium node state depends on the thermal history of the sample with quench vs slow cooling likely to impact gas binding, pore volume, and accessible catalytic sites.

How Reproducible is the Synthesis of Zr-Porphyrin Metal-Organic Frameworks? An Interlaboratory Study

Metal-organic frameworks (MOFs) are a rapidly growing class of materials that offer great promise in various applications. However, the synthesis remains challenging: for example, a range of crystal structures can often be accessed from the same building blocks, which complicates the phase selectivity. Likewise, the high sensitivity to slight changes in synthesis conditions may cause reproducibility issues. This is crucial, as it hampers the research and commercialization of affected MOFs. Here, it presents the first-ever interlaboratory study of the synthetic reproducibility of two Zr-porphyrin MOFs, PCN-222 and PCN-224, to investigate the scope of this problem. For PCN-222, only one sample out of ten was phase pure and of the correct symmetry, while for PCN-224, three are phase pure, although none of these show the spatial linker order characteristic of PCN-224. Instead, these samples resemble dPCN-224 (disordered PCN-224), which has recently been reported. The variability in thermal behavior, defect content, and surface area of the synthesised samples are also studied. The results have important ramifications for field of metal-organic frameworks and their crystallization, by highlighting the synthetic challenges associated with a multi-variable synthesis space and flat energy landscapes characteristic of MOFs.

A Critical Assessment on Calculating Vibrational Spectra in Nanostructured Materials

Vibrational spectroscopy is an omnipresent spectroscopic technique to characterize functional nanostructured materials such as zeolites, metal-organic frameworks (MOFs), and metal-halide perovskites (MHPs). The resulting experimental spectra are usually complex, with both low-frequency framework modes and high-frequency functional group vibrations. Therefore, theoretically calculated spectra are often an essential element to elucidate the vibrational fingerprint. In principle, there are two possible approaches to calculate vibrational spectra: (i) a static approach that approximates the potential energy surface (PES) as a set of independent harmonic oscillators and (ii) a dynamic approach that explicitly samples the PES around equilibrium by integrating Newton's equations of motions. The dynamic approach considers anharmonic and temperature effects and provides a more genuine representation of materials at true operating conditions; however, such simulations come at a substantially increased computational cost. This is certainly true when forces and energy evaluations are performed at the quantum mechanical level. Molecular dynamics (MD) techniques have become more established within the field of computational chemistry. Yet, for the prediction of infrared (IR) and Raman spectra of nanostructured materials, their usage has been less explored and remain restricted to some isolated successes. Therefore, it is currently not a priori clear which methodology should be used to accurately predict vibrational spectra for a given system. A comprehensive comparative study between various theoretical methods and experimental spectra for a broad set of nanostructured materials is so far lacking. To fill this gap, we herein present a concise overview on which methodology is suited to accurately predict vibrational spectra for a broad range of nanostructured materials and formulate a series of theoretical guidelines to this purpose. To this end, four different case studies are considered, each treating a particular material aspect, namely breathing in flexible MOFs, characterization of defects in the rigid MOF UiO-66, anharmonic vibrations in the metal-halide perovskite CsPbBr3, and guest adsorption on the pores of the zeolite H-SSZ-13. For all four materials, in their guest- and defect-free state and at sufficiently low temperatures, both the static and dynamic approach yield qualitatively similar spectra in agreement with experimental results. When the temperature is increased, the harmonic approximation starts to fail for CsPbBr3 due to the presence of anharmonic phonon modes. Also, the spectroscopic fingerprints of defects and guest species are insufficiently well predicted by a simple harmonic model. Both phenomena flatten the potential energy surface (PES), which facilitates the transitions between metastable states, necessitating dynamic sampling. On the basis of the four case studies treated in this Review, we can propose the following theoretical guidelines to simulate accurate vibrational spectra of functional solid-state materials: (i) For nanostructured crystalline framework materials at low temperature, insights into the lattice dynamics can be obtained using a static approach relying on a few points on the PES and an independent set of harmonic oscillators. (ii) When the material is evaluated at higher temperatures or when additional complexity enters the system, e.g., strong anharmonicity, defects, or guest species, the harmonic regime breaks down and dynamic sampling is required for a correct prediction of the phonon spectrum. These guidelines and their illustrations for prototype material classes can help experimental and theoretical researchers to enhance the knowledge obtained from a lattice dynamics study.

Physicochemical properties and density functional theory calculation of octahedral UiO-66 with Bis(Trifluoromethanesulfonyl)imide ionic liquids

In this study, the physicochemical properties and molecular interactions between zirconium-based metal-organic framework (UiO-66) and three different ionic liquids based on bis(trifluoromethanesulfonyl)imide anion (EMIM+, BMIM+ and OMIM+) was performed via a combined experimental and computational approach. The ionic liquid loaded UiO-66 or IL@UiO-66 was synthesized and characterized to understand the host-guest interaction. Density functional theory calculation was performed to analyse the electronic structure of IL@UiO-66 to provide molecular insight on the dominant interactions occurred in the hybrid material. Results showed that all ILs were successfully incorporated into the micropores of UiO-66. The 3D framework was retained even after loaded with ILs as analyzed from XRD pattern. FTIR spectrum reveals that interactions of ILs with UiO-66 influenced by the alkyl chain length of the cation. The anion has a profound affinity with the UiO-66 due to the presence of electronegative atoms. Phase transition study from DSC suggested that the incorporation of ILs has stabilized the framework of UiO-66 by shifting the endothermic peak to a higher state. These findings were further elaborated with DFT calculation. Geometrical optimizations confirmed the structural parameter changes of UiO-66 when loaded with ILs. These was mainly contributed by the non-covalent interactions which was confirmed by the reduced density gradient scattered plot. Another important findings are the strength of hydrogen bonding at the host-guest interface was influenced by the alkyl chain length. The molecular orbital analysis also shows that the size of alkyl chain influence the reactivity of the hybrid material. The present study provides fundamental insights on the molecular interaction of UiO-66 and ILs as a hybrid material, which can open new possibilities for advanced material for metal-organic framework applications in energy storage system, catalysis, gas storage and medicinal chemistry.

Challenges in modelling dynamic processes in realistic nanostructured materials at operating conditions

The question is addressed in how far current modelling strategies are capable of modelling dynamic phenomena in realistic nanostructured materials at operating conditions. Nanostructured materials used in applications are far from perfect; they possess a broad range of heterogeneities in space and time extending over several orders of magnitude. Spatial heterogeneities from the subnanometre to the micrometre scale in crystal particles with a finite size and specific morphology, impact the material's dynamics. Furthermore, the material's functional behaviour is largely determined by the operating conditions. Currently, there exists a huge length-time scale gap between attainable theoretical length-time scales and experimentally relevant scales. Within this perspective, three key challenges are highlighted within the molecular modelling chain to bridge this length-time scale gap. Methods are needed that enable (i) building structural models for realistic crystal particles having mesoscale dimensions with isolated defects, correlated nanoregions, mesoporosity, internal and external surfaces; (ii) the evaluation of interatomic forces with quantum mechanical accuracy albeit at much lower computational cost than the currently used density functional theory methods and (iii) derivation of the kinetics of phenomena taking place in a multi-length-time scale window to obtain an overall view of the dynamics of the process. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

Calculation of Self, Corrected, and Transport Diffusivities of Isopropyl Alcohol in UiO-66

The UiO-6x family of metal-organic frameworks has been extensively studied for applications in chemical warfare agent (CWA) capture and destruction. An understanding of intrinsic transport phenomena, such as diffusion, is key to understanding experimental results and designing effective materials for CWA capture. However, the relatively large size of CWAs and their simulants makes diffusion in the small-pored pristine UiO-66 very slow and hence impractical to study directly with direct molecular simulations because of the time scales required. We used isopropanol (IPA) as a surrogate for CWAs to investigate the fundamental diffusion mechanisms of a polar molecule within pristine UiO-66. IPA can form hydrogen bonds with the μ3-OH groups bound to the metal oxide clusters in UiO-66, similar to some CWAs, and can be studied by direct molecular dynamics simulations. We report self, corrected, and transport diffusivities of IPA in pristine UiO-66 as a function of loading. Our calculations highlight the importance of the accurate modeling of the hydrogen bonding interactions on diffusivities, with about an order of magnitude decrease in diffusion coefficients when the hydrogen bonding between IPA and the μ3-OH groups is included. We found that a fraction of the IPA molecules have very low mobility during the course of a simulation, while a small fraction are highly mobile, exhibiting mean square displacements far greater than the ensemble average.

UiO-66 framework with an encapsulated spin probe: synthesis and exceptional sensitivity to mechanical pressure

Probes sensitive to mechanical stress are in demand for the analysis of pressure distribution in materials, and the design of pressure sensors based on metal-organic frameworks (MOFs) is highly promising due to their structural tunability. We report a new pressure-sensing material, which is based on the UiO-66 framework with trace amounts of a spin probe (0.03 wt%) encapsulated in cavities. To obtain this material, we developed an approach for encapsulation of stable nitroxide radical TEMPO ((2,2,6,6-tetramethylpiperidin-1-yl)oxyl) into the micropores of UiO-66 during its solvothermal synthesis. Pressure read-out using electron paramagnetic resonance (EPR) spectroscopy allows monitoring the degradation of the defected MOF structure upon pressurization, where full collapse of pores occurs at as low a pressure as 0.13 GPa. The developed methodology can be used in and ex situ and provides sensitive tools for non-destructive mapping of pressure effects in various materials.

Effect of Modulation and Functionalization of UiO-66 Type MOFs on Their Surface Thermodynamic Properties and Lewis Acid–Base Behavior

In this study, we investigated the surface thermodynamic properties of four MOF structures of the UiO-66 series, by employing seven molecular models, a thermal model, and three other methods using the inverse gas chromatography (IGC) technique at infinite dilution. We first determined the effect of the modulation of UiO-66 by an acid (e.g., formic acid and acetic acid) and on the other hand, we studied the effect of the functionalization of the organic linker by an amine group (NH2) on their dispersive component of the surface energy and on their Lewis acid–base properties. We found that all the studied MOFs presented an amphoteric character with a strong acidity whose acidity/basicity ratio is greater than 1 using all the models and methods in IGC. Moreover, the introduction of a modulator such as acetic acid or formic acid in the synthesis of these MOFs increased the number of structural defects and therefore increased the acidity of these MOFs. Similarly, the functionalization of the MOF by the NH2 group leads to an increase in the basicity constant of the functionalized MOF while remaining smaller than their acidity constant. In addition, the use of acids as modulators and amine groups as functional groups resulted in an increase in the dispersive component of the surface energy of the MOFs. Finally, comparing the results obtained by the different models and methods and based on the increasing order of the acidity of each MOF, it was clear that the thermal model resulted in more exact and precise values than the others. Our findings pave the way for the design and development of new acid catalysts based on UiO-66 structures.

Computational Characterization of Zr-Oxide MOFs for Adsorption Applications

Zr-oxide secondary building units construct metal-organic framework (MOF) materials with excellent gas adsorption properties and high mechanical, thermal, and chemical stability. These attributes have led Zr-oxide MOFs to be well-recognized for a wide range of applications, including gas storage and separation, catalysis, as well as healthcare domain. Here, we report structure search methods within the Cambridge Structural Database (CSD) to create a curated subset of 102 Zr-oxide MOFs synthesized to date, bringing a unique record for all researchers working in this area. For the identified structures, we manually corrected the proton topology of hydroxyl and water molecules on the Zr-oxide nodes and characterized their textural properties, Brunauer-Emmett-Teller (BET) area, and topology. Importantly, we performed systematic periodic density functional theory (DFT) calculations comparing 25 different combinations of basis sets and functionals to calculate framework partial atomic charges for use in gas adsorption simulations. Through experimental verification of CO₂ adsorption in selected Zr-oxide MOFs, we demonstrate the sensitivity of CO₂ adsorption predictions at the Henry's regime to the choice of the DFT method for partial charge calculations. We characterized Zr-MOFs for their CO₂ adsorption performance via high-throughput grand canonical Monte Carlo (GCMC) simulations and revealed how the chemistry of the Zr-oxide node could have a significant impact on CO₂ uptake predictions. We found that the maximum CO₂ uptake is obtained for structures with the heat of adsorption values >25 kJ/mol and the largest cavity diameters of ca. 6-7 Å. Finally, we introduced augmented reality (AR) visualizations as a means to bring adsorption phenomena alive in porous adsorbents and to dynamically explore gas adsorption sites in MOFs.

Defective UiO-66-NH2 Functionalized with Stable Superoxide Radicals toward Electrocatalytic Nitrogen Reduction with High Faradaic Efficiency

The electrocatalytic nitrogen reduction reaction (NRR) to NH3 is limited by low Faradaic efficiency (FE). Herein, defective UiO-66-NH2 functionalized with quite stable superoxide radicals (O2•) is developed as a highly active NRR catalyst. The experimental and computational results show that one linker per Zr6 node is missed and two Zr atoms are exposed in the defective UiO-66-NH2. One of the two exposed Zr atoms can stably adsorb O2•, and thus, a Zr-OO• site forms during the preparations without light excitation or postoxidation, while the other Zr atom is activated as an active site. The synergistic effects of the two Zr sites in the defective UiO-66-NH2 suppress hydrogen and hydrazine evolutions considerably. They are as follows: (i) due to repulsion of the proton on the active Zr site and stabilization of the proton on the Zr-OO• site, the active Zr site is unfavorable for the adsorption of the proton with a high energy barrier, which is the HER rate-determining step (RDS); (ii) under the assistance of the OO• of the Zr-OO• site, the first hydrogenation step of *N2 (i.e., NRR RDS) on the active Zr site is promoted; and (iii) relying on the assistance of the OO• of the Zr-OO• site, the continuous hydrogenation of *NH2NH2 to produce NH3 on the active Zr site is spontaneously exothermic, whereas its desorption to hydrazine is blocked. Accordingly, an extremely high FE of ∼85.21% has been realized along with a high yield rate of NH3 (∼52.81 μg h-1 mgcat-1). To the best of our knowledge, it is the highest FE that has been achieved in recent years. Radical scavenging treatment of the defective UiO-66-NH2 and detailed investigations of two categories of control samples further verify the favorable effects of the O2• that closely correlates with the missed linkers on the performance of the NRR to NH3. This work opens a new way toward highly efficient NRR catalysts, i.e., stable radical-activating defective metal-organic frameworks.

Combined Deep Learning and Classical Potential Approach for Modeling Diffusion in UiO-66

Modeling of diffusion of adsorbates through porous materials with atomistic molecular dynamics (MD) can be a challenging task if the flexibility of the adsorbent needs to be included. This is because potentials need to be developed that accurately account for the motion of the adsorbent in response to the presence of adsorbate molecules. In this work, we show that it is possible to use accurate machine learning atomistic potentials for metal-organic frameworks in concert with classical potentials for adsorbates to accurately compute diffusivities though a hybrid potential approach. As a proof-of-concept, we have developed an accurate deep learning potential (DP) for UiO-66, a metal-organic framework, and used this DP to perform hybrid potential simulations, modeling diffusion of neon and xenon through the crystal. The adsorbate-adsorbate interactions were modeled with Lennard-Jones (LJ) potentials, the adsorbent-adsorbent interactions were described by the DP, and the adsorbent-adsorbate interactions used LJ cross-interactions. Thus, our hybrid potential allows for adsorbent-adsorbate interactions with classical potentials but models the response of the adsorbent to the presence of the adsorbate through near-DFT accuracy DPs. This hybrid approach does not require refitting the DP for new adsorbates. We calculated self-diffusion coefficients for Ne in UiO-66 from DFT-MD, our hybrid DP/LJ approach, and from two different classical potentials for UiO-66. Our DP/LJ results are in excellent agreement with DFT-MD. We modeled diffusion of Xe in UiO-66 with DP/LJ and a classical potential. Diffusion of Xe in UiO-66 is about a factor of 30 slower than that of Ne, so it is not computationally feasible to compute Xe diffusion with DFT-MD. Our hybrid DP-classical potential approach can be applied to other MOFs and other adsorbates, making it possible to use an accurate DP generated from DFT simulations of an empty adsorbent in concert with existing classical potentials for adsorbates to model adsorption and diffusion within the porous material, including adsorbate-induced changes to the framework.

Exploration of Hierarchical Metal-Organic Framework as Ultralight, High-Strength Mechanical Metamaterials

Due to the extraordinarily high surface to volume ratio and enormous structural and chemical diversities, metal-organic frameworks (MOFs) have drawn much attention in applications such as heterogeneous catalysis, gas storage separation, and drug delivery, and so on. However, the potential of MOF materials as mechanical metamaterials has not been investigated. In this work, we demonstrated that through the concerted effort of molecular construct and mesoscopic structural design, hierarchical MOFs can exhibit superb mechanical properties. With the cutting-edge in situ transmission and scanning electron microscope (TEM and SEM) techniques, the mechanical properties of hollow UiO-66 octahedron particles were quantitatively studied by compression on individual specimens. Results showed that the yield strength and Young's modulus of the hierarchical porous framework material presented a distinct "smaller is stronger and stiffer" size dependency, and the maximum yield strength and Young's modulus reached 580 ± 55 MPa and 4.3 ± 0.5 GPa, respectively. The specific strengths were measured as 0.15 ± 0.03 to 0.68 ± 0.11 GPa g^-1 cm³, which is comparable to the previously reported state-of-the-art mechanical metamaterials like glassy carbon nanolattices and pyrolytic carbon nanolattices. This work revealed that MOF materials can be made into a new class of low-density, high-strength mechanical metamaterials and provided insight into the mechanical stability of nanoscale MOFs for practical applications.

Validation of two contrasting capturing mechanisms for gaseous formaldehyde between two different types of strong metal-organic framework adsorbents

In this research, the adsorption behavior of formaldehyde (FA) onto two types of metal-organic frameworks (MOFs: MOF-199 [M199] and UiO-66-NH₂ [U6N]) is investigated against changes in the key process variables (e.g., FA partial pressure (0.5-10 Pa), temperature (30-120 °C), and relative humidity (RH: 0%, 50%, and 100%)). The results revealed that the FA adsorption behavior onto both MOFs is exothermic in nature. Besides, their relative dominance for FA uptake varies interactively with the changes in RH and FA partial pressure levels. As the FA levels increase in dry conditions, their breakthrough volumes (BTV (100% BT)) exhibit contrasting trends: The values of U6N decreased noticeably from 5232 and 3792 L·atm·g^-1, while those of M199 increased from 4152 to 5772 L·atm·g^-1. The superiority of U6N over M199 in the lower FA level (at<5 Pa) is supported by the Lewis acid-base interactions with amine groups (U6N) in line with kinetic/isotherm studies. Such superiority is also persistent at higher (10 Pa) FA level under all humid conditions in line with its higher moisture stability. However, in dry conditions, the reversal of relative dominance in which M199 exhibits enhanced efficacy for 10 Pa FA uptake (relative to U6N) should reflect its breathing effects with the potent role of pore-diffusion mechanism. This study offers valuable insights into the construction of tunable adsorbents with enhanced adsorptivity toward key targets.

Zirconium-Based Metal Organic Frameworks for the Capture of Carbon Dioxide and Ethanol Vapour. A Comparative Study

This paper reports on the comparison of three zirconium-based metal organic frameworks (MOFs) for the capture of carbon dioxide and ethanol vapour at ambient conditions. In terms of efficiency, two parameters were evaluated by experimental and modeling means, namely the nature of the ligands and the size of the cavities. We demonstrated that amongst three Zr-based MOFs, MIP-202 has the highest affinity for CO2 (-50 kJ·mol-1 at low coverage against around -20 kJ·mol-1 for MOF-801 and Muc Zr MOF), which could be related to the presence of amino functions borne by its aspartic acid ligands as well as the presence of extra-framework anions. On the other side, regardless of the ligand size, these three materials were able to adsorb similar amounts of carbon dioxide at 1 atm (between 2 and 2.5 µmol·m-2 at 298 K). These experimental findings were consistent with modeling studies, despite chemisorption effects, which could not be taken into consideration by classical Monte Carlo simulations. Ethanol adsorption confirmed these results, higher enthalpies being found at low coverage for the three materials because of stronger van der Waals interactions. Two distinct sorption processes were proposed in the case of MIP-202 to explain the shape of the enthalpic profiles.

Synthesis and Mechanical Properties of sub 5-µm PolyUiO-66 Thin Films on Gold Surfaces

Metal-organic framework (MOF) thin films currently lack the mechanical stability needed for electronic device applications. Polymer-based metal-organic frameworks (polyMOFs) have been suggested to provide mechanical advantages over MOFs, however, the mechanical properties of polyMOFs have not yet been characterized. In this work, we developed a method to synthesize continuous sub-5 µm polyUiO-66(Zr) films on Au substrates, which allowed us to undertake initial mechanical property investigations.  Comparisons between polyUiO-66 and UiO-66 thin films determined polyUiO-66 thin films exhibit a lower modulus but similar hardness to UiO-66 thin films. The initial mechanical characterization indicates that further development is needed to leverage the mechanical property advantages of polyMOFs over MOFs. Additionally, the demonstration in this work of a continuous surface-supported polyUiO-66 thin film also enables utilization of the emerging class of polyMOF materials in sensors and devices applications.

Understanding disorder and linker deficiency in porphyrinic zirconium-based metal-organic frameworks by resolving the Zr8O6 cluster conundrum in PCN-221

Porphyrin-based metal–organic frameworks (MOFs), exemplified by MOF-525, PCN-221, and PCN-224, are promising systems for catalysis, optoelectronics, and solar energy conversion. However, subtle differences between synthetic protocols for these three MOFs give rise to vast discrepancies in purported product outcomes and description of framework topologies. Here, based on a comprehensive synthetic and structural analysis spanning local and long-range length scales, we show that PCN-221 consists of Zr₆O₄(OH)₄ clusters in four distinct orientations within the unit cell, rather than Zr₈O₆ clusters as originally published, and linker vacancies at levels of around 50%, which may form in a locally correlated manner. We propose disordered PCN-224 (dPCN-224) as a unified model to understand PCN-221, MOF-525, and PCN-224 by varying the degree of orientational cluster disorder, linker conformation and vacancies, and cluster–linker binding. Our work thus introduces a new perspective on network topology and disorder in Zr-MOFs and pinpoints the structural variables that direct their functional properties.

Zirconium-based metal–organic frameworks have defective structures that are useful in catalysis and gas storage. Here, the authors study the interplay between cluster disorder and linker vacancies in PCN-221 and propose a new structure model with tilted Zr₆O₄(OH)₄ clusters rather than Zr₈O₆ clusters.

A Computational Study of Isopropyl Alcohol Adsorption and Diffusion in UiO-66 Metal-Organic Framework: The Role of Missing Linker Defect

Defect engineering leads to an effective manipulation of the physical and chemical properties of metal-organic frameworks (MOFs). Taking the common missing linker defect as an example, the defective MOF generally possesses larger pores and a greater surface area/volume ratio, both of which favor an increased amount of adsorption. When it comes to the self-diffusion of adsorbates in MOFs, however, the missing linker is a double-edged sword: the unsaturated metal sites, due to missing linkers, could interact more strongly with adsorbates and result in a slower self-diffusion. Therefore, it is of fundamental importance to evaluate the two competing factors and reveal which one is dominating, a faster self-diffusion due to larger volume or a slower self-diffusion owing to strong interactions at unsaturated sites. In this work, via Monte Carlo and molecular dynamics simulations, we investigate the behavior of isopropyl alcohol (IPA) in the Zr-based UiO-66 MOFs, with a specific focus on the missing linker effects. The results reveal that unsaturated Zr sites bind strongly with IPA molecules, which in return would significantly reduce the self-diffusion coefficient of IPA. Besides this, for the same level of missing linkers, the location of defective sites also makes a difference. We expect such a theoretical study will provide an in-depth understanding of self-diffusion under confinement, inspire better defect engineering strategics, and promote MOF based materials toward challenging real-life applications.

Probing the Influence of Defects, Hydration, and Composition on Prussian Blue Analogues with Pressure

The vast compositional space of Prussian blue analogues (PBAs), formula A_xM[M′(CN)₆]_y·nH₂O, allows for a diverse range of functionality. Yet, the interplay between composition and physical properties—e.g., flexibility and propensity for phase transitions—is still largely unknown, despite its fundamental and industrial relevance. Here we use variable-pressure X-ray and neutron diffraction to explore how key structural features, i.e., defects, hydration, and composition, influence the compressibility and phase behavior of PBAs. Defects enhance the flexibility, manifesting as a remarkably low bulk modulus (B₀ ≈ 6 GPa) for defective PBAs. Interstitial water increases B₀ and enables a pressure-induced phase transition in defective systems. Conversely, hydration does not alter the compressibility of stoichiometric MnPt(CN)₆, but changes the high-pressure phase transitions, suggesting an interplay between low-energy distortions. AMnCo(CN)₆ (A^I = Rb, Cs) transition from F4̅3m to P4̅n2 upon compression due to octahedral tilting, and the critical pressure can be tuned by the A-site cation. At 1 GPa, the symmetry of Rb_0.87Mn[Co(CN)₆]_0.91 is further lowered to the polar space group Pn by an improper ferroelectric mechanism. These fundamental insights aim to facilitate the rational design of PBAs for applications within a wide range of fields.

Water adsorption fingerprinting of structural defects/capping functions in Zr-fumarate MOFs: a hybrid computational-experimental approach

Engineering structural defects in MOFs has been intensively applied to modulate their adsorption-related properties. Zr-fumarate MOF (also known as MOF-801) is a prototypical defective MOF with proven versatile adsorption/separation performances depending on the synthetic conditions, however the relationship between the nature/concentration of both structure defects/capping functions and its adsorption features is still far from being fully understood. In this work, we first present a systematic theoretical exploration of the individual contributions of linker and cluster defects as well as of the capping functions to the overall water adsorption profile of the MOF-801 framework. This computational effort based on the construction of defective structure models and the use of Grand Canonical Monte Carlo simulations further enabled the identification of the overarching defective structure for two MOF-801 samples based on their experimental adsorption isotherms reported previously. An experimental effort was then deployed to synthesize two Zr-fumarate MOF samples with controlled nature and concentration of structural defects as well as capping functions. This computational-experimental hybrid strategy revealed the water adsorption isotherm as a fingerprint of the nature and concentration of structural defect/capping groups exhibited by the MOF adsorbent. We expect this study to deliver meaningful insights to further design MOFs with target adsorption features through a rational engineering of structural defects.

The micromechanical model to computationally investigate cooperative and correlated phenomena in metal-organic frameworks

Computational insight into the impact of cooperative phenomena and correlated spatial disorder on the macroscopic behaviour of metal-organic frameworks (MOFs) is essential in order to consciously engineer these phenomena for targeted applications. However, the spatial extent of these effects, ranging over hundreds of nanometres, limits the applicability of current state-of-the-art computational tools in this field. To obtain a fundamental understanding of these long-range effects, the micromechanical model is introduced here. This model overcomes the challenges associated with conventional coarse-graining techniques by exploiting the natural partitioning of a MOF material into unit cells. By adopting the elastic deformation energy as the central quantity, the micromechanical model hierarchically builds on experimentally accessible input parameters that are obtained from atomistic quantum mechanical or force field simulations. As a result, the here derived micromechanical equations of motion can be adopted to shed light on the effect of long-range cooperative phenomena and correlated spatial disorder on the performance of mesoscale MOF materials.

Defect Engineering of Copper Paddlewheel-Based Metal-Organic Frameworks of Type NOTT-100: Implementing Truncated Linkers and Its Effect on Catalytic Properties

A series of new defect-engineered metal-organic frameworks (DEMOFs) were synthesized by framework doping with truncated linkers employing the mixed-linker approach. Two tritopic defective (truncated) linkers, biphenyl-3,3',5-tricarboxylates (L_H) lacking a ligating group and 5-(5-carboxypyridin-3-yl)isophthalates (L_Py) bearing a weaker interacting ligator site, were integrated into the framework of Cu₂(BPTC) (NOTT-100, BPTC = biphenyl-3,3',5,5'-tetracarboxylates). Incorporating L_H into the framework mainly generates missing metal node defects, thereby obtaining dangling COOH groups in the framework. However, introducing L_Py forms more modified metal nodes featuring reduced and more accessible Cu sites. In comparison with the pristine NOTT-100, the defect-engineered NOTT-100 (DE-NOTT-100) samples show two unique features: (i) functional groups (the protonated carboxylate groups as the Brønsted acid sites or the pyridyl N atoms as the Lewis basic sites), which can act as second active sites, are incorporated into the MOF frameworks, and (ii) more modified paddlewheels, which provided extra coordinatively unsaturated sites, are generated. The cooperative functioning of the above characteristics enhances the catalytic performance of certain types of reactions. For a proof of concept, two exemplary reactions, namely, the cycloaddition of CO₂ with propylene oxide to propylene carbonate and the cyclopropanation of styrene, were carried out to evaluate the catalytic activities of those DE-NOTT-100 materials depending on the defect structure.

Experimental Demonstration of Dynamic Temperature-Dependent Behavior of UiO-66 Metal-Organic Framework: Compaction of Hydroxylated and Dehydroxylated Forms of UiO-66 for High-Pressure Hydrogen Storage

High-pressure (700 MPa or ∼100 000 psi) compaction of dehydroxylated and hydroxylated UiO-66 for H₂ storage applications is reported. The dehydroxylation reaction was found to occur between 150 and 300 °C. The H₂ uptake capacity of powdered hydroxylated UiO-66 reaches 4.6 wt % at 77 K and 100 bar, which is 21% higher than that of dehydroxylated UiO-66 (3.8 wt %). On compaction, the H₂ uptake capacity of dehydroxylated UiO-66 pellets reduces by 66% from 3.8 to 1.3 wt %, while for hydroxylated UiO-66 the pellets show only a 9% reduction in capacity from 4.6 to 4.2 wt %. This implies that the H₂ uptake capacity of compacted hydroxylated UiO-66 is at least three times higher than that of dehydroxylated UiO-66, and therefore, hydroxylated UiO-66 is more promising for hydrogen storage applications. The H₂ uptake capacity is closely related to compaction-induced changes in the porosity of UiO-66. The effect of compaction is greatest in partially dehydroxylated UiO-66 samples that are thermally treated at 200 and 290 °C. These compacted samples exhibit XRD patterns indicative of an amorphous material, low porosity (surface area reduces from between 700 and 1300 m²/g to ca. 200 m²/g and pore volume from between 0.4 and 0.6 cm³/g to 0.1 and 0.15 cm³/g), and very low hydrogen uptake (0.7-0.9 wt % at 77 K and 100 bar). The observed activation-temperature-induced dynamic behavior of UiO-66 is unusual for metal-organic frameworks (MOFs) and has previously only been reported in computational studies. After compaction at 700 MPa, the structural properties and H₂ uptake of hydroxylated UiO-66 remain relatively unchanged but are extremely compromised upon compaction of dehydroxylated UiO-66. Therefore, UiO-66 responds in a dynamic manner to changes in activation temperature within the range in which it has hitherto been considered stable.

Charting the Metal-Dependent High-Pressure Stability of Bimetallic UiO-66 Materials

In theory, bimetallic UiO-66(Zr:Ce) and UiO-66(Zr:Hf) metal-organic frameworks (MOFs) are extremely versatile and attractive nanoporous materials as they combine the high catalytic activity of UiO-66(Ce) or UiO-66(Hf) with the outstanding stability of UiO-66(Zr). Using in situ high-pressure powder X-ray diffraction, however, we observe that this expected mechanical stability is not achieved when incorporating cerium or hafnium in UiO-66(Zr). This observation is akin to the earlier observed reduced thermal stability of UiO-66(Zr:Ce) compounds. To elucidate the atomic origin of this phenomenon, we chart the loss-of-crystallinity pressures of 22 monometallic and bimetallic UiO-66 materials and systematically isolate their intrinsic mechanical stability from their defect-induced weakening. This complementary experimental/computational approach reveals that the intrinsic mechanical stability of these bimetallic MOFs decreases nonlinearly upon cerium incorporation but remains unaffected by the zirconium: hafnium ratio. Additionally, all experimental samples suffer from defect-induced weakening, a synthesis-controlled effect that is observed to be independent of their intrinsic stability.

Engineering a Highly Defective Stable UiO-66 with Tunable Lewis- Brønsted Acidity: The Role of the Hemilabile Linker

The stability of metal-organic frameworks (MOFs) typically decreases with an increasing number of defects, limiting the number of defects that can be created and limiting catalytic and other applications. Herein, we use a hemilabile (Hl) linker to create up to a maximum of six defects per cluster in UiO-66. We synthesized hemilabile UiO-66 (Hl-UiO-66) using benzene dicarboxylate (BDC) as linker and 4-sulfonatobenzoate (PSBA) as the hemilabile linker. The PSBA acts not only as a modulator to create defects but also as a coligand that enhances the stability of the resulting defective framework. Furthermore, upon a postsynthetic treatment in H₂SO₄, the average number of defects increases to the optimum of six missing BDC linkers per cluster (three per formula unit), leaving the Zr-nodes on average sixfold coordinated. Remarkably, the thermal stability of the materials further increases upon this treatment. Periodic density functional theory calculations confirm that the hemilabile ligands strengthen this highly defective structure by several stabilizing interactions. Finally, the catalytic activity of the obtained materials is evaluated in the acid-catalyzed isomerization of α-pinene oxide. This reaction is particularly sensitive to the Brønsted or Lewis acid sites in the catalyst. In comparison to the pristine UiO-66, which mainly possesses Brønsted acid sites, the Hl-UiO-66 and the postsynthetically treated Hl-UiO-66 structures exhibited a higher Lewis acidity and an enhanced activity and selectivity. This is further explored by CD₃CN spectroscopic sorption experiments. We have shown that by tuning the number of defects in UiO-66 using PSBA as the hemilabile linker, one can achieve highly defective and stable MOFs and easily control the Brønsted to Lewis acid ratio in the materials and thus their catalytic activity and selectivity.

Recent Innovation of Metal-Organic Frameworks for Carbon Dioxide Photocatalytic Reduction

The accumulation of carbon dioxide (CO2) pollutants in the atmosphere begets global warming, forcing us to face tangible catastrophes worldwide. Environmental affability, affordability, and efficient CO2 metamorphotic capacity are critical factors for photocatalysts; metal-organic frameworks (MOFs) are one of the best candidates. MOFs, as hybrid organic ligand and inorganic nodal metal with tailorable morphological texture and adaptable electronic structure, are contemporary artificial photocatalysts. The semiconducting nature and porous topology of MOFs, respectively, assists with photogenerated multi-exciton injection and adsorption of substrate proximate to void cavities, thereby converting CO2. The vitality of the employment of MOFs in CO2 photolytic reaction has emerged from the fact that they are not only an inherently eco-friendly weapon for pollutant extermination, but also a potential tool for alleviating foreseeable fuel crises. The excellent synergistic interaction between the central metal and organic linker allows decisive implementation for the design, integration, and application of the catalytic bundle. In this review, we presented recent MOF headway focusing on reports of the last three years, exhaustively categorized based on central metal-type, and novel discussion, from material preparation to photocatalytic, simulated performance recordings of respective as-synthesized materials. The selective CO2 reduction capacities into syngas or formate of standalone or composite MOFs with definite photocatalytic reaction conditions was considered and compared.

Dynamic Interplay between Defective UiO-66 and Protic Solvents in Activated Processes

UiO‐66, composed by Zr‐oxide inorganic bricks [Zr₆(μ₃‐O)₄(μ₃‐OH)₄] and organic terephthalate linkers, is one of the most studied metal–organic frameworks (MOFs) due to its exceptional thermal, chemical, and mechanical stability. Thanks to its high connectivity, the material can withstand structural deformations during activation processes such as linker exchange, dehydration, and defect formation. These processes do alter the zirconium coordination number in a dynamic way, creating open metal sites for catalysis and thus are able to tune the catalytic properties. In this work, it is shown, by means of first‐principle molecular‐dynamics simulations at operating conditions, how protic solvents may facilitate such changes in the metal coordination. Solvent can induce structural rearrangements in the material that can lead to undercoordinated but also overcoordinated metal sites. This is demonstrated by simulating activation processes along well‐chosen collective variables. Such enhanced MD simulations are able to track the intrinsic dynamics of the framework at realistic conditions.

Unraveling the thermodynamic criteria for size-dependent spontaneous phase separation in soft porous crystals

Soft porous crystals (SPCs) harbor a great potential as functional nanoporous materials owing to their stimuli-induced and tuneable morphing between different crystalline phases. These large-amplitude phase transitions are often assumed to occur cooperatively throughout the whole material, which thereby retains its perfect crystalline order. Here, we disprove this paradigm through mesoscale first-principles based molecular dynamics simulations, demonstrating that morphological transitions do induce spatial disorder under the form of interfacial defects and give rise to yet unidentified phase coexistence within a given sample. We hypothesize that this phase coexistence can be stabilized by carefully tuning the experimental control variables through, e.g., temperature or pressure quenching. The observed spatial disorder helps to rationalize yet elusive phenomena in SPCs, such as the impact of crystal downsizing on their flexible nature, thereby identifying the crystal size as a crucial design parameter for stimuli-responsive devices based on SPC nanoparticles and thin films.

Soft porous crystals hold big promise as functional nanoporous materials due to their stimuli responsive flexibility. Here, molecular dynamics simulations reveal a new type of spatial disorder in mesoscale crystals that helps to understand the size-dependency of their phase transition behavior.

Novel metal–organic framework materials: blends, liquids, glasses and crystal–glass composites

This Feature Article reviews a range of amorphisation mechanisms of Metal–Organic Frameworks (MOFs) and presents recent advances to produce novel MOF materials including porous MOF glasses, MOF crystal–glass composites, flux melted MOF glasses and blended zeolitic imidazolate framework glasses.

Defective Metal-Organic Frameworks

The targeted incorporation of defects into crystalline matter allows for the manipulation of many properties and has led to relevant discoveries for optimized and even novel technological applications of materials. It is therefore exciting to see that defects are now recognized to be similarly useful in tailoring properties of metal-organic frameworks (MOFs). For instance, heterogeneous catalysis crucially depends on the number of active catalytic sites as well as on diffusion limitations. By the incorporation of missing linker and missing node defects into MOFs, both parameters can be accessed, improving the catalytic properties. Furthermore, the creation of defects allows for adding properties such as electronic conductivity, which are inherently absent in the parent MOFs. Herein, progress of the rapidly evolving field of the past two years is overviewed, putting a focus on properties that are altered by the incorporation and even tailoring of defects in MOFs. A brief account is also given on the emerging quantitative understanding of defects and heterogeneity in MOFs based on scale-bridging computational modeling and simulations.

Mechanical Properties in Metal-Organic Frameworks: Emerging Opportunities and Challenges for Device Functionality and Technological Applications

Some of the most remarkable recent developments in metal-organic framework (MOF) performance properties can only be rationalized by the mechanical properties endowed by their hybrid inorganic-organic nanoporous structures. While these characteristics create intriguing application prospects, the same attributes also present challenges that will need to be overcome to enable the integration of MOFs with technologies where these promising traits can be exploited. In this review, emerging opportunities and challenges are identified for MOF-enabled device functionality and technological applications that arise from their fascinating mechanical properties. This is discussed not only in the context of their more well-studied gas storage and separation applications, but also for instances where MOFs serve as components of functional nanodevices. Recent advances in understanding MOF mechanical structure-property relationships due to attributes such as defects and interpenetration are highlighted, and open questions related to state-of-the-art computational approaches for quantifying their mechanical properties are critically discussed.

Extension of the QuickFF force field protocol for an improved accuracy of structural, vibrational, mechanical and thermal properties of metal-organic frameworks

QuickFF was originally launched in 2015 to derive accurate force fields for isolated and complex molecular systems in a quick and easy way. Apart from the general applicability, the functionality was especially tested for metal-organic frameworks (MOFs), a class of hybrid materials consisting of organic and inorganic building blocks. Herein, we launch a new release of the QuickFF protocol which includes new major features to predict structural, vibrational, mechanical and thermal properties with greater accuracy, without compromising its robustness and transparent workflow. First, the ab initio data necessary for the fitting procedure may now also be derived from periodic models for the molecular system, as opposed to the earlier cluster-based models. This is essential for an accurate description of MOFs with one-dimensional metal-oxide chains. Second, cross terms that couple internal coordinates (ICs) and anharmonic contributions for bond and bend terms are implemented. These features are essential for a proper description of vibrational and thermal properties. Third, the fitting scheme was modified to improve robustness and accuracy. The new features are tested on MIL-53(Al), MOF-5, CAU-13 and NOTT-300. As expected, periodic input data are proven to be essential for a correct description of structural, vibrational and thermodynamic properties of MIL-53(Al). Bulk moduli and thermal expansion coefficients of MOF-5 are very accurately reproduced by static and dynamic simulations using the newly derived force fields which include cross terms and anharmonic corrections. For the flexible materials CAU-13 and NOTT-300, the transition pressure is accurately predicted provided cross terms are taken into account. © 2018 Wiley Periodicals, Inc.

Dual Role of Water in Heterogeneous Catalytic Hydrolysis of Sarin by Zirconium-Based Metal-Organic Frameworks

Recent experimental studies on Zr^IV-based metal-organic frameworks (MOFs) have shown the extraordinary effectiveness of these porous materials for the detoxification of phosphorus-based chemical warfare agents (CWAs). However, pressing challenges remain with respect to characterizing these catalytic processes both at the molecular and crystalline levels. We here use theory to compare the reactivity of different zirconium-based MOFs for the catalytic hydrolysis of the CWA sarin, using both periodic and cluster modeling. We consider both hydrated and dehydrated secondary building units, as well as linker functionalized MOFs, to more fully understand and rationalize available experimental findings as well as to enable concrete predictions for achieving higher activities for the decomposition of CWAs.

Vacancy defect configurations in the metal-organic framework UiO-66: energetics and electronic structure

Vacancy lattice sites in the metal-organic framework UiO-66 are known to have a profound effect on the material properties. Here we use density functional theory to compare the energies of defect arrangements containing missing linkers and missing metal clusters for different choices of charge compensation. Our results show that the preference for missing metal clusters or missing linker defects depends on the charge compensation as well as the overall concentration of defects in the crystal. Both regimes can be experimentally accessible depending on the synthesis conditions. We investigate the electronic structure of the different types of defects, showing that, despite some changes in the localisation of the frontier orbitals, the electronic energy levels are only weakly affected by the presence of point defects.

On the intrinsic dynamic nature of the rigid UiO-66 metal-organic framework

UiO-66 is a showcase example of an extremely stable metal-organic framework, which maintains its structural integrity during activation processes such as linker exchange and dehydration. The framework can even accommodate a substantial number of defects without compromising its stability. These observations point to an intrinsic dynamic flexibility of the framework, related to changes in the coordination number of the zirconium atoms. Herein we follow the dynamics of the framework in situ, by means of enhanced sampling molecular dynamics simulations such as umbrella sampling, during an activation process, where the coordination number of the bridging hydroxyl groups capped in the inorganic Zr6(μ3-O)4(μ3-OH)4 brick is reduced from three to one. Such a reduction in the coordination number occurs during the dehydration process and in other processes where defects are formed. We observe a remarkable fast response of the system upon structural changes of the hydroxyl group. Internal deformation modes are detected, which point to linker decoordination and recoordination. Detached linkers may be stabilized by hydrogen bonds with hydroxyl groups of the inorganic brick, which gives evidence for an intrinsic dynamic acidity even in the absence of protic guest molecules. Our observations yield a major step forward in the understanding on the molecular level of activation processes realized experimentally but that is hard to track on a purely experimental basis.

Influence of a Confined Methanol Solvent on the Reactivity of Active Sites in UiO-66

UiO-66, composed of Zr-oxide bricks and terephthalate linkers, is currently one of the most studied metal-organic frameworks due to its exceptional stability. Defects can be introduced in the structure, creating undercoordinated Zr atoms which are Lewis acid sites. Here, additional Brønsted sites can be generated by coordinated protic species from the solvent. In this Article, a multilevel modeling approach was applied to unravel the effect of a confined methanol solvent on the active sites in UiO-66. First, active sites were explored with static periodic density functional theory calculations to investigate adsorption of water and methanol. Solvent was then introduced in the pores with grand canonical Monte Carlo simulations, followed by a series of molecular dynamics simulations at operating conditions. A hydrogen-bonded network of methanol molecules is formed, allowing the protons to shuttle between solvent methanol, adsorbed water, and the inorganic brick. Upon deprotonation of an active site, the methanol solvent aids the transfer of protons and stabilizes charged configurations via hydrogen bonding, which could be crucial in stabilizing reactive intermediates. The multilevel modeling approach adopted here sheds light on the important role of a confined solvent on the active sites in the UiO-66 material, introducing dynamic acidity in the system at finite temperatures by which protons may be easily shuttled from various positions at the active sites.

Reliably Modeling the Mechanical Stability of Rigid and Flexible Metal-Organic Frameworks

Over the past two decades, metal-organic frameworks (MOFs) have matured from interesting academic peculiarities toward a continuously expanding class of hybrid, nanoporous materials tuned for targeted technological applications such as gas storage and heterogeneous catalysis. These oft-times crystalline materials, composed of inorganic moieties interconnected by organic ligands, can be endowed with desired structural and chemical features by judiciously functionalizing or substituting these building blocks. As a result of this reticular synthesis, MOF research is situated at the intriguing intersection between chemistry and physics, and the building block approach could pave the way toward the construction of an almost infinite number of possible crystalline structures, provided that they exhibit stability under the desired operational conditions. However, this enormous potential is largely untapped to date, as MOFs have not yet found a major breakthrough in technological applications. One of the remaining challenges for this scale-up is the densification of MOF powders, which is generally achieved by subjecting the material to a pressurization step. However, application of an external pressure may substantially alter the chemical and physical properties of the material. A reliable theoretical guidance that can presynthetically identify the most stable materials could help overcome this technological challenge. In this Account, we describe the recent research the progress on computational characterization of the mechanical stability of MOFs. So far, three complementary approaches have been proposed, focusing on different aspects of mechanical stability: (i) the Born stability criteria, (ii) the anisotropy in mechanical moduli such as the Young and shear moduli, and (iii) the pressure-versus-volume equations of state. As these three methods are grounded in distinct computational approaches, it is expected that their accuracy and efficiency will vary. To date, however, it is unclear which set of properties are suited and reliable for a given application, as a comprehensive comparison for a broad variety of MOFs is absent, impeding the widespread use of these theoretical frameworks. Herein, we fill this gap by critically assessing the performance of the three computational models on a broad set of MOFs that are representative for current applications. These materials encompass the mechanically rigid UiO-66(Zr) and MOF-5(Zn) as well as the flexible MIL-47(V) and MIL-53(Al), which undergo pressure-induced phase transitions. It is observed that the Born stability criteria and pressure-versus-volume equations of state give complementary insight into the macroscopic and microscopic origins of instability, respectively. However, interpretation of the Born stability criteria becomes increasingly difficult when less symmetric materials are considered. Moreover, pressure fluctuations during the simulations hamper their accuracy for flexible materials. In contrast, the pressure-versus-volume equations of state are determined in a thermodynamic ensemble specifically targeted to mitigate the effects of these instantaneous fluctuations, yielding more accurate results. The critical Account presented here paves the way toward a solid computational framework for an extensive presynthetic screening of MOFs to select those that are mechanically stable and can be postsynthetically densified before their use in targeted applications.