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 single component adsorption of alkanes in carbon slit pores was studied using configurational-biased grand canonical Monte Carlo simulations. Wide ranges of temperature, pressure, alkane chain length, and slit height were studied to evaluate their effects on adsorption. Adsorption isotherms and density and orientation profiles were calculated. The behavior of long alkanes at high temperatures was found to be similar to short alkanes at lower temperatures. This suggests that the isotherms may be related through the Polanyi potential theory. [ABSTRACT FROM AUTHOR]

The boundary conditions of intraparticle diffusion in nanoporous materials may be chosen to approach the limiting cases of either absorbing or reflecting boundaries, depending on the host–guest system under study and the temperature of measurement. Pulsed field gradient nuclear magnetic resonance is applied to monitor molecular diffusion of n-hexane and of an n-hexane–tetrafluoromethane mixture adsorbed in zeolite crystallites of type NaX under either of these limiting conditions. Taking advantage of the thus-established peculiarities of mass transfer at the interface between the zeolite bulk phase and the surrounding atmosphere, three independent routes for probing the crystal size are compared. These techniques are based on (i) the measurement of the effective diffusivity under complete confinement, (ii) the application of the so-called NMR tracer desorption technique, and (iii) an analysis of the time dependence of the effective diffusivity in the short-time limit where, by an appropriate variation of the adsorbate and the measuring conditions, the limiting cases of reflecting and adsorbing boundaries could be considered. All these techniques are found to yield coinciding results, which are in excellent agreement with the crystal sizes determined by microscopy. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Presents a method to simulate diffusion and reaction in a single-file system. Construction of a deterministic model consisting of (M+1)[sup N] variables for the system; Reduction of the order of the system by considering only pairs of adjacent cells; Ability of the model to capture correlations between cells when the dominant mode of transport is through single-site hops.

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]

Abstract: The three-point model and molecular simulation are used to understand the enantioselective adsorption of a homologous series of chiral compounds on a common chiral selector that is incorporated on a solid surface. We demonstrate how the enantioselectivity of homologous chiral compounds may go down, go up, remain the same, or even reverse. Conditions that lead to consistent enantioselectivity for homologous chiral compounds are also proposed. [Copyright &y& Elsevier]

Abstract: Metal-organic frameworks are promising materials for the capture of carbon dioxide. Finding the best metal-organic frameworks from the vast number of possibilities could be greatly accelerated by efficient computational screening techniques. We have previously reported an effective screening protocol for predicting carbon dioxide adsorption performance in metal-organic frameworks that uses grand canonical Monte Carlo simulations of gas adsorption. In the model, molecules interact via van der Waals and electrostatic interactions with each other and the framework. However, the method requires single-point quantum mechanics calculations for the estimation of atomic partial charges. In this study we investigate the feasibility of a modified protocol that bypasses the computationally expensive quantum mechanics calculations by applying instead the charge equilibration method. We compare the results of both protocols directly on fourteen metal-organic frameworks and conclude that the new protocol is sufficiently accurate for screening purposes and is significantly faster. [Copyright &y& Elsevier]

Hybrid quantum mechanics/molecular mechanics (QM/MM) methods were employed to study electronic effects of 5,5′-substituents on the enantioselectivity of (salen)Mn catalysts. These methods were shown to effectively capture electronic effects of substituents if the substituents were treated with a simple electrostatic coupling using point charges. Correlations between calculated properties of (salen)MnO indicative of its reactivity and the modified Hammett constant σpara쇪 applied to study electronic framework effects on enantioselectivity for (salen)Mn heterogenized in a metal−organic framework (MOF). The results suggest that these effects are too small to account for all of the decrease in enantioselectivity observed experimentally and that steric constraints of the framework likely play a role. The reactivity correlations also enabled fast ligand screening for the prediction of highly enantioselective (salen)Mn catalysts that could be used in MOFs. A modified version of the DREIDING force field was also developed to allow QM/MM calculations of (salen)Mn species. [ABSTRACT FROM AUTHOR]

