Search

Designing nanoporous catalysts to destroy chemical warfare agents (CWAs) and environmental contaminants requires consideration of both intrinsic catalytic activity and the mass transfer of molecules in and out of the pores. Polar adsorbates such as CWAs experience a heterogeneous environment in many metal-organic frameworks (MOFs) due to the arrangement of the metal nodes and organic linkers of the MOF. However, quantitative relationships between the pore architecture and the resulting diffusion properties of polar molecules have not been established. We used molecular dynamics simulations to calculate the diffusion coefficients of the CWA simulant dimethyl methyl phosphonate (DMMP) in a diverse set of 776 MOFs with Zr

Industrial-scale thermal separation processes have contributed greatly to the rise in carbon dioxide emissions. Porous materials, such as metal-organic frameworks (MOFs), can potentially reduce these emissions by achieving nonthermal chemical separations through the physical adsorption of targeted species with high selectivity. Here, we report the synthesis of the channel-based MOFs

The direct and selective conversion of methane to methanol can be considered a holy grail for catalysis research. In this work, we study a metal-organic framework known as MFU-4l, modified by design to include highly reactive iron-oxo species for the catalytic C-H bond activation of methane. We investigate the oxidation of methane and the further potential oxidation of the product methanol using N₂O as an oxygen source and map the potential energy landscape of these reactions using density functional theory calculations. We show that the highest energy barrier encountered during the methane oxidation process is not the C-H bond breaking, but the activation of the iron center by N₂O. Furthermore, the potential energy landscape for the C-H bond activation exhibits a large, high-energy plateau region instead of a sharp transition state, thus differing from the traditional radical rebound mechanism. This insight offers interesting potential routes to enhance the catalytic activity of the catalyst, to hinder unwanted deactivation pathways, and to reduce the activity towards the over-oxidation of the product.

In the past two decades, various microporous materials have been developed as useful adsorbents for gas adsorption for a wide range of industries. Considerable efforts have been made to regulate the pore accessibility in microporous materials for the manipulation of guest molecules' admission and release. It has long been known that some microporous adsorbents suddenly become highly accessible to guest molecules at specific conditions

An iron–triazolate metal–organic framework (MOF) is computationally investigated for the catalytic oxidation of strong C–H bonds. The MOF is predicted to form reactive iron-oxo active sites, and design rules to guide future experiments are discussed.

Although mass transfer of molecules in and out of porous materials such as zeolites and metal–organic frameworks impacts many applications, the fast and reproducible measurement of intracrystalline...

Hydrogen is an attractive option for energy storage because it can be produced from renewable sources and produces environmentally benign byproducts. However, the volumetric energy density of molecular hydrogen at ambient conditions is low compared to other storage methods like batteries, so it must be compressed to attain a viable energy density for applications such as transportation. Nanoporous materials have attracted significant interest for gas storage because they can attain high storage density at lower pressure than conventional compression. In this work, we examine how to improve the cryogenic hydrogen storage capacity of a series of porous aromatic frameworks (PAFs) by controlling the pore size and increasing the surface area by adding functional groups. We also explore tradeoffs in gravimetric and volumetric measures of the hydrogen storage capacity and the effects of temperature swings using grand canonical Monte Carlo simulations. We also consider the effects of adding functional groups to the metal–organic framework NU-1000 to improve its hydrogen storage capacity. We find that highly flexible alkane chains do not improve the hydrogen storage capacity in NU-1000 because they do not extend into the pores; however, rigid chains containing alkyne groups do increase the surface area and hydrogen storage capacity. Finally, we demonstrate that the deliverable capacity of hydrogen in NU-1000 can be increased from 40.0 to 45.3 g/L (at storage conditions of 100 bar and 77 K and desorption conditions of 5 bar and 160 K) by adding long, rigid alkyne chains into the pores.

Raw input and output data.

