In silico prediction of MOFs with high deliverable capacity or internal surface area

Computational Simulations of Metal-Organic Frameworks to Enhance Adsorption Applications

Metal-organic frameworks (MOFs), renowned for their exceptional porosity and crystalline structure, stand at the forefront of gas adsorption and separation applications. Shortly after their discovery through experimental synthesis, computational simulations quickly become an important method in broadening the use of MOFs by offering deep insights into their structural, functional, and performance properties. This review specifically addresses the pivotal role of molecular simulations in enlarging the molecular understanding of MOFs and enhancing their applications, particularly for gas adsorption. After reviewing the historical development and implementation of molecular simulation methods in the field of MOFs, high-throughput computational screening (HTCS) studies used to unlock the potential of MOFs in CO2 capture, CH4 storage, H2 storage, and water harvesting are visited and recent advancements in these adsorption applications are highlighted. The transformative impact of integrating artificial intelligence with HTCS on the prediction of MOFs' performance and directing the experimental efforts on promising materials is addressed. An outlook on current opportunities and challenges in the field to accelerate the adsorption applications of MOFs is finally provided.

Coarse-Grained Crystal Graph Neural Networks for Reticular Materials Design

Reticular materials, including metal-organic frameworks and covalent organic frameworks, combine the relative ease of synthesis and an impressive range of applications in various fields from gas storage to biomedicine. Diverse properties arise from the variation of building units─metal centers and organic linkers─in almost infinite chemical space. Such variation substantially complicates the experimental design and promotes the use of computational methods. In particular, the most successful artificial intelligence algorithms for predicting the properties of reticular materials are atomic-level graph neural networks, which optionally incorporate domain knowledge. Nonetheless, the data-driven inverse design involving these models suffers from the incorporation of irrelevant and redundant features such as a full atomistic graph and network topology. In this study, we propose a new way of representing materials, aiming to overcome the limitations of existing methods; the message passing is performed on a coarse-grained crystal graph that comprises molecular building units. To highlight the merits of our approach, we assessed the predictive performance and energy efficiency of neural networks built on different materials representations, including composition-based and crystal-structure-aware models. Coarse-grained crystal graph neural networks showed decent accuracy at low computational costs, making them a valuable alternative to omnipresent atomic-level algorithms. Moreover, the presented models can be successfully integrated into an inverse materials design pipeline as estimators of the objective function. Overall, the coarse-grained crystal graph framework is aimed at challenging the prevailing atom-centric perspective on reticular materials design.

Computational Identification and Experimental Demonstration of High-Performance Methane Sorbents

Remarkable methane uptake is demonstrated experimentally in three metal-organic frameworks (MOFs) identified by computational screening: UTSA-76, UMCM-152 and DUT-23-Cu. These MOFs outperform the benchmark sorbent, HKUST-1, both volumetrically and gravimetrically, under a pressure swing of 80 to 5 bar at 298 K. Although high uptake at elevated pressure is critical for achieving this performance, a low density of high-affinity sites (coordinatively unsaturated metal centers) also contributes to a more complete release of stored gas at low pressure. The identification of these MOFs facilitates the efficient storage of natural gas via adsorption and provides further evidence of the utility of computational screening in identifying overlooked sorbents.

In Silico Evolution of High-Performing Metal Organic Frameworks for Methane Adsorption

The increased use of transition fuels, such as natural gas, and the resulting increase in methane emissions have resulted in a need for novel methane storage materials. Metal-organic frameworks (MOFs) have shown promise as efficient storage materials. A virtually limitless number of potential MOFs can be hypothesized, which exhibit a wide variety of different structural and chemical characteristics. Because of the numerous possibilities, identification of the best MOF for methane storage can be a potentially challenging problem. In this work, determination of the best such MOF was cast as an inverse function problem. The function, a random forest (RF) model using 12 structural and chemical descriptors, was trained on 10% of a data set consisting of 130 398 hypothetical MOFs (hMOFs) to predict simulated methane uptake. The RF model was tested on the remaining 90% of the data. After validation, a genetic algorithm (GA) was used to evolve in silico the best MOFs for methane adsorption. The RF model was imbedded into the GA as the fitness function to predict the methane uptake of the evolved MOFs (eMOFs). The best 15 eMOFs matched hMOFs found in the top 1% of the database. Nine of the 15 eMOFs were found in the top 0.1%. More impressively, two of the eMOFs matched the top two hypothetical MOFs with the highest methane uptake values out of the entire database of 130 398 MOFs. Further, by leveraging the ensemble nature of the GA, it was possible to characterize the importance of the different material properties for methane adsorption, providing fundamental insight for future material design strategies.

Predicting hydrogen storage in MOFs via machine learning