Abstract: Hydrogen adsorption isotherms in 20 metal–organic frameworks (MOFs) were obtained from grand canonical Monte Carlo simulations and used to develop quantitative correlations to quickly estimate the H2 uptake of MOFs at 120bar from the free volume of the adsorbent. These correlations were then verified with a diverse set of 22 MOFs that were not used in developing the correlations. The simulated high-pressure H2 isotherms in five representative MOFs show that extremely high pressures (∼80,000bar) are required to completely fill the MOF pores with H2 at 298K. Nevertheless, much lower pressures (several 100bar) are sufficient to attain the H2 storage targets for MOFs with large free volumes. For IRMOF-16 with a large free volume (4.322cm3/g), only modest pressures (170 and 280bar) are required to achieve the 2010 and 2015 gravimetric targets. However, considerably higher pressures are required for attaining the volumetric targets. In addition, three interesting pressure regions, which are observed in the H2 densities within the MOFs (ρad ) up to very high pressures, are analyzed by radial distribution functions between the adsorbed H2 molecules as well as those between H2 molecules and metal corners. [ABSTRACT FROM AUTHOR]

Abstract: Hydrogen storage and delivery in eight representative metal–organic frameworks (MOFs) was investigated using grand canonical Monte Carlo (GCMC) simulations. The simulations demonstrate that the optimal isosteric heat of adsorption (Q st) for maximum hydrogen delivery using MOFs is approximately 20kJ/mol. The results also suggest that increasing the Q st for MOFs with large surface areas (>4800m2/g) is required to attain current hydrogen storage targets in terms of deliverable capacity. [Copyright &y& Elsevier]

Abstract: In this short communication we present a study of xylene adsorption equilibria in distorted and non-distorted AlPO4-11 (AEL) materials using grand canonical Monte Carlo (GCMC) simulations. A previous energy minimization study showed that the AlPO4-11 crystal distorts upon xylene adsorption, but the effect of this distortion on the adsorption loadings was not known. In this work, we compare predictions from GCMC simulations in the distorted and non-distorted frameworks with experiments. The simulations show that the distorted structure adsorbs smaller amounts of xylenes than the non-distorted structure, and simulations in both structures predict ortho-selectivity, as seen in experiment. The lower adsorption capacity of the distorted structure can be related to changes in the adsorption sites. The results indicate that, besides pore blocking, structural deformations can be an important cause of lower adsorption capacity in molecular sieves. [Copyright &y& Elsevier]

Abstract: A cesium-exchanged mordenite zeolite was tested as a hydrocarbon trap for cold-start automotive applications using a modified temperature programmed desorption (TPD) procedure. In the modified TPD procedure, a step increase in feed concentration was performed simultaneously with the start of the temperature ramp to more closely mimic cold start conditions. The feed stream contained mixtures of toluene, propane, and water. The cesium-exchanged mordenite trapped propane under these dynamic conditions, demonstrated improved propane trapping when toluene was coadsorbed, and retained its trapping capacity after nearly 40h of hydrothermal treatment. [Copyright &y& Elsevier]

Abstract: The adsorption of xylene isomers in AlPO4-11 (AEL network) was investigated using biased grand canonical Monte Carlo (GCMC) simulations. Preferential o-xylene adsorption was predicted by the simulations, in agreement with previously reported experimental data. In AlPO4-11 the selective adsorption behavior comes from the smaller length of the o-xylene molecule along the crystallographic c-axis compared to p-xylene. This is in contrast to AlPO4-5 and AlPO4-8, where the ortho-selectivity is caused by the characteristic face-to-face positioning of o-xylene. Energy minimization studies were also performed in a flexible AlPO4-11 lattice to study the structural changes upon xylene adsorption. The energy minimization study showed that the AlPO4-11 crystal distorts upon p- and o-xylene adsorption. The distortion mechanism is related to the strong interaction between xylene methyl groups and the sieve oxygen atoms in the O3 position in the wide region of the pore. [Copyright &y& Elsevier]

