Enabling Computational Design of ZIFs Using ReaxFF

Deciphering the controlling factors for phase transitions in zeolitic imidazolate frameworks

Zeolitic imidazolate frameworks (ZIFs) feature complex phase transitions, including polymorphism, melting, vitrification, and polyamorphism. Experimentally probing their structural evolution during transitions involving amorphous phases is a significant challenge, especially at the medium-range length scale. To overcome this challenge, here we first train a deep learning-based force field to identify the structural characteristics of both crystalline and non-crystalline ZIF phases. This allows us to reproduce the structural evolution trend during the melting of crystals and formation of ZIF glasses at various length scales with an accuracy comparable to that of ab initio molecular dynamics, yet at a much lower computational cost. Based on this approach, we propose a new structural descriptor, namely, the ring orientation index, to capture the propensity for crystallization of ZIF-4 (Zn(Im)2, Im = C3H3N2-) glasses, as well as for the formation of ZIF-zni (Zn(Im)2) out of the high-density amorphous phase. This crystal formation process is a result of the reorientation of imidazole rings by sacrificing the order of the structure around the zinc-centered tetrahedra. The outcomes of this work are useful for studying phase transitions in other metal-organic frameworks (MOFs) and may thus guide the development of MOF glasses.

Sub-Micrometer Phonon Mean Free Paths in Metal-Organic Frameworks Revealed by Machine Learning Molecular Dynamics Simulations

Metal-organic frameworks (MOFs) are a family of materials that have high porosity and structural tunability and hold great potential in various applications, many of which require a proper understanding of the thermal transport properties. Molecular dynamics (MD) simulations play an important role in characterizing the thermal transport properties of various materials. However, due to the complexity of the structures, it is difficult to construct accurate empirical interatomic potentials for reliable MD simulations of MOFs. To this end, we develop a set of accurate yet highly efficient machine-learned potentials for three typical MOFs, including MOF-5, HKUST-1, and ZIF-8, using the neuroevolution potential approach as implemented in the GPUMD package, and perform extensive MD simulations to study thermal transport in the three MOFs. Although the lattice thermal conductivity values of the three MOFs are all predicted to be smaller than 1 W/(m K) at room temperature, the phonon mean free paths (MFPs) are found to reach the sub-micrometer scale in the low-frequency region. As a consequence, the apparent thermal conductivity only converges to the diffusive limit for micrometer single crystals, which means that the thermal conductivity is heavily reduced in nanocrystalline MOFs. The sub-micrometer phonon MFPs are also found to be correlated with a moderate temperature dependence of thermal conductivity between those in typical crystalline and amorphous materials. Both the large phonon MFPs and the moderate temperature dependence of thermal conductivity fundamentally change our understanding of thermal transport in MOFs.

Aqueous Stability of Metal-Organic Frameworks Using ReaxFF-Based Metadynamics Simulations

Aqueous stability is a critical property for the application of metal-organic framework (MOF) materials in humid conditions. The sampling of the free energy surface for a water reaction is challenging due to a lack of a reactive force field. Here, we developed a ReaxFF force field for simulating the reaction of zeolitic imidazole frameworks (ZIFs) with water. We carried out metadynamics simulations based on ReaxFF to study the reaction of water with a few different types of MOFs. We also conducted an experimental water immersion test and characterized the XRD, TG, and gas adsorption properties of the MOFs before and after the immersion test. By considering the energy barrier for a hydrolysis reaction, the simulation results are in good agreement with the experiments. MOFs with open structures and large pores are found to be unstable in metadynamics simulations, where the water molecule can attack or bond with the metallic node relatively easily. In contrast, it is more difficult for water to attack the Zn atom in the ZnN₄ tetrahedral structure of ZIFs. We also found that ZIFs with the -NO₂ functional groups have higher water stability. Discrepancies between the metadynamics simulation and gas adsorption experiments have been accounted for from the phase/crystallinity change of the structure reflected in the X-ray diffraction and thermogravimetry analysis of the MOF samples.

Brittle-to-ductile transition and theoretical strength in a metal-organic framework glass