The H2 capacities of a diverse set of 918,734 metal-organic frameworks (MOFs) sourced from 19 databases is predicted via machine learning (ML). Using only 7 structural features as input, ML identifies 8,282 MOFs with the potential to exceed the capacities of state-of-the-art materials. The identified MOFs are predominantly hypothetical compounds having low densities (<0.31 g cm-3) in combination with high surface areas (>5,300 m2 g-1), void fractions (∼0.90), and pore volumes (>3.3 cm3 g-1). The relative importance of the input features are characterized, and dependencies on the ML algorithm and training set size are quantified. The most important features for predicting H2 uptake are pore volume (for gravimetric capacity) and void fraction (for volumetric capacity). The ML models are available on the web, allowing for rapid and accurate predictions of the hydrogen capacities of MOFs from limited structural data; the simplest models require only a single crystallographic feature.

Machine Learning Enabled Tailor-Made Design of Application-Specific Metal-Organic Frameworks

In the development of advanced nanoporous materials, one clear and unavoidable challenge in hand is the sheer size (in principle, infinite) of the materials space to be explored. While high-throughput screening techniques allow us to narrow down the enormous-scale database of nanoporous materials, there are still practical limitations stemming from a costly molecular simulation in estimating a material's performance and the necessity of a sophisticated descriptor identifying materials. With an attempt to transition away from the screening-based approaches, this paper presents a computational approach combining the Monte Carlo tree search and recurrent neural networks for the tailor-made design of metal-organic frameworks toward the desired target applications. In the demonstration cases for methane-storage and carbon-capture applications, our approach showed significant efficiency in designing promising and novel metal-organic frameworks. We expect that this approach would easily be extended to other applications by simply adjusting the reward function according to the target performance property.

Inverse design of porous materials using artificial neural networks

Generating optimal nanomaterials using artificial neural networks can potentially lead to a notable revolution in future materials design. Although progress has been made in creating small and simple molecules, complex materials such as crystalline porous materials have yet to be generated using any of the neural networks. Here, we have implemented a generative adversarial network that uses a training set of 31,713 known zeolites to produce 121 crystalline porous materials. Our neural network takes in inputs in the form of energy and material dimensions, and we show that zeolites with a user-desired range of 4 kJ/mol methane heat of adsorption can be reliably produced using our neural network. The fine-tuning of user-desired capability can potentially accelerate materials development as it demonstrates a successful case of inverse design of porous materials.

Spherical Sandwich Au@Pd@UIO-67/Pt@UIO- n ( n = 66, 67, 69) Core-Shell Catalysts: Zr-Based Metal-Organic Frameworks for Effectively Regulating the Reverse Water-Gas Shift Reaction

In this study, spherical sandwich Au@Pd@UIO-67/Pt@UIO- n ( n = 66, 67, 69) core-shell catalysts were assembled. Au nanoparticles (NPs) were used as the core for the epitaxial growth of Pd shells, and Au@Pd core-shell NPs were successfully encapsulated in the center of monodispersed Au@Pd@UIO-67 nanospheres. Pt NPs were fully fixed onto the nanosphere surfaces to obtain Au@Pd@UIO-67/Pt composites; further coating with UIO- n led to Au@Pd@UIO-67/Pt@UIO- n, in which Pt NPs are sandwiched between the Au@Pd@UIO-67 core and the UIO- n shell. The Au@Pd core-shell NPs efficiently controlled the morphology and structure of UIO-67 and enhanced the CO selectivity of the catalyst. Pt NPs increased the CO₂ conversion, and the UIO- n component effectively regulated the reverse water-gas shift reaction.

Automated Multiscale Approach To Predict Self-Diffusion from a Potential Energy Field

For large-scale screening studies there is a need to estimate the diffusion of gas molecules in nanoporous materials more efficiently than (brute force) molecular dynamics. In particular for systems with low diffusion coefficients molecular dynamics can be prohibitively expensive. An alternative is to compute the hopping rates between adsorption sites using transition state theory. For large-scale screening this requires the automatic detection of the transition states between the adsorption sites along the different diffusion paths. Here an algorithm is presented that analyzes energy grids for the moving particles. It detects the energies at which diffusion paths are formed, together with their directions. This allows for easy identification of nondiffusive systems. For diffusive systems, it partitions the grid coordinates assigned to energy basins and transitions states, permitting a transition state theory based analysis of the diffusion. We test our method on CH₄ diffusion in zeolites, using a standard kinetic Monte Carlo simulation based on the output of our grid analysis. We find that it is accurate, fast, and rigorous without limitations to the geometries of the diffusion tunnels or transition states.

Exceptional hydrogen storage achieved by screening nearly half a million metal-organic frameworks

Few hydrogen adsorbents balance high usable volumetric and gravimetric capacities. Although metal-organic frameworks (MOFs) have recently demonstrated progress in closing this gap, the large number of MOFs has hindered the identification of optimal materials. Here, a systematic assessment of published databases of real and hypothetical MOFs is presented. Nearly 500,000 compounds were screened computationally, and the most promising were assessed experimentally. Three MOFs with capacities surpassing that of IRMOF-20, the record-holder for balanced hydrogen capacity, are demonstrated: SNU-70, UMCM-9, and PCN-610/NU-100. Analysis of trends reveals the existence of a volumetric ceiling at ∼40 g H2 L-1. Surpassing this ceiling is proposed as a new capacity target for hydrogen adsorbents. Counter to earlier studies of total hydrogen uptake in MOFs, usable capacities in the highest-capacity materials are negatively correlated with density and volumetric surface area. Instead, capacity is maximized by increasing gravimetric surface area and porosity. This suggests that property/performance trends for total capacities may not translate to usable capacities.