We report uptake capacities and transport diffusivities, D, for each of eight linear alkanes (ranging from C5 to C16) in quartz crystal-supported films of solvent-evacuated ZIF-8. Analyses of the alkane uptake profiles revealed that the transport dynamics are governed by guest diffusion through metal-organic framework (MOF) (ZIF-8) crystallites rather than by rates of entry into films at the MOF/vapor interface. The obtained diffusivities range from just over 10-18 m2/s to just under 10-14 m2/s. Notably, minimum cross-sectional widths for all guests exceed the crystallographically measured width of ZIF-8's largest apertures and imply consistently with previous experimental and computational studies that apertures expand to accommodate guest uptake. On average, each additional carbon decreases the transport diffusivity of an alkane by twofold. Closer examination, however, reveals an odd-even effect such that linear alkanes having even numbers of carbons diffuse more rapidly than alkanes featuring one more or one less carbon atom. Thus, ZIF-8's differentiation of transport diffusivities for pairs of alkanes differing in length by only one carbon atom can be significantly greater than the aforementioned factor of 2. Elucidation of the microscopic basis for the odd-even behavior, however, awaits the outcome of molecular dynamics calculations that are beyond the scope of the present study. For compact, solvothermally prepared films, guest transport is dominated by 1D diffusion from the film/vapor interface and toward the underlying quartz crystal. For much lower density, electrophoretically deposited (EPD) films, crystallites behave nearly independently, and guest transport can be adequately modeled by assuming rapid permeation of macroscopic voids between crystallites, followed by entry and rate-limiting radial diffusion into isolated crystallites. One consequence is that EPD films can be much more rapidly infiltrated by molecular guests than can compact, solvothermally grown films. The combined results have potentially favorable implications for the development of kinetic separation schemes for closely related analytes.

Atomically precise cerium oxo clusters offer a platform to investigate structure-property relationships that are much more complex in the ill-defined bulk material cerium dioxide. We investigated the activity of the MCe

Summary The modular nature of metal–organic frameworks (MOFs) enables synthetic control over their physical and chemical properties, but it can be difficult to know which MOFs would be optimal for a given application. High-throughput computational screening and machine learning are promising routes to efficiently navigate the vast chemical space of MOFs but have rarely been used for the prediction of properties that need to be calculated by quantum mechanical methods. Here, we introduce the Quantum MOF (QMOF) database, a publicly available database of computed quantum-chemical properties for more than 14,000 experimentally synthesized MOFs. Throughout this study, we demonstrate how machine learning models trained on the QMOF database can be used to rapidly discover MOFs with targeted electronic structure properties, using the prediction of theoretically computed band gaps as a representative example. We conclude by highlighting several MOFs predicted to have low band gaps, a challenging task given the electronically insulating nature of most MOFs.

Reticular frameworks are crystalline porous materials that form via the self-assembly of molecular building blocks in different topologies, with many having desirable properties for gas storage, separation, catalysis, biomedical applications and so on. The notable variety of building blocks makes reticular chemistry both promising and challenging for prospective materials design. Here we propose an automated nanoporous materials discovery platform powered by a supramolecular variational autoencoder for the generative design of reticular materials. We demonstrate the automated design process with a class of metal–organic framework (MOF) structures and the goal of separating carbon dioxide from natural gas or flue gas. Our model shows high fidelity in capturing MOF structural features. We show that the autoencoder has a promising optimization capability when jointly trained with multiple top adsorbent candidates identified for superior gas separation. MOFs discovered here are strongly competitive against some of the best-performing MOFs/zeolites ever reported. Reticular frameworks are crystalline porous materials with desirable properties such as gas separation, but their large design space presents a challenge. An automated nanoporous materials discovery platform powered by a supramolecular variational autoencoder can efficiently explore this space.