Metal-organic framework (MOF) glasses, a new type of melt-quenched glass, show great promise to deal with the alleviation of greenhouse effects, energy storage and conversion. However, the mechanical behavior of MOF glasses, which is of critical importance given the need for long-term stability, is not well understood. Using both micro- and nanoscale loadings, we find that pillars of a zeolitic imidazolate framework (ZIF) glass have a compressive strength falling within the theoretical strength limit of ≥E/10, a value which is thought to be unreachable in amorphous materials. Pillars with a diameter larger than 500 nm exhibited brittle failure with deformation mechanisms including shear bands and nearly vertical cracks, while pillars with a diameter below 500 nm could carry large plastic strains of ≥20% in a ductile manner with enhanced strength. We report this room-temperature brittle-to-ductile transition in ZIF-62 glass for the first time and demonstrate that theoretical strength and large ductility can be simultaneously achieved in ZIF-62 glass at the nanoscale. Large-scale molecular dynamics simulations have identified that microstructural densification and atomistic rearrangement, i.e., breaking and reconnection of inter-atomistic bonds, were responsible for the exceptional ductility. The insights gained from this study provide a way to manufacture ultra-strong and ductile MOF glasses and may facilitate their processing toward real-world applications.

Nucleation of zeolitic imidazolate frameworks: from molecules to nanoparticles

We have studied the clusters involved in the initial stages of nucleation of Zeolitic Imidazolate Frameworks, employing a wide range of computational techniques. In the pre-nucleating solution, the prevalent cluster is the ZnIm₄ cluster (formed by a zinc cation, Zn²⁺, and four imidazolate anions, Im^-), although clusters such as ZnIm₃, Zn₂Im₇, Zn₂Im₇, Zn₃Im₉, Zn₃Im₁₀, or Zn₄Im₁₂ have energies that are not much higher, so they would also be present in solution at appreciable quantities. All these species, except ZnIm₃, have a tetrahedrally coordinated Zn²⁺ cation. Small Zn_xIm_y clusters are less stable than the ZnIm₄ cluster. The first cluster that is found to be more stable than ZnIm₄ is the Zn₄₁Im₈₈ cluster, which is a disordered cluster with glassy structure. Bulk-like clusters do not begin to be more stable than glassy clusters until much larger sizes, since the larger cluster we have studied (Zn₁₄₄Im₂₈₈) is still less stable than the glassy Zn₄₁Im₈₈ cluster, suggesting that Ostwald's rule (the less stable polymorph crystallizes first) could be fulfilled, not for kinetic, but for thermodynamic reasons. Our results suggest that the first clusters formed in the nucleation process would be glassy clusters, which then undergo transformation to any of the various crystal structures possible, depending on the kinetic routes provided by the synthesis conditions. Our study helps elucidate the way in which the various species present in solution interact, leading to nucleation and crystal growth.

Modeling the Effect of Defects and Disorder in Amorphous Metal-Organic Frameworks

Amorphous metal-organic frameworks (aMOFs) are a class of disordered framework materials with a defined local order given by the connectivity between inorganic nodes and organic linkers, but absent long-range order. The rational development of function for aMOFs is hindered by our limited understanding of the underlying structure-property relationships in these systems, a consequence of the absence of long-range order, which makes experimental characterization particularly challenging. Here, we use a versatile modeling approach to generate in silico structural models for an aMOF based on Fe trimers and 1,3,5-benzenetricarboxylate (BTC) linkers, Fe-BTC. We build a phase space for this material that includes nine amorphous phases with different degrees of defects and local order. These models are analyzed through a combination of structural analysis, pore analysis, and pair distribution functions. Therefore, we are able to systematically explore the effects of the variation of each of these features, both in isolation and combined, for a disordered MOF system, something that would not be possible through experiment alone. We find that the degree of local order has a greater impact on structure and properties than the degree of defects. The approach presented here is versatile and allows for the study of different structural features and MOF chemistries, enabling the derivation of design rules for the rational development of aMOFs.

Beyond the Average: Spatial and Temporal Fluctuations in Oxide Glass-Forming Systems

Atomic structure dictates the performance of all materials systems; the characteristic of disordered materials is the significance of spatial and temporal fluctuations on composition-structure-property-performance relationships. Glass has a disordered atomic arrangement, which induces localized distributions in physical properties that are conventionally defined by average values. Quantifying these statistical distributions (including variances, fluctuations, and heterogeneities) is necessary to describe the complexity of glass-forming systems. Only recently have rigorous theories been developed to predict heterogeneities to manipulate and optimize glass properties. This article provides a comprehensive review of experimental, computational, and theoretical approaches to characterize and demonstrate the effects of short-, medium-, and long-range statistical fluctuations on physical properties (e.g., thermodynamic, kinetic, mechanical, and optical) and processes (e.g., relaxation, crystallization, and phase separation), focusing primarily on commercially relevant oxide glasses. Rigorous investigations of fluctuations enable researchers to improve the fundamental understanding of the chemistry and physics governing glass-forming systems and optimize structure-property-performance relationships for next-generation technological applications of glass, including damage-resistant electronic displays, safer pharmaceutical vials to store and transport vaccines, and lower-attenuation fiber optics. We invite the reader to join us in exploring what can be discovered by going beyond the average.

