Search

Recent experiments show that the diffusivities of C8 and C10 n‐alkanes in ZIF‐8 are higher than those of C7 and C9, respectively. We investigated this unusual 'odd‐even' effect by simulating the adsorption of C7‐C14 n‐alkanes in ZIF‐8 using hybrid Monte Carlo molecular simulations. The resultant adsorption isotherms, guest‐host energies, isosteric heats, and chain length distributions are analyzed for trends among the n‐alkanes. ZIF‐8 cages filled with n‐alkanes are characterized using a combination of image processing and data science to quantify the differences in packing that occur with chain length. Results indicate that packing changes drastically from C7 to C8 and from C9 to C10, which is consistent with the diffusion trends previously reported in experiments. [ABSTRACT FROM AUTHOR]

This study introduces recommendations for conducting molecular simulations of vapor adsorption, with an emphasis on enhancing the accuracy, reproducibility, and comparability of results. The first aspect we address is consistency in the implementation of some details of typical molecular models, including tail corrections and cutoff distances, due to their significant influence on generated data. We highlight the importance of explicitly calculating the saturation pressures at relevant temperatures using methods such as Gibbs ensemble Monte Carlo simulations and illustrate some pitfalls in extrapolating saturation pressures using this method. For grand canonical Monte Carlo (GCMC) simulations, the input fugacity is usually calculated using an equation of state, which often requires the critical parameters of the fluid. We show the importance of using critical parameters derived from the simulation with the same model to ensure internal consistency between the simulated explicit adsorbate phase and the implicit bulk phase in GCMC. We show the advantages of presenting isotherms on a relative pressure scale to facilitate easier comparison among models and with experiment. Extending these guidelines to a practical case study, we evaluate the performance of various isoreticular metal–organic frameworks (MOFs) in adsorption cooling applications. This includes examining the advantages of using propane and isobutane as working fluids and identifying MOFs with a superior performance. [ABSTRACT FROM AUTHOR]

We investigated eight representative metal–organic frameworks for methane storage using molecular simulation. Validated force fields were used to calculate the amount adsorbed for pure methane and its mixtures with CO2 and H2O at 5.8 and 65 bar at 298 K within the composition limits specified for natural gas. Within the analyzed concentrations, the MOFs without open metal sites were minimally influenced by the presence of CO2 and H2O. However, for the MOFs with open metal sites, the presence of these species proved to be harmful. We found that concentrations as low as 25 ppm of water can reduce the delivered volume of methane by more than 20%. A detailed analysis of the adsorption mechanisms leading to this poisoning is presented. These results highlight the possible limitations of MOFs with open metal sites for use in natural gas storage. [ABSTRACT FROM AUTHOR]

Hydrogen is an appealing energy storage solution for electric vehicles due to its low environmental impact and faster recharge times compared to batteries. However, there are many engineering challenges involved in safely storing a sufficient amount of hydrogen onboard a vehicle with a reasonable volumetric density. Nanoporous materials such as metal–organic frameworks (MOFs) have the potential to store hydrogen at high density and only moderate pressure. Considerable research has been devoted to finding new MOFs for hydrogen storage in recent years; however, a MOF that provides sufficient hydrogen density and is suitable to commercial applications has not yet been found. Much of this research makes use of molecular modelling to screen thousands of materials in a high-throughput way. Computational screening can be an effective tool for gaining insight into structure-performance relationships as well as finding specific candidates for an application. Recently, some research groups have also used machine learning to analyze data more effectively and accelerate the screening process. In this review, we discuss some recent advances in using molecular modelling and machine learning to find materials for hydrogen storage. We also discuss and compare some popular models for the hydrogen molecule and the accuracy of different equations of state, which are important considerations for accurate molecular simulations. [ABSTRACT FROM AUTHOR]

The effective detoxification of chemical warfare agents, specifically nerve agents, is a pressing issue in the modern world. Due to the high toxicity of these molecules, simulants are often used in experiments as substitutes for the agents. However, there is little reason to believe that the current simulants used in the literature are optimal predictors of nerve agent reactivity. Density functional theory calculations were performed on the alkaline hydrolysis of over 100 organophosphate molecules to identify improved simulants for the G‐series nerve agents soman and sarin, based on low toxicity and similarity to nerve agent hydrolysis energetics and degradation mechanism. This screening highlighted 5 molecules that have nearly identical reaction barriers to the actual agents, while being far less toxic. Quantitative structure–activity relationship (QSAR) models were also derived to determine the most significant molecular descriptors for describing the hydrolysis free energy barriers of these reactions. The optimal QSAR model was subjected to a thorough statistical analysis and validation procedure to confirm its predictive capacity, showing excellent quantitative and ranking accuracy. It was further shown that the model trained on G‐series agents can reliably predict energetics for other organophosphate classes as well, including VX. Through these computational insights, experimentalists may be aided in accurately and safely studying these reactions with less toxic simulants. [ABSTRACT FROM AUTHOR]

Density functional theory was used to study interactions between phosphine and 17 functional groups, revealing that only the functional groups containing metal ions bind phosphine strongly. This led to a detailed examination of the adsorption of the pnictogen hydride gases ammonia, phosphine, and arsine on 33 metal catecholates in order to examine periodic trends in the binding strength for these adsorbates. Phosphine and arsine adsorption is primarily driven by donation of electron density from the adsorbate molecule to the metal atom, while ammonia binding involves both electron donation and a significant Coulomb attraction between the negatively charged N atom and the positive metal atom. This Coulomb effect results in notably different binding behavior for ammonia than for phosphine and arsine, which are quite similar to each other. Generally, for metals on the left side of the periodic table, the Coulomb effect dominates in ammonia adsorption, and metals on the right side of the periodic table have lower positive charge and accept greater amounts of electron density, making the electron-sharing contribution more important to the binding strength. We find that more electronegative metals, particularly Pd, Ir, Pt, and Au, which accept more electron density from the adsorbate, yield the best selectivity for the target molecules over water because water adsorption is mostly due to Coulomb attraction. We also analyze changes in the d orbital occupancy upon ammonia and phosphine adsorption. [ABSTRACT FROM AUTHOR]