The surface area is one of the most important quantities for characterizing novel porous materials. The BET analysis is the standard method for determining surface areas from nitrogen adsorption isotherms and was originally derived for multilayer gas adsorption onto flat surfaces. Metal-organic frameworks (MOFs) are a relatively new class of crystalline, porous materials that have been shown to exhibit very large BET surface areas. These materials are microporous and possess surfaces that are far from flat. In some MOFs, adsorption occurs through a pore-filling mechanism rather than by layer formation. Thus, it is unclear whether BET surface area numbers reported for these materials are truly meaningful. Given the standard practice of reporting BET surface areas for novel porous materials, a critical test of the BET method is much needed. In this work, grand canonical Monte Carlo simulations were used to predict adsorption isotherms for nitrogen in a series of MOFs. The predicted isotherms were used as pseudoexperimental data to test the applicability of the BET theory for obtaining surface areas of microporous MOFs. BET surface areas calculated from the simulated isotherms agree very well with the accessible surface areas calculated directly from the crystal structures in a geometric fashion. In addition, the surface areas agree well with experimental reports in the literature. These results provide a strong validation that the BET theory can be used to obtain surface areas of MOFs. [ABSTRACT FROM AUTHOR]

The extraframework locations of many metal cations in pentasil zeolites, such as ferrierite, are not known, but new infrared bands in the 800–1000 cm−1 region are experimentally observed upon exchange. To investigate the origin of these bands and to better understand cation siting, density functional theory calculations were performed on ferrierite cluster models for seven divalent metals: magnesium, manganese, iron, cobalt, nickel, copper, and zinc. The calculated properties, such as the metal–zeolite binding energy, follow a U-shaped trend with respect to the number of metal d electrons. The calculated infrared bands are consistent with the occupation of the B and G ferrierite sites by the metal cations. However, a definitive assignment of the published experimental infrared bands to particular extraframework sites was not possible. Despite this limitation, all seven metals were shown to perturb the ferrierite framework similarly. All of the cations coordinate to four framework oxygen atoms and the framework T–O bonds associated with these oxygen atoms are elongated by the metal. This elongation is shown to be proportional to the wavenumber of the infrared bands. [Copyright &y& Elsevier]

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]

Two simulation methods are presented for generating a grand canonical ensemble of molecules with internal degrees of freedom. Molecular dynamics integration is used as a Monte Carlo perturbation move in these methods. The integration move helps to explore the intramolecular degrees of freedom with high acceptance ratios. The proofs of the acceptance criteria for both methods are presented. Tests are described to show that the simulation produces results matching an analytical solution for a simple test system. The method is further tested by predicting adsorption isotherms of n-butane in the zeolite silicalite. An example using these algorithms to generate adsorption isotherms of more complex chiral molecules in chiral zeolitelike materials is also presented. [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 article discusses the use of a combined experimental and modeling approach to screen metal-organic frameworks (MOFs) for carbon dioxide (CO2) capture derived from flue gas in line with the continuous attempt of greenhouse gas mitigation. It highlights the significance of understanding MOFs in the fundamentals of developing MOF-based processes and discusses the methods of molecular modeling strategy for (CO2). The advantage of molecular approach is also mentioned.

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]

Highlights: [•] Calculated hydrogen Maxwell–Stefan diffusivity in Mg-functionalized MOFs. [•] Highlighted relationship between hydrogen heat of adsorption and M–S diffusivity. [•] Showed M–S diffusivity decreases exponentially with increasing heat of adsorption. [•] Lessons from this study can be applied to diffusion in other host/guest systems. [ABSTRACT FROM AUTHOR]

Highlights: [•] MOFs with Mg alkoxide functionalizations were studied. [•] Hydrogen diffusion coefficients were obtained from MD simulations. [•] Results suggest interesting differences in diffusivities. [•] Mass transfer should be rapid in hydrogen storage applications of these materials. [ABSTRACT FROM AUTHOR]

Abstract: Molecular simulations were used to screen a diverse collection of eight homochiral metal–organic frameworks (MOFs) for their ability to separate 19 chiral compounds by enantioselective adsorption. The simulation model was validated by comparison with available experimental data. It was found that high enantioselectivity is strongly correlated with a close match between the size of the pore and the size of the chiral sorbate molecule. However, there is also a possibility of no enantioselectivity even when the size of the pore matches with the size of the chiral sorbate molecule. A four-point model was used to explain this observation, and a solution to promote high enantioselectivity has been proposed. [Copyright &y& Elsevier]