Mixed metal node effect in zeolitic imidazolate frameworks

We synthesized two series of bimetallic (zinc and cobalt) zeolitic imidazolate frameworks (ZIF-62) under different solvothermal conditions. It is found that the structure of the derived ZIF crystals is highly sensitive to synthesis conditions. One series possesses the standard ZIF-62 structure, whereas the other has a mixed structure composed of both the standard structure and an unknown one. The standard series exhibits a slight negative deviation from linearity of melting temperature (T _m) and glass transition temperature (T _g) with the substitution of Co for Zn. In contrast, the new series displays a stronger negative deviation. These negative deviations from linearity indicate the mixed metal node effect in bimetallic ZIF-62 due to the structural mismatch between Co²⁺ and Zn²⁺ and to the difference in their electronic configurations. The new series involves both cobalt-rich and zinc-rich phases, whereas the standard one shows one homogeneous phase. Density functional theory calculations predict that the substitution of Co for Zn increases the bulk modulus of the ZIF crystals. This work indicates that the structure, melting behaviour, and mechanical properties of ZIFs can be tuned by metal node substitution and by varying the synthetic conditions. Both series of ZIFs have higher glass forming abilities due to their higher T _g/T _m ratios (0.77-0.84) compared to most good glass formers.

Towards Predicting the Sequential Appearance of Zeolitic Imidazolate Frameworks Synthesized by Mechanochemistry

We use computational materials methods to study the sequential appearance of zinc-based zeolitic imidazolate frameworks (ZIFs) generated in the mechanochemical conversion process. We consider nine ZIF topologies, namely RHO, ANA, QTZ, SOD, KAT, DIA, NEB, CAG and GIS, combined with the two ligands 2-methylimidazolate and 2-ethylimidazolate. Of the 18 combinations obtained, only six (three for each ligand) were actually observed during the mechanosynthesis process. Energy and porosity calculations based on density functional theory, in combination with the Ostwald rule of stages, were found to be insufficient to distinguish the experimentally observed ZIFs. We then show, using classical molecular dynamics, that only ZIFs withstanding quasi-hydrostatic pressure P ≥ 0.3 GPa without being destroyed were observed in the laboratory. This finding, along with the requirement that successive ZIFs be generated with decreasing porosity and/or energy, provides heuristic rules for predicting the sequences of mechanically generated ZIFs for the two ligands considered.

Capturing Free-Radical Polymerization by Synergetic Ab Initio Calculations and Topological Reactive Molecular Dynamics

Photocurable polymers are used ubiquitously in 3D printing, coatings, adhesives, and composite fillers. In the present work, the free radical polymerization of photocurable compounds is studied using reactive classical molecular dynamics combined with a dynamical approach of the nonequilibrium molecular dynamics (D-NEMD). Different concentrations of radicals and reaction velocities are considered. The mechanical properties of the polymer resulting from 1,6-hexanediol dimethacrylate systems are characterized in terms of viscosity, diffusion constant, and activation energy, whereas the topological ones through the number of cycles (polymer loops) and cyclomatic complexity. Effects like volume shrinkage and delaying of the gel point for increasing monomer concentration are also predicted, as well as the stress-strain curve and Young's modulus. Combining ab initio, reactive molecular dynamics, and the D-NEMD method might lead to a novel and powerful tool to describe photopolymerization processes and to original routes to optimize additive manufacturing methods relying on photosensitive macromolecular systems.

ResponZIF Structures: Zeolitic Imidazolate Frameworks as Stimuli-Responsive Materials

Zeolitic imidazolate frameworks (ZIFs) have long been recognized as a prominent subset of the metal-organic framework (MOF) family, in part because of their ease of synthesis and good thermal and chemical stability, alongside attractive properties for diverse potential applications. Prototypical ZIFs like ZIF-8 have become embodiments of the significant promise held by porous coordination polymers as next-generation designer materials. At the same time, their intriguing property of experiencing significant structural changes upon the application of external stimuli such as temperature, mechanical pressure, guest adsorption, or electromagnetic fields, among others, has placed this family of MOFs squarely under the umbrella of stimuli-responsive materials. In this review, we provide an overview of the current understanding of the triggered structural and electronic responses observed in ZIFs (linker and bond dynamics, crystalline and amorphous phase changes, luminescence, etc.). We then describe the state-of-the-art experimental and computational methodology capable of shedding light on these complex phenomena, followed by a comprehensive summary of the stimuli-responsive nature of four prototypical ZIFs: ZIF-8, ZIF-7, ZIF-4, and ZIF-zni. We further expose the relevant challenges for the characterization and fundamental understanding of responsive ZIFs, including how to take advantage of their flexible properties for new application avenues.