To investigate postcombustion capture of CO2 in the presence of water, we developed a pressure/vacuum swing adsorption (P/VSA) cycle model consisting of a system of partial differential algebraic equations incorporating mass and energy balances, the Ergun equation for pressure changes, competitive Langmuir isotherms, and the linear driving force model. Four potential adsorbents, zeolites 13X and 5A and the MOFs HKUST-1 and Ni-MOF-74, are investigated, evaluated, and compared. Using this simulation, a two-stage Skarstrom cycle, coupled with an upstream dehydration unit and a downstream compression unit, is optimized using a nondominant sorting genetic algorithm, NSGA-II, to minimize the overall cost of capturing 90% of CO2 from flue gas at a purity of 90% and compressing it for pipeline transportation at 110 bar. The results show that under dry flue gas conditions, zeolite 13X is the best performing adsorbent with an overall cost of $32.1/ton of CO2. Under humid flue gas conditions, zeolites 13X and 5A performed equally well with overall costs of capturing CO2 of approximately $34.1/ton of CO2. [ABSTRACT FROM AUTHOR]

Adsorption hysteresis is a phenomenon related to phase transitions that can impact applications such as gas storage and separations in porous materials. Computational approaches can greatly facilitate the understanding of phase transitions and phase equilibria in porous materials. In this work, adsorption isotherms for methane, ethane, propane, and n-hexane were calculated from atomistic grand canonical Monte Carlo (GCMC) simulations in a metal-organic framework having both micropores and mesopores to better understand hysteresis and phase equilibria between connected pores of different size and the external bulk fluid. At low temperatures, the calculated isotherms exhibit sharp steps accompanied by hysteresis. As a complementary simulation method, canonical (NVT) ensemble simulations with Widom test particle insertions are demonstrated to provide additional information about these systems. The NVT+Widom simulations provide the full van der Waals loop associated with the sharp steps and hysteresis, including the locations of the spinodal points and points within the metastable and unstable regions that are inaccessible to GCMC simulations. The simulations provide molecular-level insight into pore filling and equilibria between high- and low-density states within individual pores. The effect of framework flexibility on adsorption hysteresis is also investigated for methane in IRMOF-1. Developing computational tools to study the phase behaviour of adsorbate molecules in complex pore architectures can greatly facilitate our understanding of phase transitions and phase equilibria in porous materials. Here, molecular simulations are used to study the hysteresis and phase equilibria of n-alkane adsorbate molecules in a metal–organic framework with both micropores and mesopores, and simulations in the canonical ensemble with Widom insertions are shown to be a powerful complementary approach to grand canonical Monte Carlo simulations. [ABSTRACT FROM AUTHOR]

There is an almost unlimited number of metal-organic frameworks (MOFs). This creates exciting opportunities but also poses a problem: how do we quickly find the best MOFs for a given application? Molecular simulations have advanced sufficiently that many MOF properties - especially structural and gas adsorption properties - can be predicted computationally, and molecular modeling techniques are now used increasingly to guide the synthesis of new MOFs. With increasing computational power and improved simulation algorithms, it has become possible to conduct high-throughput computational screening to identify promising MOF structures and uncover structure-property relations. We review these efforts and discuss future directions in this new field. [ABSTRACT FROM AUTHOR]

Metal–organic frameworks (MOFs) with open metal sites have been widely investigated for the selective adsorption of small molecules via redox mechanisms where charge transfer can take place between the binding site and the adsorbate of interest. Quantum-chemical screening methods based on density functional theory have emerged as a promising route to accelerate the discovery of MOFs with enhanced binding affinities toward various adsorbates. However, the success of this approach is linked to the accuracy of the underlying density functional approximations (DFAs). In this work, we compare commonly used generalized gradient approximation (GGA), GGA+U, and meta-GGA exchange-correlation functionals in modeling redox-dependent binding at open metal sites in MOFs using O2 and N2 as representative small molecules. We find that the self-interaction error inherent to the widely used Perdew, Burke, and Ernzerhof (PBE) GGA predicts metal sites that are artificially redox-active, as evidenced by their strong binding affinities, short metal–adsorbate bond distances, and large degree of charge transfer. The incorporation of metal-specific, empirical Hubbard U corrections based on the transition metal oxide literature systematically reduces the redox activity of the open metal sites, often improving agreement with experiment. Additionally, the binding behavior shifts from strong chemisorption to weaker physisorption as a function of U. The M06-L meta-GGA typically predicts binding energies between those of PBE-D3(BJ) and PBE-D3(BJ)+U when using empirically derived U values from the transition metal oxide literature. Despite the strong sensitivity of the binding affinities toward a given DFA, the GGA, GGA+U, and meta-GGA approaches often yield the same qualitative trends and structure–property relationships. [ABSTRACT FROM AUTHOR]

Although technologically promising, the reduction of carbon dioxide (CO2) to produce carbon monoxide (CO) remains economically challenging owing to the lack of an inexpensive, active, highly selective, and stable catalyst. We show that nanocrystalline cubic molybdenum carbide (a-Mo2C), prepared through a facile and scalable route, offers 100% selectivity for CO2 reduction to CO while maintaining its initial equilibrium conversion at high space velocity after more than 500 hours of exposure to harsh reaction conditions at 600°C. The combination of operando and postreaction characterization of the catalyst revealed that its high activity, selectivity, and stability are attributable to crystallographic phase purity, weak CO-Mo2C interactions, and interstitial oxygen atoms, respectively. Mechanistic studies and density functional theory (DFT) calculations provided evidence that the reaction proceeds through an H2-aided redox mechanism. [ABSTRACT FROM AUTHOR]

Highlights from the article: Through quantitative structure-activity relationship modeling and statistical analysis, the most significant molecular descriptors for describing the energetics of alkaline organophosphate hydrolysis reactions are also determined. DFT, hydrolysis, nerve agents, QSAR, simulants.