Abstract: Experimental results in the literature show that two metal-organic frameworks (MOFs) containing copper and cobalt nodes are active for hydroperoxide decomposition, which is an important reaction in auto-oxidation processes. Density functional theory (DFT) calculations reported here for these systems suggest that the metal sites in the interior of these MOFs are not the active sites for this type of reaction due to the steric effects of the adjacent linkers. This implies that the experimental catalysis observed may occur on the exterior surface of the MOF crystals. Additional calculations with a copper paddlewheel node show that, despite being able to form complexes with hydroperoxides, the metal sites in copper paddlewheels do not catalyze hydroperoxide decomposition. Preliminary calculations involving undercoordinated metal atoms as a model for metal sites on the MOF exterior crystal surface suggest that these sites could be catalytically active. [Copyright &y& Elsevier]

The potential application of a homochiral metal-organic framework (MOF) made from cadmium corners and BINOL-type linkers for enantioselective separations of chiral hydrocarbons was investigated using grand canonical Monte Carlo simulations. Adsorption of racemic mixtures of (R,S)-1,3-dimethyl-1,2-propadiene, (R,S)-1,2-dimethylcyclobutane and (R,S)-1,2-dimethylcyclopropane was simulated. The results show that an enantiomeric excess (ee) of over 50% can be achieved for 1,3-dimethyl-1,2-propadiene, while more moderate ee values were observed for the cyclic compounds. The adsorption sites, diastereomeric complexes and adsorption energies were analysed to understand the enantioselective adsorption. It is shown that the small zigzag pores of the homochiral MOF account primarily for the enantioselective adsorption, whereas the larger helical pores show almost no enantioselectivity. The results also demonstrate that both intrinsic enantioselectivity and accessibility of the adsorption sites are important for the overall enantioselective separation. [ABSTRACT FROM AUTHOR]

Abstract: Quantum chemical calculations were employed to develop a better understanding of the adsorption properties of BaNaY interacting with molecules relevant to deNO x catalysis. First, various basis sets and levels of theory were tested for barium-containing species and gas-phase reactions, and it was shown that the choice of the basis set for barium is critical. Density functional theory (DFT) with the B3LYP functional and SDD as the basis set was selected based on its combination of relative accuracy and speed. This level of theory was then used to calculate energies, geometries, and frequencies for acetaldehyde, acetic acid, nitromethane, and water adsorbed and in the gas phase. The predicted properties were compared to experiment where available and reasonable agreement was found. To study the effect of the zeolite framework on the adsorption properties, the size of the zeolite cluster was increased from 6T to 36T to 96T using the embedded ONIOM method. Inclusion of van der Waals interactions with increasing cluster size only changed the adsorption enthalpy by a small amount for all of the adsorbates studied. It was found that the interaction between the empty 6s orbital of Ba and the lone pair orbital of the oxygen atom in the adsorbed state, revealed by natural bond orbital analysis, correlates with the adsorption enthalpy and the gas-phase charge of the oxygen atom that interacts with Ba. [Copyright &y& Elsevier]

Knudsen dynamics simulations were performed in cylindrical pores with multiple sections of different diameters. The number and location of hits between particles and the pore walls can be related to the known distributions in the uniform diameter cases. Both sweep-gas and pass-through modes of operation were examined. The contact of the particles traveling inside the pores with the pore walls can be regulated by adjusting the aspect ratio (section length/section diameter) of the pore sections and the ratio of their cross-sectional areas. The results were found to be independent of the transition angle between sections except for very smooth transitions. Analytical expressions were developed to relate the transmission probability in pores of multiple sections to the transmission probabilities of the constituent sections. [Copyright &y& Elsevier]

Grand canonical Monte Carlo and molecular dynamics simulations have been performed for methane, n-alkanes, cyclohexane and benzene in two novel nanoporous metal-organic materials. The first material, bipyridine molecular squares, consists of discrete square molecules with corners formed by rhenium complexes and edges formed by bipyridine links, giving a small cavity within each square. The material is considered in both its crystalline form and an amorphous packing of squares. The second material is IRMOF-1, a periodic, crystalline structure also with metal corners and organic bridging units. Adsorption isotherms and self-diffusion coefficients are reported and provide insight into molecular interactions in these materials. [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]