This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (NanoMOFdeli), ERC-2016-COG 726380, Innovate UK (104384) and EPSRC IAA (IAA/RG85685). N.R. acknowledges the support of the Cambridge International Scholarship and the TrinityHenry Barlow Scholarship (Honorary). O.K.F. and R.Q.S. acknowledge funding from the U.S. Department of Energy (DE-FG02-08ER15967). R.S.F. and D.B. acknowledge funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (SCoTMOF), ERC-2015-StG 677289. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia, LLC., a wholly owned subsidiary of Honeywell International, Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA-0003525. The authors gratefully acknowledge funding from the U.S. Department of Energy, Office of Energy Efficiency and Renewable Energy, Hydrogen and Fuel Cell Technologies Office, through the Hydrogen Storage Materials Advanced Research Consortium (HyMARC). This paper describes objective technical results and analysis. Any subjective views or opinions that might be expressed in the paper do not necessarily represent the views of the U.S. Department of Energy or the United States Government. J.D.E. acknowledges the support of the Alexander von Humboldt Foundation and the Center for Information Services and High Performance Computing (ZIH) at TU Dresden. S.K.G. and S.M. acknowledge SERB (Project No. CRG/2019/000906), India for financial support. K.K. and R.K. acknowledge Active Co. Research Grant for funding. S.K. acknowledges funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (COSMOS), ERC-2017-StG 756489. N.L. and J.G.M acknowledge funding from the European Commission through the H2020-MSCA-RISE-2019 program (ZEOBIOCHEM -872102) and the Spanish MICINN and AEI/FEDER (RTI2018-099504-B-C21). N.L. thanks the University of Alicante for funding (UATALENTO17-05). ICN2 is supported by the Severo Ochoa program from the Spanish MINECO (Grant No. SEV-2017-0706) S.M.J.R. and A.L. wish to thank the Fund for Scientific Research Flanders (FWO), under grant nos. 12T3519N and 11D2220N. L.S. was supported by the EPSRC Cambridge NanoDTC EP/L015978/1. C.T.Y. and T.S.N. acknowledges funds from the National Research Foundation of Korea, NRF-2017M3A7B4042140 and NRF-2017M3A7B4042235. P.F. and H. Y. acknowledge US Department of Energy, Office of Basic Energy Sciences, Materials Sciences and Engineering Division under Award No. DE-SC0010596 (P.F.). R.O. would like to acknowledge funding support during his Ph.D. study from Indonesian Endowment Fund for Education-LPDP with the contract No. 202002220216006. Daniel Siderius: Official contribution of the National Institute of Standards and Technology (NIST), not subject to copyright in the United States of America. Daniel Siderius: Certain commercially available items may be identified in this paper. This identification does not imply recommendation by NIST, nor does it imply that it is the best available for the purposes described. B.V.L, S.T.E and A.M.P acknowledge funding from the European Research Council (ERC) under the European Union's Horizon 2020 Research and Innovation Program (Grant agreement no. 639233, COFLeaf)., Porosity and surface area analysis play a prominent role in modern materials science. At the heart of this sits the Brunauer–Emmett–Teller (BET) theory, which has been a remarkably successful contribution to the field of materials science. The BET method was developed in the 1930s for open surfaces but is now the most widely used metric for the estimation of surface areas of microand mesoporous materials. Despite its widespread use, the calculation of BET surface areas causes a spread in reported areas, resulting in reproducibility problems in both academia and industry. To prove this, for this analysis, 18 already-measured raw adsorption isotherms were provided to sixty-one labs, who were asked to calculate the corresponding BET areas. This roundrobin exercise resulted in a wide range of values. Here, the reproducibility of BET area determination from identical isotherms is demonstrated to be a largely ignored issue, raising critical concerns over the reliability of reported BET areas. To solve this major issue, a new computational approach to accurately and systematically determine the BET area of nanoporous materials is developed. The software, called “BET surface identification” (BETSI), expands on the well-known Rouquerol criteria and makes an unambiguous BET area assignment possible., European Research Council (ERC) ERC-2016-COG 726380 ERC-2015-StG 677289 ERC-2017-StG 756489 639233, UK Research & Innovation (UKRI) Innovate UK 104384 UK Research & Innovation (UKRI), Engineering & Physical Sciences Research Council (EPSRC) IAA/RG85685, Cambridge International Scholarship TrinityHenry Barlow Scholarship, United States Department of Energy (DOE) DE-FG02-08ER15967, National Nuclear Security Administration DE-NA-0003525, United States Department of Energy (DOE), Alexander von Humboldt Foundation, Center for Information Services and High Performance Computing (ZIH) at TU Dresden, Department of Science & Technology (India), Science Engineering Research Board (SERB), India CRG/2019/000906, Active Co. Research Grant, European Commission through the H2020-MSCA-RISE-2019 program ZEOBIOCHEM -872102, Spanish MICINN and AEI/FEDER RTI2018-099504-B-C21, University of Alicante UATALENTO17-05, Spanish Government SEV-2017-0706 FWO 12T3519N 11D2220N, UK Research & Innovation (UKRI), Engineering & Physical Sciences Research Council (EPSRC) EP/L015978/1, National Research Foundation of Korea NRF-2017M3A7B4042140 NRF-2017M3A7B4042235, United States Department of Energy (DOE) DE-SC0010596, Indonesian Endowment Fund for Education-LPDP 202002220216006