Through quantum‐chemical calculations, we investigate a family of metal–organic frameworks (MOFs) containing triazolate linkers, M2X2(BBTA) (M=metal, X=bridging anion, H2BBTA=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. [ABSTRACT FROM AUTHOR]

Through quantum‐chemical calculations, we investigate a family of metal–organic frameworks (MOFs) containing triazolate linkers, M2X2(BBTA) (M=metal, X=bridging anion, H2BBTA=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. [ABSTRACT FROM AUTHOR]

• VO x and CeO 2 surface speciation on Al 2 O 3 depends on the ALD growth sequence. • The presence of CeO 2 does not impact the electronic structure of the V O bond. • VO x and CeO 2 on Al 2 O 3 migrate at elevated temperatures under reducing conditions. • CeO 2 clusters on Al 2 O 3 and ceria surface sites in VO x /CeO 2 favor total oxidation. • Different active sites on CeO 2 show varying reactivity towards alkane and alkene. Alumina-supported VO x -CeO 2 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 VO x and CeO 2 surface sites to the cyclohexane ODH reaction mechanism are investigated through a comparison of the catalytic activity of the supported mixed metal oxide materials, VO x /Al 2 O 3 , VO x /CeO 2 and a bare CeO 2 support. The highest activity is observed for CeO 2 , 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 CeO 2 on Al 2 O 3 favor total oxidation, similar to exposed CeO 2 sites in VO x /CeO 2. [ABSTRACT FROM AUTHOR]

This review spotlights the role of atomic‐level modeling in research on metal‐organic frameworks (MOFs), especially the key methodologies of density functional theory (DFT), Monte Carlo (MC) simulations, and molecular dynamics (MD) simulations. The discussion focuses on how periodic and cluster‐based DFT calculations can provide novel insights into MOF properties, with a focus on predicting structural transformations, understanding thermodynamic properties and catalysis, and providing information or properties that are fed into classical simulations such as force field parameters or partial charges. Classical simulation methods, highlighting force field selection, databases of MOFs for high‐throughput screening, and the synergistic nature of MC and MD simulations, are described. By predicting equilibrium thermodynamic and dynamic properties, these methods offer a wide perspective on MOF behavior and mechanisms. Additionally, the incorporation of machine learning (ML) techniques into quantum and classical simulations is discussed. These methods can enhance accuracy, expedite simulation setup, reduce computational costs, as well as predict key parameters, optimize geometries, and estimate MOF stability. By charting the growth and promise of computational research in the MOF field, the aim is to provide insights and recommendations to facilitate the incorporation of computational modeling more broadly into MOF research. [ABSTRACT FROM AUTHOR]

A machine learning (ML) methodology that uses a histogram of interaction energies has been applied to predict gas adsorption in metal–organic frameworks (MOFs) using results from atomistic grand canonical Monte Carlo (GCMC) simulations as training and test data. In this work, the method is first extended to binary mixtures of spherical species, in particular, Xe and Kr. In addition, it is shown that single-component adsorption of ethane and propane can be predicted in good agreement with GCMC simulation using a histogram of the adsorption energies felt by a methyl probe in conjunction with the random forest ML method. The results for propane can be improved by including a small number of MOF textural properties as descriptors. We also discuss the most significant features, which provides physical insight into the most beneficial adsorption energy sites for a given application. [ABSTRACT FROM AUTHOR]

Metal–organic frameworks (MOFs) are a class of nanoporous materials with highly tunable structures in terms of both chemical composition and topology. Due to their tunable nature, high‐throughput computational screening is a particularly appealing method to reduce the time‐to‐discovery of MOFs with desirable physical and chemical properties. In this work, a fully automated, high‐throughput periodic density functional theory (DFT) workflow for screening promising MOF candidates was developed and benchmarked, with a specific focus on applications in catalysis. As a proof‐of‐concept, we use the high‐throughput workflow to screen MOFs containing open metal sites (OMSs) from the Computation‐Ready, Experimental MOF database for the oxidative C—H bond activation of methane. The results from the screening process suggest that, despite the strong C—H bond strength of methane, the main challenge from a screening standpoint is the identification of MOFs with OMSs that can be readily oxidized at moderate reaction conditions. © 2019 Wiley Periodicals, Inc. In this work, a high‐throughput periodic density functional theory (DFT) workflow is developed for the automated screening of metal–organic frameworks for applications in catalysis. As a proof‐of‐concept, the DFT workflow is then used to screen MOFs with open metal sites from the Computation‐Ready, Experimental MOF database to identify promising candidates for methane activation. [ABSTRACT FROM AUTHOR]

The interactions between uranium and non‐innocent organic species are an essential component of fundamental uranium redox chemistry. However, they have seldom been explored in the context of multidimensional, porous materials. Uranium‐based metal–organic frameworks (MOFs) offer a new angle to study these interactions, as these self‐assembled species stabilize uranium species through immobilization by organic linkers within a crystalline framework, while potentially providing a method for adjusting metal oxidation state through coordination of non‐innocent linkers. We report the synthesis of the MOF NU‐1700, assembled from U4+‐paddlewheel nodes and catecholate‐based linkers. We propose this highly unusual structure, which contains two U4+ ions in a paddlewheel built from four linkers—a first among uranium materials—as a result of extensive characterization via powder X‐ray diffraction (PXRD), sorption, transmission electron microscopy (TEM), and thermogravimetric analysis (TGA), in addition to density functional theory (DFT) calculations. [ABSTRACT FROM AUTHOR]

The interactions between uranium and non‐innocent organic species are an essential component of fundamental uranium redox chemistry. However, they have seldom been explored in the context of multidimensional, porous materials. Uranium‐based metal–organic frameworks (MOFs) offer a new angle to study these interactions, as these self‐assembled species stabilize uranium species through immobilization by organic linkers within a crystalline framework, while potentially providing a method for adjusting metal oxidation state through coordination of non‐innocent linkers. We report the synthesis of the MOF NU‐1700, assembled from U4+‐paddlewheel nodes and catecholate‐based linkers. We propose this highly unusual structure, which contains two U4+ ions in a paddlewheel built from four linkers—a first among uranium materials—as a result of extensive characterization via powder X‐ray diffraction (PXRD), sorption, transmission electron microscopy (TEM), and thermogravimetric analysis (TGA), in addition to density functional theory (DFT) calculations. [ABSTRACT FROM AUTHOR]

Metal-organic frameworks (MOFs) are promising materials for a range of energy and environmental applications. Here we describe in detail a computational algorithm and code to generate MOFs based on edge-transitive topological nets for subsequent evaluation via molecular simulation. This algorithm has been previously used by us to construct and evaluate 13 512 MOFs of 41 different topologies for cryo-adsorbed hydrogen storage. Grand canonical Monte Carlo simulations are used here to evaluate the 13 512 structures for the storage of gaseous fuels such as hydrogen and methane and nondistillative separation of xenon/krypton mixtures at various operating conditions. MOF performance for both gaseous fuel storage and xenon/krypton separation is influenced by topology. Simulation data suggest that gaseous fuel storage performance is topology-dependent due to MOF properties such as void fraction and surface area combining differently in different topologies, whereas xenon/krypton separation performance is topology-dependent due to how topology constrains the pore size distribution. [ABSTRACT FROM AUTHOR]

Competitive coadsorption of water is a major problem in the deployment of adsorption-based CO2 capture. Water molecules may compete for adsorption sites, reducing the capacity of the material, and dehumidification prior to separating CO2 from N2 increases process complexity and cost. The development of adsorbent materials that can selectively adsorb CO2 in the presence of water would be a major step forward in the deployment of CO2 capture materials in practice. In this study, large-scale computational screening was carried out to search for metal-organic frameworks (MOFs) with high selectivity toward CO2 over H2O. Calculating framework charges for thousands of MOFs is a significant challenge, so initial screening used a fast, but approximate, charge calculation method. On the basis of the initial screening, 15 MOFs were selected, and Monte Carlo simulations were carried out to compute the adsorption isotherms for these MOFs using more accurate framework charges calculated by density functional theory. A detailed investigation was performed on the effect of using different methods for calculating partial charges, and it was found that electrostatic interactions contribute the majority of the adsorption energy of H2O in the selected MOFs. [ABSTRACT FROM AUTHOR]

The continuous fractional component Monte Carlo method, which was designed to overcome difficulties with insertions and deletions of molecules, is modified to include configurational bias Monte Carlo methods and is further extended to binary systems. The modified method is shown to correctly predict adsorption of Ar in silicalite, Xe and Kr in HKUST-1, and enantiomers in a homochiral metal–organic framework. The modified method is also found to be approximately an order of magnitude more efficient in inserting and deleting molecules than traditional configurational bias grand canonical Monte Carlo simulations in dense systems. [ABSTRACT FROM AUTHOR]

Diffusion of methane, ethane, propane and n-butane was studied within the micropores of several metal-organic frameworks (MOFs) of varying topologies, including the MOFs PCN-14, NU-125, NU-1100 and DUT-49. Diffusion coefficients of the pure components, as well as methane/ethane, methane/propane and methane/butane binary mixtures, were calculated using molecular dynamics simulations to understand the effect of the longer alkanes on uptake of natural gas in MOFs. The calculated self-diffusion coefficients of all four components are on the order of 10 −8 m 2 /s. The diffusion coefficients of the pure components decrease as a function of chain length in all of the MOFs studied and show different behaviour as a function of loading in different MOFs. The self-diffusivities follow the trend D PCN-14 NU-125≈NU-1100>DUT-49. By comparing the diffusion of pure methane and methane mixtures with the higher alkanes, it is observed that the diffusivity of methane is unaffected by the presence of the higher alkanes in the MOFs considered, indicating that the diffusion path of methane is not blocked by the higher alkanes present in natural gas. [ABSTRACT FROM AUTHOR]

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, e.g., above a threshold pressure or temperature. This anomalous adsorption behavior results from a gating effect, where a structural variation of the adsorbent leads to an abrupt change in the gas admission. This review summarizes the mechanisms of the gating effect, which can be a result of the deformation of the framework (e.g., expansion, contraction, reorientation, and sliding of the unit cells), the vibration of the pore-keeping groups (e.g., rotation, swing, and collapse of organic linkers), and the oscillation of the pore-keeping ions (e.g. cesium, potassium, etc.). These structural variations are induced either by the host–guest interaction or by an external stimulus, such as temperature or light, and account for the gating effect at a threshold value of the stimulus. Emphasis is given to the temperature-regulated gating effect, where the critical admission temperature is dictated by the combined effect of the gate opening and thermodynamic factors and plays a key role in regulating guest admission. Molecular simulations can improve our understanding of the gate opening/closing transitions at the atomic scale and enable the construction of quantitative models to describe the gated adsorption behaviour at the macroscale level. The gating effect in porous materials has been widely applied in highly selective gas separation and offers great potential for gas storage and sensing. [ABSTRACT FROM AUTHOR]

The ubiquity of metal–organic frameworks in recent scientific literature underscores their highly versatile nature. MOFs have been developed for use in a wide array of applications, including: sensors, catalysis, separations, drug delivery, and electrochemical processes. Often overlooked in the discussion of MOF-based materials is the mass transport of guest molecules within the pores and channels. Given the wide distribution of pore sizes, linker functionalization, and crystal sizes, molecular diffusion within MOFs can be highly dependent on the MOF-guest system. In this review, we discuss the major factors that govern the mass transport of molecules through MOFs at both the intracrystalline and intercrystalline scale; provide an overview of the experimental and computational methods used to measure guest diffusivity within MOFs; and highlight the relevance of mass transfer in the applications of MOFs in electrochemical systems, separations, and heterogeneous catalysis. [ABSTRACT FROM AUTHOR]

Nanoporous solids are ubiquitous in chemical, energy, and environmental processes, where controlled transport of molecules through the pores plays a crucial role. They are used as sorbents, chromatographic or membrane materials for separations, and as catalysts and catalyst supports. Defined as materials where confinement effects lead to substantial deviations from bulk diffusion, nanoporous materials include crystalline microporous zeotypes and metal–organic frameworks (MOFs), and a number of semi-crystalline and amorphous mesoporous solids, as well as hierarchically structured materials, containing both nanopores and wider meso- or macropores to facilitate transport over macroscopic distances. The ranges of pore sizes, shapes, and topologies spanned by these materials represent a considerable challenge for predicting molecular diffusivities, but fundamental understanding also provides an opportunity to guide the design of new nanoporous materials to increase the performance of transport limited processes. Remarkable progress in synthesis increasingly allows these designs to be put into practice. Molecular simulation techniques have been used in conjunction with experimental measurements to examine in detail the fundamental diffusion processes within nanoporous solids, to provide insight into the free energy landscape navigated by adsorbates, and to better understand nano-confinement effects. Pore network models, discrete particle models and synthesis-mimicking atomistic models allow to tackle diffusion in mesoporous and hierarchically structured porous materials, where multiscale approaches benefit from ever cheaper parallel computing and higher resolution imaging. Here, we discuss synergistic combinations of simulation and experiment to showcase theoretical progress and computational techniques that have been successful in predicting guest diffusion and providing insights. We also outline where new fundamental developments and experimental techniques are needed to enable more accurate predictions for complex systems. [ABSTRACT FROM AUTHOR]

Simulated adsorption isotherms for water in UiO-66 illustrate that defects in the form of missing linkers make this MOF more hydrophilic. Heats of adsorption and density plots further confirm the effect of defects on adsorption of water in UiO-66 at low loadings. In addition, water and CO2 isotherms indicate that not only the amount of defects but their locations within the material affect the loading of guest molecules. [ABSTRACT FROM AUTHOR]

Polymorphism in metal − organic frameworks (MOFs) means that the same chemical building blocks (nodes and linkers) can be used to construct isomeric MOFs with different topological networks. The choice of topology can substantially impact the pore network of the MOF, changing the sizes and shapes of the pores, which has implications for adsorption and separation applications. In this work, we look at the influence of topology in 38 polymorphic MOFs on the separation of linear and branched C4–C6 alkane isomers, a separation of great importance to the petrochemical industry. We find that the MOF Cu 2 (1,4-benzenedicarboxylate) in nbo topology (nbo-Cu 2 BDC) has particularly high affinity for linear alkanes due to its small pore size, which excludes the branched isomers. Upon studying this MOF in further detail, we find that it can take either of two conformations: a cubic conformation, which is typical of nbo MOFs, and a unique star conformation that contains 1D triangular and hexagonal channels. The determination of which conformation the MOF will adopt depends on steric effects between the nodes and linkers. Image 1 [ABSTRACT FROM AUTHOR]

We present a fast and accurate, semi-analytical method for predicting hydrogen adsorption in nanoporous materials. For any temperature and pressure, the adsorbed amount is calculated as an integral over the energy density of adsorption sites (guest-host interactions) plus an average guest-guest term. The guest-host interaction energy is calculated using a classical force field with hydrogen modelled as a single-site probe. The guest-guest interaction energy is approximated using an average coordination number, which is regressed using Gaussian Process Regression (GPR). Local adsorption at each site is then modelled using a Langmuir isotherm, which when weighted with its probability density gives an accurate description of hydrogen adsorption. The method is tested on 933 metal-organic frameworks (MOFs) from the Computation-Ready Experimental (CoRE) MOF database at 77 K from 10 − 5 to 100 bar, and the results are compared against GCMC predictions. To demonstrate the utility of the method, we calculated hydrogen adsorption isotherms for 12,914 existing MOF structures, at two different temperatures at a speed about 100 times that of GCMC simulations and analyzed the results. We found 13 MOFs with predicted deliverable capacities exceeding the DOE target of 50 g/L for adsorption at 100 bar, 77 K and desorption at 5 bar, 160 K. [ABSTRACT FROM AUTHOR]

Six first-row transition metal cations (Mn 2+ , Fe 2+ , Co 2+ , Ni 2+ , Cu 2+ , Zn 2+ ) were evaluated as catalysts for ethene dimerization to 1-butene. This is an important reaction in the chemistry of C C bond formation and in the conversion of natural gas to higher hydrocarbons. Two related classes of transition metal cation catalysts were investigated: 1) single transition metal cations supported on zirconium oxide nodes of the metal–organic framework NU-1000 and 2) small metal hydroxide clusters with two metal atoms (M 2 ) that could be grown by atomic layer deposition on a support exhibiting isolated hydroxyl groups. Using scaling relations, the free energies of co-adsorbed hydrogen and ethene (i.e., (H/C 2 H 4 )*) and adsorbed ethyl (i.e., C 2 H 5 *) were identified as descriptors for ethene dimerization catalysis. Using degree of rate control analysis, it was determined that the rate controlling steps are either ethene insertion (C C bond forming) or β-hydride elimination (C H bond breaking), depending on the metal. Using degree of catalyst control analysis, it was determined that activity on all the catalysts studied could be improved by tuning the free energy of C 2 H 5 *. [ABSTRACT FROM AUTHOR]