Molecular Dynamics and Machine Learning in Catalysts

Given the importance of catalysts in the chemical industry, they have been extensively investigated by experimental and numerical methods. With the development of computational algorithms and computer hardware, large-scale simulations have enabled influential studies with more atomic details reflecting microscopic mechanisms. This review provides a comprehensive summary of recent developments in molecular dynamics, including ab initio molecular dynamics and reaction force-field molecular dynamics. Recent research on both approaches to catalyst calculations is reviewed, including growth, dehydrogenation, hydrogenation, oxidation reactions, bias, and recombination of carbon materials that can guide catalyst calculations. Machine learning has attracted increasing interest in recent years, and its combination with the field of catalysts has inspired promising development approaches. Its applications in machine learning potential, catalyst design, performance prediction, structure optimization, and classification have been summarized in detail. This review hopes to shed light and perspective on ML approaches in catalysts.

Molecular Interactions and Layer Stacking Dictate Covalent Organic Framework Effective Pore Size

Interactions among ions, molecules, and confining solid surfaces are universally challenging and intriguing topics. Lacking a molecular-level understanding of such interactions in complex organic solvents perpetuates the intractable challenge of simultaneously achieving high permeance and selectivity in selectively permeable barriers. Two-dimensional covalent organic frameworks (COFs) have demonstrated ultrahigh permeance, high selectivity, and stability in organic solvents. Using reactive force field molecular dynamics modeling and direct experimental comparisons of an imine-linked carboxylated COF (C-COF), we demonstrate that unprecedented organic solvent nanofiltration separation performance can be accomplished by the well-aligned, highly crystalline pores. Furthermore, we show that the effective, as opposed to designed, pore size and solvated solute radii can change dramatically with the solvent environment, providing insights into complex molecular interactions and enabling future application-specific material design and synthesis.

Mixed hierarchical local structure in a disordered metal-organic framework

Amorphous metal-organic frameworks (MOFs) are an emerging class of materials. However, their structural characterisation represents a significant challenge. Fe-BTC, and the commercial equivalent Basolite® F300, are MOFs with incredibly diverse catalytic ability, yet their disordered structures remain poorly understood. Here, we use advanced electron microscopy to identify a nanocomposite structure of Fe-BTC where nanocrystalline domains are embedded within an amorphous matrix, whilst synchrotron total scattering measurements reveal the extent of local atomic order within Fe-BTC. We use a polymerisation-based algorithm to generate an atomistic structure for Fe-BTC, the first example of this methodology applied to the amorphous MOF field outside the well-studied zeolitic imidazolate framework family. This demonstrates the applicability of this computational approach towards the modelling of other amorphous MOF systems with potential generality towards all MOF chemistries and connectivities. We find that the structures of Fe-BTC and Basolite® F300 can be represented by models containing a mixture of short- and medium-range order with a greater proportion of medium-range order in Basolite® F300 than in Fe-BTC. We conclude by discussing how our approach may allow for high-throughput computational discovery of functional, amorphous MOFs.

Effect of diol isomer/water mixtures on the stability of Zn-MOF-74

The stability of metal-organic frameworks is a key factor in many applications in some fields that require working under harsh conditions. It is known that a large number of MOFs are vulnerable to humid air. It means that when they are exposed to water, a structural collapse of the crystal happens. In this work, Molecular Dynamics simulations using a reactive force field have been performed to study the stability of MOF-74 against the adsorption of catechol, resorcinol and hydroquinone in the presence of water. We reproduced the water instability of Zn-MOF-74 and we studied the resistance of the structure. Our simulations showed that the three isomers generate a volume change in the framework but the structural collapse does not happen. In contrast, for water-isomer mixtures, there is structural collapse. Not only do catechol, resorcinol and hydroquinone not behave as stabilizing agents but they do enhance the hydration effect on the structure.

Atomistic Mechanisms of Thermal Transformation in a Zr-Metal Organic Framework, MIL-140C

To understand the mechanisms responsible for thermal decomposition of a Zr-MOF (MIL-140C), we perform atomistic-scale molecular dynamics (MD) simulations and discuss the simulation data in comparison with the TEM images obtained for the decomposed Zr-MOF. First, we introduce the ReaxFF parameters suitable for the Zr/C/H/O chemistry and then apply them to investigate the thermal stability and morphological changes in the MIL-140C during heating. Based on the performed simulations we propose an atomic mechanism for the collapse of the MIL-140C and the molecular pathways for carbon monoxide formation, the main product of the MIL-140C thermal degradation. We also determine that the oxidation state of the ZrO_x clusters, evolved due to the thermal degradation, approximates the tetragonal phase of ZrO₂. Both simulations and experiments show a distribution of very small ZrO_x clusters embedded in the disrupted organic sheet that could contribute to the unusual high catalytic activity of the decomposed MIL-140C.