Zirconium-based metal–organic frameworks (Zr-MOFs) have been increasingly studied over the past two decades as heterogeneous catalysts due to their synthetic tunability, well-defined nature, and ch...

Effective permeation into, and diffusive mass transport within, solvent-filled metal-organic frameworks (MOFs) is critical in applications such as MOF-based chemical catalysis of condensed-phase reactions. In this work, we studied the entry from solution of a luminescent probe molecule, 1,3,5,7-tetramethyl-4,4-difluoroboradiazaindacene (BODIPY), into the 1D channel-type, zirconium-based MOF NU-1008 and subsequent transport of the probe through the MOF. Measurements were accomplished via in situ confocal fluorescence microscopy of individual crystallites, where the evolution of the fluorescence response from the crystallite was followed as functions of both time and location within the crystallite. From the confocal data, intracrystalline transport of BODIPY is well-described by one-dimensional diffusion along the channel direction. Varying the chemical identity of the solvent revealed an inverse dependence of probe-molecule diffusivity on bulk-solvent viscosity, qualitatively consistent with expectations from the Stokes-Einstein equation for molecular diffusion. At a more quantitative level, however, measured diffusion coefficients are about 100-fold smaller than expected from Stokes-Einstein, pointing to substantial channel-confinement effects. Evaluation of the confocal data also reveals a non-negligible mass transport resistance, i.e., surface barrier, associated with the probe molecule leaving the solution and permeating the exterior surface of the MOF. Permeation by the probe entails displacement of solvent from the MOF channels. The magnitude of the resistance increases with the size of the solvent molecule. This work draws attention to the importance of MOF structure, external-surface barriers, and solvent molecule identity to the overall transport process in MOFs, which should assist in understanding the performance of MOFs in applications such as condensed-phase heterogeneous catalysis.

Recent studies have shown that metal-organic frameworks (MOFs) built from hexanuclear M(IV) oxide cluster nodes are effective catalysts for nerve agent hydrolysis, where the properties of the active sites on the nodes can strongly influence the reaction energetics. The connectivity and metal identity of these M6 nodes can be easily tuned, offering extensive opportunities for computational screening to predict promising new materials. Thus, we used density functional theory (DFT) to examine the effects of node topology, connectivity, and metal identity on the binding energies of multiple nerve agents and their corresponding hydrolysis products. By computing an optimization metric based on the relative binding strengths of key hydrolysis reaction species (water, agent, and bidentate-bound products), we predicted optimal M6 nodes for hydrolyzing specific nerve agent and simulant molecules, where our results are in qualitative agreement with observed experimental trends. This analysis highlighted the notion that no single metal or node topology is optimal for all possible organophosphates, suggesting that MOFs should be selected based on the agent of interest. Using the large amount of data generated from our DFT calculations, we then derived quantitative structure-activity relationship (QSAR) models to help explain the complex trends observed in the binding energies. Through linear regression, we identified the most important descriptors for describing the binding of nerve agents and their hydrolysis products to M6 nodes. These results suggested that both molecular and node properties, including both structural and chemical features, collectively contribute to the binding energetics. By performing a thorough statistical analysis, we showed that our QSAR models are capable of making quantitatively accurate binding energy predictions for nerve agents and their hydrolysis products in a wide variety of M(IV)-MOFs. The insights gained herein can be used to guide future experiments for the synthesis of MOFs with enhanced catalytic activity for organophosphate hydrolysis.