The prospect of using copper oxide nanoclusters grown by atomic layer deposition on a porphyrin support for selective oxidation of methane to methanol was examined by means of density functional theory (DFT) calculations. Ab initio thermodynamic analysis indicates that an active site in the form of Cu(μ-O)Cu can be stabilized by activation in O 2 at 465 K. Furthermore, a moderate methane activation energy barrier (E a  = 54 kJ/mol) is predicted, and the hydrogen abstraction activity of the active site could be attributed to the radical character of the bridging oxygen. Methanol extraction in this system is limited by a thermodynamic barrier to desorption of ΔG = 57 kJ/mol at 473 K; however, desorption can be facilitated by the addition of water in a “stepped conversion” process. Overall, our results indicate similar activity between porphyrin-supported copper oxide nanoclusters and existing Cu-exchanged zeolites and provide a computational proof-of-concept for utilizing functionalized organic linkers in metal-organic frameworks (MOFs) for selective oxidation of methane to methanol. [ABSTRACT FROM AUTHOR]

Abstract: DNA G‐quadruplex structures were recently discovered to provide reliable scaffolding for two‐dimensional organic frameworks due to the strong hydrogen‐bonding ability of guanine. Herein, 2,7‐diaryl pyrene building blocks with high HOMO energies and large optical gaps are incorporated into G‐quadruplex organic frameworks. The adjustable substitution on the aryl groups provides an opportunity to elucidate the framework formation mechanism; molecular non‐planarity is found to be beneficial for restricting interlayer slippage, and the framework crystallinity is highest when intermolecular interaction and non‐planarity strike a fine balance. When guanine‐functionalized pyrenes are co‐crystallized with naphthalene diimide, charge‐transfer (CT) complexes are obtained. The photophysical properties of the pyrene‐only and CT frameworks are characterized by UV/Vis and steady‐state and time‐resolved photoluminescence spectroscopies, and by EPR spectroscopy for the CT complex frameworks. [ABSTRACT FROM AUTHOR]