Efficiently Exploring Adsorption Space to Identify Privileged Adsorbents for Chemical Separations of a Diverse Set of Molecules

Although computational models have been used to predict adsorption of molecules in large libraries of porous adsorbents, previous work of this kind has focused on a small number of molecules as potential adsorbates. In this study, molecular simulations were used to consider the adsorption of a diverse range of molecules in a large collection of metal-organic framework (MOF) materials. Specifically, 11 304 isotherms were obtained from molecular simulations of 24 different adsorbates in 471 MOFs. This information provides insight into several interesting questions that could not be addressed with previously available data. Highly computationally efficient methods are introduced that can predict isotherms for a wide range of adsorbing molecules with far less computation than traditional molecular simulations. By characterizing the 276 binary mixtures defined by the molecules considered, "privileged" adsorbents are shown to exist, which are effective for separating many different molecular mixtures. Finally, correlations that were developed previously to predict molecular solubility in polymers are found to be surprisingly effective in predicting the average properties of molecules adsorbing in MOFs.

Materials Genome in Action: Identifying the Performance Limits of Physical Hydrogen Storage

The Materials Genome is in action: the molecular codes for millions of materials have been sequenced, predictive models have been developed, and now the challenge of hydrogen storage is targeted. Renewably generated hydrogen is an attractive transportation fuel with zero carbon emissions, but its storage remains a significant challenge. Nanoporous adsorbents have shown promising physical adsorption of hydrogen approaching targeted capacities, but the scope of studies has remained limited. Here the Nanoporous Materials Genome, containing over 850 000 materials, is analyzed with a variety of computational tools to explore the limits of hydrogen storage. Optimal features that maximize net capacity at room temperature include pore sizes of around 6 Å and void fractions of 0.1, while at cryogenic temperatures pore sizes of 10 Å and void fractions of 0.5 are optimal. Our top candidates are found to be commercially attractive as "cryo-adsorbents", with promising storage capacities at 77 K and 100 bar with 30% enhancement to 40 g/L, a promising alternative to liquefaction at 20 K and compression at 700 bar.

Materials design by evolutionary optimization of functional groups in metal-organic frameworks

A genetic algorithm that efficiently optimizes a desired physical or functional property in metal-organic frameworks (MOFs) by evolving the functional groups within the pores has been developed. The approach has been used to optimize the CO₂ uptake capacity of 141 experimentally characterized MOFs under conditions relevant for postcombustion CO₂ capture. A total search space of 1.65 trillion structures was screened, and 1035 derivatives of 23 different parent MOFs were identified as having exceptional CO₂ uptakes of >3.0 mmol/g (at 0.15 atm and 298 K). Many well-known MOF platforms were optimized, with some, such as MIL-47, having their CO₂ adsorption increase by more than 400%. The structures of the high-performing MOFs are provided as potential targets for synthesis.

In silico discovery of metal-organic frameworks for precombustion CO2 capture using a genetic algorithm

Discovery of new adsorbent materials with a high CO2 working capacity could help reduce CO2 emissions from newly commissioned power plants using precombustion carbon capture. High-throughput computational screening efforts can accelerate the discovery of new adsorbents but sometimes require significant computational resources to explore the large space of possible materials. We report the in silico discovery of high-performing adsorbents for precombustion CO2 capture by applying a genetic algorithm to efficiently search a large database of metal-organic frameworks (MOFs) for top candidates. High-performing MOFs identified from the in silico search were synthesized and activated and show a high CO2 working capacity and a high CO2/H2 selectivity. One of the synthesized MOFs shows a higher CO2 working capacity than any MOF reported in the literature under the operating conditions investigated here.

Towards accurate porosity descriptors based on guest-host interactions

For nanoporous materials at the characterization level, geometry-based approaches have become the methods of choice to provide information, often encoded in numerical descriptors, about the pores and the channels of a porous material. Examples of most common descriptors of the latter are pore limiting diameters, accessible surface area and accessible volume. The geometry-based methods exploit hard-sphere approximation for atoms, which (1) reduces costly computations of the interatomic interactions between the probe guest molecule and the porous material framework atoms, (2) effectively exploit applied mathematics methods such as Voronoi decomposition to represent and characterize porosity. In this work, we revisit and quantify the shortcoming of the geometry-based approaches. To do so, we have developed a series of algorithms to calculate pore descriptors such as void fraction, accessible surface area, pore limiting diameters (largest included sphere, and largest free sphere) based on a classical force field model of interactions between the guest and the framework atoms. Our resulting energy-based methods are tested on diverse sets of metal-organic frameworks and zeolite structures and comparisons against results obtained from geometric-based method indicate deviations in the cases for structures with small pore sizes. The method provides both high accuracy and performance making it suitable when screening a large database of materials.