Density functional theory is used to study the tunability of trigonal prismatic SBUs found in metal–organic frameworks (MOFs) such as MIL-100, MIL-101, and PCN-250/MIL-127 of chemical composition M...

Through quantum-chemical calculations, we investigate a family of metal-organic frameworks (MOFs) containing triazolate linkers, M2 X2 (BBTA) (M=metal, X=bridging anion, H2 BBTA=1H,5H-benzo(1,2-d:4,5-d')bistriazole), for their ability to form terminal metal-oxo sites and subsequently activate the C-H bond of methane. By varying the metal and bridging anion in the framework, we show how to significantly tune the reactivity of this series of MOFs. The electronic structure of the metal-oxo active site is analyzed for each combination of metal and bridging ligand, and we find that spin density localized on the oxo ligand is not an inherent requirement for low C-H activation barriers. For the Mn- and Fe-containing frameworks, a transition from ferromagnetic to antiferromagnetic coupling between the metal binding site and terminal oxo ligand during the C-H activation process can greatly reduce the kinetic barrier, a unique case of two-state reactivity without a change in the net spin multiplicity.

Alumina-supported VOx-CeO2 materials are synthesized by atomic layer deposition (ALD) and evaluated in cyclohexane oxidative dehydrogenation (ODH). The order of deposition of the two metal oxides by ALD influences the size, geometry and oxidation state of the resulting surface species as indicated by Raman spectroscopy coupled with density functional theory (DFT) calculations and XPS. The individual contributions of VOx and CeO2 surface sites to the cyclohexane ODH reaction mechanism are investigated through a comparison of the catalytic activity of the supported mixed metal oxide materials, VOx/Al2O3, VOx/CeO2 and a bare CeO2 support. The highest activity is observed for CeO2, and different surface active sites are distinguished by introducing cyclohexane and cyclohexene in the reactant feed mixture. The catalytic behavior is dependent on the nature of V O S (S, Support) bonds, and small clusters of CeO2 on Al2O3 favor total oxidation, similar to exposed CeO2 sites in VOx/CeO2.

Metal-organic frameworks (MOFs) with coordinatively unsaturated metal sites are appealing as adsorbent materials due to their tunable functionality and ability to selectively bind small molecules. Through the use of computational screening methods based on periodic density functional theory, we investigate O2 and N2 adsorption at the coordinatively unsaturated metal sites of several MOF families. A variety of design handles are identified that can be used to modify the redox activity of the metal centers, including changing the functionalization of the linkers (replacing oxido donors with sulfido donors), anion exchange of bridging ligands (considering μ-Br-, μ-Cl-, μ-F-, μ-SH-, or μ-OH- groups), and altering the formal oxidation state of the metal. As a result, we show that it is possible to tune the O2 affinity at the open metal sites of MOFs for applications involving the strong and/or selective binding of O2. In contrast with O2 adsorption, N2 adsorption at open metal sites is predicted to be relatively weak across the MOF dataset, with the exception of MOFs containing synthetically elusive V2+ open metal sites. As one example from the screening study, we predicted that exchanging the μ-Cl- ligands of M2Cl2(BBTA) (H2BBTA = 1H,5H-benzo(1,2-d:4,5-d')bistriazole) with μ-OH- groups would significantly enhance the strength of O2 adsorption at the open metal sites without a corresponding increase in the N2 affinity. Experimental investigation of Co2Cl2(BBTA) and Co2(OH)2(BBTA) confirms that the former exhibits weak physisorption of both N2 and O2, whereas the latter is capable of chemisorbing O2 at room temperature in a highly selective manner. The O2 chemisorption behavior is attributed to the greater electron-donating character of the μ-OH- ligands and the presence of H-bonding interactions between the μ-OH- bridging ligands and the reduced O2 adsorbate.

The efficient preparation of single-crystalline ionic polymers and fundamental understanding of their structure-property relationships at the molecular level remains a challenge in chemistry and materials science. Here, we describe the single-crystal structure of a highly ordered polycationic polymer (polyelectrolyte) and its proton conductivity. The polyelectrolyte single crystals can be prepared on a gram-scale in quantitative yield, by taking advantage of an ultraviolet/sunlight-induced topochemical polymerization, from a tricationic monomer-a self-complementary building block possessing a preorganized conformation. A single-crystal-to-single-crystal photopolymerization was revealed unambiguously by in situ single-crystal X-ray diffraction analysis, which was also employed to follow the progression of molecular structure from the monomer, to a partially polymerized intermediate, and, finally, to the polymer itself. Collinear polymer chains are held together tightly by multiple Coulombic interactions involving counterions to form two-dimensional lamellar sheets (1 nm in height) with sub-nanometer pores (5 Å). The polymer is extremely stable under 254 nm light irradiation and high temperature (above 500 K). The extraordinary mechanical strength and environmental stability-in combination with its impressive proton conductivity (∼3 × 10

Using density functional theory (DFT), we studied the catalytic activity of iron oxide nanoclusters that mimic the structure of the active site in the soluble form of methane monooxygenase (sMMO) for the partial oxidation of methane to methanol. Using N

Zirconium-based metal-organic frameworks (Zr-MOFs) have been reported to be efficient catalysts for the hydrolysis of organophosphate chemical warfare agents (CWAs) in buffered solutions. However, for the gas-phase reaction, which is more relevant to the situation in a battlefield gas mask application, the kinetics of Zr-MOF catalysts may be severely hindered by strong product inhibition. To improve the catalytic performance, we computationally screened a series of synthetically accessible Zr-MOF-supported bimetallic metal-oxo clusters in which the metal-oxygen-metal active motif is preserved, aiming to find catalysts that have lower binding affinities to the hydrolysis product. For the promising catalyst Al

Many metal–organic framework (MOF) materials have been reported in the literature as promising for carbon capture applications based on isotherm data or simple adsorbent metrics. However, adsorption process conditions are often neglected in these evaluations. In this study, we performed process-level simulation and optimization of pressure swing adsorption processes on a set of promising MOFs reported in the literature for post-combustion carbon capture. Zeolite 13X was also included as a benchmark material. We examined the ability of the MOFs to achieve the Department of Energy goals of 90% CO2 purity and 90% CO2 recovery by employing process-level optimization using three different cycle configurations: a modified Skarstrom cycle, a five-step cycle, and a fractionated vacuum swing adsorption cycle. Then, we ranked the MOFs based on their economic performance by looking at the productivity and energy requirements for each cycle. We compared this ranking of the MOFs with the rankings provided by other metrics and found that the adsorbent rankings suggested by simplified metrics may differ significantly from the rankings predicted by detailed process optimization. The economic optimization analysis suggests that the best performing MOFs from those analyzed here are UTSA-16, Cu-TDPAT, Zn-MOF-74, Ti-MIL-91, and SIFSIX-3-Ni. Looking at the connection between process performance and material properties, we found that high CO2 working capacity, small pore size, and large difference between the heat of adsorption of CO2 and N2 promote CO2 capture ability based on this small data set. We synthesized one of the top performing MOFs, SIFSIX-3-Ni, and measured its CO2 and N2 adsorption isotherms. The measured isotherms allowed us to estimate the N2 heat of adsorption for SIFSIX-3-Ni, which was not previously available and was required for the process-level modelling.

The rational design and synthesis of robust metal-organic frameworks (MOFs) based on novel organic building blocks are fundamental aspects of reticular chemistry. Beyond simply fabricating new organic linkers, however, it is important to elucidate structure-property relationships at the molecular level to develop high-performing materials. In this work, we successfully targeted a highly porous and robust cage-type MOF (NU-200) with an