Abstract: DNA G‐quadruplex structures were recently discovered to provide reliable scaffolding for two‐dimensional organic frameworks due to the strong hydrogen‐bonding ability of guanine. Herein, 2,7‐diaryl pyrene building blocks with high HOMO energies and large optical gaps are incorporated into G‐quadruplex organic frameworks. The adjustable substitution on the aryl groups provides an opportunity to elucidate the framework formation mechanism; molecular non‐planarity is found to be beneficial for restricting interlayer slippage, and the framework crystallinity is highest when intermolecular interaction and non‐planarity strike a fine balance. When guanine‐functionalized pyrenes are co‐crystallized with naphthalene diimide, charge‐transfer (CT) complexes are obtained. The photophysical properties of the pyrene‐only and CT frameworks are characterized by UV/Vis and steady‐state and time‐resolved photoluminescence spectroscopies, and by EPR spectroscopy for the CT complex frameworks. [ABSTRACT FROM AUTHOR]

Molecular diffusion is commonly found to slow down with increasing molecular size. Deviations from this pattern occur in some host materials with pore sizes approaching the diameters of the guest molecules. A variety of theoretical models have been suggested to explain deviations from this pattern, but robust experimental data are scarcely available. Here, we present such data, obtained by monitoring the chain length dependence of the uptake of n-alkanes in the zeolitic imidazolate framework ZIF-4. A monotonic decrease in diffusivity from ethane to n-butane was observed, followed by an increase for n-pentane, and another decrease for n-hexane. This observation was confirmed by uptake measurements with n-butane/n-pentane mixtures, which yield faster uptake of n-pentane. Further evidence is provided by the observation of overshooting effects, i.e., by transient n-pentane concentrations exceeding the (eventually attained) equilibrium value. Accompanying grand canonical Monte Carlo simulations reveal, for the larger n-alkanes, significant differences between the adsorbed and gas phase molecular configurations, indicating strong confinement effects within ZIF-4, which, with increasing chain length, may be expected to give rise to configurational shifts facilitating molecular propagation at particular chain lengths. [ABSTRACT FROM AUTHOR]