Bond switching is responsible for nanoductility in zeolitic imidazolate framework glasses

The fracture mechanism of zeolitic imidazolate framework (ZIF) glasses is revealed to be associated with bond switching of organic linkers around central Zn nodes. The bond switching is more pronounced for ZIF glasses with smaller organic linkers.

Electronic Structure Modeling of Metal-Organic Frameworks

Owing to their molecular building blocks, yet highly crystalline nature, metal-organic frameworks (MOFs) sit at the interface between molecule and material. Their diverse structures and compositions enable them to be useful materials as catalysts in heterogeneous reactions, electrical conductors in energy storage and transfer applications, chromophores in photoenabled chemical transformations, and beyond. In all cases, density functional theory (DFT) and higher-level methods for electronic structure determination provide valuable quantitative information about the electronic properties that underpin the functions of these frameworks. However, there are only two general modeling approaches in conventional electronic structure software packages: those that treat materials as extended, periodic solids, and those that treat materials as discrete molecules. Each approach has features and benefits; both have been widely employed to understand the emergent chemistry that arises from the formation of the metal-organic interface. This Review canvases these approaches to date, with emphasis placed on the application of electronic structure theory to explore reactivity and electron transfer using periodic, molecular, and embedded models. This includes (i) computational chemistry considerations such as how functional, k-grid, and other model variables are selected to enable insights into MOF properties, (ii) extended solid models that treat MOFs as materials rather than molecules, (iii) the mechanics of cluster extraction and subsequent chemistry enabled by these molecular models, (iv) catalytic studies using both solids and clusters thereof, and (v) embedded, mixed-method approaches, which simulate a fraction of the material using one level of theory and the remainder of the material using another dissimilar theoretical implementation.

Fracture toughness of a metal-organic framework glass

Metal-organic framework glasses feature unique thermal, structural, and chemical properties compared to traditional metallic, organic, and oxide glasses. So far, there is a lack of knowledge of their mechanical properties, especially toughness and strength, owing to the challenge in preparing large bulk glass samples for mechanical testing. However, a recently developed melting method enables fabrication of large bulk glass samples (>25 mm³) from zeolitic imidazolate frameworks. Here, fracture toughness (K_Ic) of a representative glass, namely ZIF-62 glass (Zn(C₃H₃N₂)_1.75(C₇H₅N₂)_0.25), is measured using single-edge precracked beam method and simulated using reactive molecular dynamics. K_Ic is determined to be ~0.1 MPa m^0.5, which is even lower than that of brittle oxide glasses due to the preferential breakage of the weak coordinative bonds (Zn-N). The glass is found to exhibit an anomalous brittle-to-ductile transition behavior, considering its low fracture surface energy despite similar Poisson’s ratio to that of many ductile metallic and organic glasses.

Metal-organic framework glasses are gaining interest, but large samples are difficult to fabricate and mechanical properties are not well understood. Here, the authors use experiments and simulations to assess fracture toughness and flexural strength of a zeolitic imidazolate framework glass.

Metal-Organic Framework Glasses Possess Higher Thermal Conductivity than Their Crystalline Counterparts

The thermal conductivity (κ) of glasses is known to always be lower than that of their corresponding crystals due to the stronger phonon-phonon scattering in the former. However, it is unknown whether this relation holds for metal-organic frameworks. Here, we report our discovery of an inverse relation in κ between glass and crystal for two zeolitic imidazolate frameworks (ZIFs), ZIF-4 and ZIF-62, that is, melt-quenched ZIF-4 and ZIF-62 glasses possess higher thermal conductivities than their crystalline counterparts. We find that the ZIF crystal pellets exhibit ultralow κ (∼0.1 W m^-1 K^-1) and that the higher κ of the ZIF glasses is due to the collapse of internal cavities and higher atomic number density in the latter. For other systems like oxides, vitrification causes higher free volume, but the opposite is found for the ZIFs, that is, lower free volume owing to the partial collapse of the crystalline framework upon melting.

ZIF-8 nanoparticles thin film at an oil–water interface as an electrocatalyst for the methanol oxidation reaction without the application of noble metals

Monodispersed and nano-sized ZIF-8 was synthesized at an oil–water interface for the first time and applied as an electrocatalyst for the methanol oxidation reaction in an alkaline medium.