Computational high-throughput screening was carried out to assess a large number of experimentally reported metal-organic frameworks (MOFs) and zeolites for their utility in hexane isomer separation. Through the work, we identified many MOFs and zeolites with high selectivity (SL+M > 10) for the group of n-hexane, 2-methylpentane, and 3-methylpentane (linear and monobranched isomers) versus 2,2-dimethylbutane and 2,3-dimethylbutane (dibranched isomers). This group of selective sorbents includes VICDOC (Fe2(BDP)3), a MOF with triangular pores that is known to exhibit high isomer selectivity and capacity. For three of these structures, the adsorption isotherms for a 10-component mixture of hexane and heptane isomers were calculated. Subsequent simulations of column breakthrough curves showed that the DEYVUA MOF exhibits a longer process cycle time than VICDOC MOF or MRE zeolite, which are previously reported, high-performing materials, illustrating the importance of capacity in designing MOFs for practical applications. Among the identified candidates, we synthesized and characterized a MOF in a new copper form with high predicted adsorbent capacity (qL+M > 1.2 mol/L) and moderately high selectivity (SL+M ≈ 10). Finally, we examined the role of pore shape in hexane isomer separations, especially of triangular-shaped pores. We show through the potential energy surface and three-dimensional siting analyses that linear alkanes do not populate the corners of narrow triangular channels and that structures with nontriangular pores can efficiently separate hexane isomers. Detailed thermodynamic analysis illustrates how differences in the free energy of adsorption contribute to shape-selective separation in nanoporous materials. [ABSTRACT FROM AUTHOR]

Owing to the vast diversity of linkers, nodes, and topologies, metal–organic frameworks can be tailored for specific tasks, such as chemical separations or catalysis. Accordingly, these materials have attracted significant interest for capture and/or detoxification of toxic industrial chemicals and chemical warfare agents. In this paper, we review recent experimental and computational work pertaining to the capture of several industrially-relevant toxic chemicals, including NH3, SO2, NO2, H2S, and some volatile organic compounds, with particular emphasis on the challenging issue of designing materials that selectively adsorb these chemicals in the presence of water. We also examine recent research on the capture and catalytic degradation of chemical warfare agents such as sarin and sulfur mustard using metal–organic frameworks. [ABSTRACT FROM AUTHOR]

Efforts to computationally characterize large numbers of nanoporous materials often rely on databases of experimentally resolved crystal structures. The accuracy of experimental crystal structures used in such calculations has a significant impact on the reliability of the results. In this work, we report structures optimized using periodic density functional theory (DFT) for more than 800 experimentally synthesized metal-organic frameworks (MOFs). Many MOFs changed significantly upon structural optimization, particularly materials that were crystallographically resolved in their solvated form. For each MOF, we simulated the adsorption of CH4 and CO2 using grand canonical Monte Carlo both before and after DFT optimization. The DFT optimization has a large impact on simulated gas adsorption in some cases. For example, CO2 loading at 1 bar changed by more than 25% in over 25% of the MOFs we considered. [ABSTRACT FROM AUTHOR]

The methane deliverable capacity of adsorbent materials is a critical performance metric that will determine the viability of using adsorbed natural gas (ANG) technology in vehicular applications. ARPA-E recently set a target deliverable capacity of 315 cc(STP)/cc that a viable adsorbent material should achieve to yield a driving range competitive with incumbent fuels. However, recent computational screening of hundreds of thousands of materials suggests that the target is unattainable. In this work, we aim to determine whether the observed limits in deliverable capacity (~200 cc(STP)/cc) are fundamental limits arising from thermodynamic or material design constraints. Our efforts focus on simulating methane adsorption isotherms in a large number of systems, resulting in a broad exploration of different combinations of spatial distributions and energetics of adsorption sites. All systems were classified into five adsorption scenarios with varying degrees of realism in the manner that adsorption sites are created and endowed with energetics. The scenarios range from methane adsorption on discrete idealized lattice sites to adsorption in metal–organic frameworks with coordinatively unsaturated sites (CUS) provided by metalated catechol groups. Our findings strongly suggest that the ARPA-E target is unattainable, although not due to thermodynamic constraints but due to material design constraints. On the other hand, we also find that the currently observed deliverable capacity limits may be moderately surpassed. For instance, incorporation of CUS in IRMOF-10 is predicted to yield a 217 cc(STP)/cc deliverable capacity. The modified material has a ~0.85 void fraction and a heat of adsorption of ~15 kJ/mol. This suggests that similar, moderate improvements over existing materials could be achieved as long as CUS incorporation still maintains a relatively large void fraction. Nonetheless, we conclude that more significant improvements in deliverable capacity will require changes in the currently proposed operation conditions. [ABSTRACT FROM AUTHOR]

Aiming to elucidate guest-induced structural changes in the coordination polymer CPL-2, grand canonical Monte Carlo (GCMC) simulations were used to predict CO2 loadings in this material, and the results were compared with experimental isotherms. Our calculations suggest that CPL-2 exhibits more pronounced CO2-induced structural changes than previously reported. As the initial evidence, the isotherm simulated in the previously reported CPL-2 structure (experimentally resolved from X-ray diffraction in the “as-synthesized” CPL-2) underestimated the measured CO2 loadings at high pressure, indicating that CPL-2 might undergo structural changes that enable higher pore volumes at high pressure. GCMC simulations in CPL-2 structures considering moderate unit cell expansions reported in the literature still underestimated high-pressure experimental loadings. However, considering an incremental rotation of the CPL-2 bipyridyl pillars with increasing CO2 pressure, we were able to trace the measured isotherm with the simulation data. Computational analysis shows that ligand rotation in CPL-2 enables higher pore volumes, which, in turn, accommodate more CO2 as the gas pressure increases. Desorption measurements suggest that hysteresis in the CO2 isotherm of CPL-2 may also be linked to ligand rotation, and the measured adsorption/desorption cycles show that the rotation is reversible. Based on our simulations for CPL-4 and CPL-5 and previously reported experimental data, it is likely that these materials, which differ from CPL-2 in the bipyridyl ligand, behave similarly in the presence of CO2. Our results help understand the behavior of these materials, which present the kind of structural changes that could be potentially exploited to enhance the CO2 working capacity of ultra-microporous materials for carbon capture applications. [ABSTRACT FROM AUTHOR]

Catalytic activity for acyl transfer from Nacetylimidazole (NAI) to different pyridylcarbinol (PC) regioisomers (2-PC, 3-PC, and 4-PC) is demonstrated for a set of topologically diverse, zirconium-based (porphinato)zinc metal-organic frameworks (MOFs). The MOFs studied are PCN-222, MOF-525, and NU-902, which are based on the esq, ftw, and scu topologies, respectively. The experimentally obtained reaction kinetics are discussed in light of molecular modeling results. The catalytic activity is shown to vary across the series of MOFs due to the different extent to which each topology facilitates reactant preconcentration and alignment of PC and NAI via coordination to framework porphyrin sites (orientation effects). Trends of experimental initial reaction rates with MOF topology and PC regioisomer are consistent with preconcentration effects, which depend on the number of porphyrin sites per volume of MOF, as well as with orientation effects, which depend on the number of porphyrin pairs per volume of MOF that bind PC and NAI in such a way that they are primed to form the required transition state. The present work shows how the proper alignment of catalytic linkers can enhance reaction rates in MOFs. [ABSTRACT FROM AUTHOR]

We describe the incorporation of a bistable mechanically interlocked molecule (MIM) into a robust Zr-based metal--organic framework (MOF), NU-1000, by employing a post-synthetic functionalization protocol. On average, close to two bistable [2]catenanes can be incorporated per repeating unit of the hexagonal channels of NU-1000. The reversible redox-switching of the bistable [2]catenanes is retained inside the MOF, as evidenced by solid-state UV-vis-NIR reflectance spectroscopy and cyclic voltammetry. This research demonstrates that bistable MIMs are capable of exhibiting robust dynamics inside the nanopores of a MOF. [ABSTRACT FROM AUTHOR]

Cu-MOF-74 (also known as Cu-CPO-27) was identified as a sorbent having one of the highest densities of Cu(ii) sites per unit volume. Given that Cu(ii) in the framework can be thermally activated to yield a five-coordinate Cu(ii) species, we identified this MOF as a potential candidate for maximal volumetric uptake of ammonia. To that end, the kinetic breakthrough of ammonia in Cu-MOF-74/Cu-CPO-27 was examined under both dry and humid conditions. Under dry conditions the MOF exhibited a respectable performance (2.6 vs. 2.9 NH3 per nm3 for the current record holder HKUST-1), and under 80% relative humidity, the MOF outperformed HKUST-1 (5.9 vs. 3.9 NH3 per nm3, respectively). [ABSTRACT FROM AUTHOR]

Porous metal-organic frameworks (MOFs) have been studied in the context of a wide variety of applications, particularly in relation to molecular storage and separation sciences. Recently, we reported a green, renewable framework material composed ofγ-cyclodextrin (γ-CD) and alkali metal salts--namely, CD-MOF. This porous material has been shown to facilitate the separation of mixtures of alkylaromatic compounds, including the BTEX mixture (benzene, toluene, ethylbenzene, and the regioisomers of xylene), into their pure components, in both the liquid and gas phases, in an energy-efficient manner which could have implications for the petrochemical industry. Here, we report the ability of CD-MOF to separate a wide variety of mixtures, including ethylbenzene from styrene, haloaromatics, terpinenes, pinenes and other chiral compounds. CD-MOF retains saturated compounds to a greater extent than their unsaturated analogues. Also, the location of a double bond within a molecule influences its retention within the extended framework, as revealed in the case of the structural isomers of pinene and terpinine, where the isomers with exocyclic double bonds are more highly retained than those with endocyclic double bonds. The ability of CD-MOF to separate various mono- and disubstituted haloaromatic compounds appears to be controlled by both the size of the halogen substituents and the strength of the noncovalent bonding interactions between the analyte and the framework, an observation which has been confirmed by molecular simulations. Since CD-MOF is a homochiral framework, it is also able to resolve the enantiomers of chiral analytes, including those of limonene and 1-phenylethanol. These findings could lead to cheaper and easier-to-prepare stationary phases for HPLC separations when compared with other chiral stationary phases, such as CD-bonded silica particles. [ABSTRACT FROM AUTHOR]

A new software package, RASPA, for simulating adsorption and diffusion of molecules in flexible nanoporous materials is presented. The code implements the latest state-of-the-art algorithms for molecular dynamics and Monte Carlo (MC) in various ensembles including symplectic/measure-preserving integrators, Ewald summation, configurational-bias MC, continuous fractional component MC, reactive MC and Baker's minimisation. We show example applications of RASPA in computing coexistence properties, adsorption isotherms for single and multiple components, self- and collective diffusivities, reaction systems and visualisation. The software is released under the GNU General Public License. [ABSTRACT FROM PUBLISHER]

Metal-organic frameworks (MOFs) can exhibit exceptionally high surface areas, which are experimentally estimated by applying the BET theory to measured nitrogen isotherms. The Brunauer, Emmett, and Teller (BET)-estimated nitrogen monolayer loading is thus converted to a "BET area," but the meaning of MOF BET areas remains under debate. Recent emphasis has been placed on the usage of four so-called "BET consistency criteria." Using these criteria and simulated nitrogen isotherms for perfect crystals, we calculated BET areas for graphene and 25 MOFs having different pore-size distributions. BET areas were compared with their corresponding geometrically calculated, nitrogen-accessible surface areas (NASAs). Analysis of simulation snapshots elucidated the contributions of "pore-filling" and "monolayer-formation" to the nitrogen adsorption loadings in different MOF pores, revealing the origin of inaccuracies in BET-calculated monolayer loadings, which largely explain discrepancies between BET areas and NASAs. We also find that even if all consistency criteria are satisfied, the BET calculation can significantly overestimate the true monolayer loading, especially in MOFs combining mesopores (d ≥ 20 Å) and large micropores (d = 10-20 Å), due to the overlap of pore-filling and monolayer-formation regimes of these two kinds of pores. While it is not always possible to satisfy all consistency criteria, it is critical to minimize the deviation from these criteria during BET range selection to consistently compare BET areas of different MOFs and for comparing simulated and experimental BET areas of a given MOF. To accurately assess the quality of a MOF sample, it is best to compare experimental BET areas with simulated BET areas rather than with calculated NASAs. [ABSTRACT FROM AUTHOR]

"Breathing" metal-organic frameworks (MOFs) are an emerging class of soft porous crystals (SPCs) with potential for high working capacity for gas storage applications. However, most breathing MOFs have low stability and/or low surface area. Here we report a water-stable, high surface area, breathing MOF of ftw topology, NU-1105. While Zr6-oxo clusters as nodes introduce water stability in NU-1105, its high surface area and breathing character stem from its pyrene-based tetracarboxylate (Py-FP) linkers, in which the fluorene units (F) in the FP "arms" play a key role in promoting breathing behavior. During gas sorption studies, the "closed pore" (cp) ↔ "open pore" (op) transition of NU-1105 occurs at a propane pressure of ~3 bar. At 1 bar, NU-1105 is in its cp form and adsorbs less propane than it would in its op form, highlighting improved working capacity. In situ powder X-ray diffraction during propane sorption was used to track the cp ↔ op transition, and molecular modeling was used to elucidate the structure of the op and cp forms of NU-1105. According to TD-DFT calculations, the proposed conformations of the Py-FP linkers in the op and cp forms are consistent with the measured excitation and emission spectra of the op and cp forms of NU-1105. Similar structural transitions are also observed in the porphyrinic MOF NU-1104 depending on the identity of the porphyrin core; we observed breathing behavior if the constituent Por-PTP linker is nonmetalated. [ABSTRACT FROM AUTHOR]

CO2 capture is essential for reducing the carbon footprint of coal-fired power plants. Here we show, both experimentally and computationally, a new design strategy for capturing CO2 in nanoporous adsorbents. The approach involves ‘complementary organic motifs’ (COMs), which have a precise alignment of charge densities that is complementary to the CO2 quadrupole. Two promising COMs were post-synthetically incorporated into a robust metal–organic framework (MOF) material using solvent-assisted ligand incorporation (SALI). We demonstrate that these COM-functionalized MOFs exhibit high capacity and selectivity for CO2 relative to other reported motifs. [ABSTRACT FROM AUTHOR]

A kinetic analysis of gas-phase cyclohexene epoxidation by H 2 O 2 over mesoporous TS-1 was performed. The production of cyclohexene oxide was very stable with high selectivity. Based on the kinetic analysis, the gas-phase mechanism is proposed to be similar to that of the liquid-phase reaction: an Eley–Rideal type mechanism, in which the reaction between a Ti–OOH intermediate and the physisorbed alkene is the rate-determining step. When the partial pressure of water or H 2 O 2 was varied, a compensation effect was observed. Based on the kinetic model, the compensation effect is attributed to variations in the surface coverage of intermediates, specifically the competitive adsorption of water and H 2 O 2 at the Ti active sites. A meaningful activation energy can only be obtained at high surface coverages of H 2 O 2 and was determined to be 40 ± 2 kJ/mol. [ABSTRACT FROM AUTHOR]