Engineering topics
Snurr, Randall Q.
Publications and source records attributed to Snurr, Randall Q..
Fingerprinting diverse nanoporous materials for optimal hydrogen storage conditions using meta-learning
A meta-learning model predicts the optimal hydrogen storage conditions and selects from a diverse set of nanoporous materials.
Machine learning using host/guest energy histograms to predict adsorption in metal–organic frameworks: Application to short alkanes and Xe/Kr mixtures
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. Here, 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.
Molecular fingerprint and machine learning to accelerate design of high-performance homochiral metal–organic frameworks
In this report computational screening was employed to calculate the enantioseparation capabilities of 45 functionalized homochiral metal–organic frameworks (FHMOFs), and machine learning (ML) and molecular fingerprint (MF) techniques were used to find new FHMOFs with high performance. With increasing temperature, the enantioselectivities for (R,S)-1,3-dimethyl-1,2-propadiene are improved. The “glove effect” in the chiral pockets was proposed to explain the correlations between the steric effect of functional groups and performance of FHMOFs. Moreover, the neighborhood component analysis and RDKit/MACCS MFs show the highest predictive effect on enantioselectivities among the four ML classification algorithms with nine MFs that were tested. Based on the importance of MF, 85 new FHMOFs were designed, and a newly designed FHMOF, NO 2 -NHOH-FHMOF, with high similarity to the optimal MFs achieved improved chiral separation performance, with enantioselectivities of 85%. The design principles and new chiral pockets obtained by ML and MFs could facilitate the development of new materials for chiral separation.
Selective Photodimerization in a Cyclodextrin Metal–Organic Framework
For the most part, enzymes contain one active site wherein they catalyze in a serial manner chemical reactions between substrates both efficiently and rapidly. Imagine if a situation could be created within a chiral porous crystal containing trillions of active sites where substrates can reside in vast numbers before being converted in parallel into products. Here, we report how it is possible to incorporate 1-anthracenecarboxylate (1-AC – ) as a substrate into a γ-cyclodextrin-containing metal–organic framework (CD-MOF-1), where the metals are K + cations, prior to carrying out [4+4] photodimerizations between pairs of substrate molecules, affording selectively one of four possible regioisomers. One of the high-yielding regioisomers exhibits optical activity as a result of the presence of an 8:1 ratio of the two enantiomers following separation by high-performance liquid chromatography. The solid-state superstructure of 1-anthracenecarboxylate potassium salt (1-ACK), which is co-crystallized with γ-cyclodextrin, reveals that pairs of substrate molecules are not only packed inside tunnels between spherical cavities present in CD-MOF-1, but also stabilized—in addition to hydrogen-bonding to the C-2 and C-3 hydroxyl groups on the d-glucopyranosyl residues present in the γ-cyclodextrin tori—by combinations of hydrophobic and electrostatic interactions between the carboxyl groups in 1-AC– and four K + cations on the waistline between the two γ-cyclodextrin tori in the tunnels. Here, these non-covalent bonding interactions result in preferred co-conformations that account for the highly regio- and enantioselective [4+4] cycloaddition during photoirradiation. Theoretical calculations, in conjunction with crystallography, support the regio- and stereochemical outcome of the photodimerization.
Zirconium Metal–Organic Frameworks Integrating Chloride Ions for Ammonia Capture and/or Chemical Separation
Not provided.
Machine learning the quantum-chemical properties of metal–organic frameworks for accelerated materials discovery
Not Available
Ammonia Capture within Zirconium Metal–Organic Frameworks: Reversible and Irreversible Uptake
Not provided.
Fast and Accurate Machine Learning Strategy for Calculating Partial Atomic Charges in Metal–Organic Frameworks
Computational high-throughput screening using molecular simulations is a powerful tool for identifying top-performing metal–organic frameworks (MOFs) for gas storage and separation applications. Accurate partial atomic charges are often required to model the electrostatic interactions between the MOF and the adsorbate, especially when the adsorption involves molecules with dipole or quadrupole moments such as water and CO 2 . Although ab initio methods can be used to calculate accurate partial atomic charges, these methods are impractical for screening large material databases because of the high computational cost. We developed a random forest machine learning model to predict the partial atomic charges in MOFs using a small yet meaningful set of features that represent both the elemental properties and the local environment of each atom. The model was trained and tested on a collection of about 320 000 density-derived electrostatic and chemical (DDEC) atomic charges calculated on a subset of the Computation-Ready Experimental Metal–Organic Framework (CoRE MOF-2019) database and separately on charge model 5 (CM5) charges. The model predicts accurate atomic charges for MOFs at a fraction of the computational cost of periodic density functional theory (DFT) and is found to be transferable to other porous molecular crystals and zeolites. In conclusion, a strong correlation is observed between the partial atomic charge and the average electronegativity difference between the central atom and its bonded neighbors.
Inverse design of nanoporous crystalline reticular materials with deep generative models
Reticular frameworks are crystalline porous materials that form via the self-assembly of molecular building blocks in different topologies, with many having desirable properties for gas storage, separation, catalysis, biomedical applications and so on. The notable variety of building blocks makes reticular chemistry both promising and challenging for prospective materials design. Here we propose an automated nanoporous materials discovery platform powered by a supramolecular variational autoencoder for the generative design of reticular materials. We demonstrate the automated design process with a class of metal–organic framework (MOF) structures and the goal of separating carbon dioxide from natural gas or flue gas. Our model shows high fidelity in capturing MOF structural features. Here, we show that the autoencoder has a promising optimization capability when jointly trained with multiple top adsorbent candidates identified for superior gas separation. MOFs discovered here are strongly competitive against some of the best-performing MOFs/zeolites ever reported.
Supramolecular Porous Assemblies of Atomically Precise Catalytically Active Cerium-Based Clusters
Abstract not provided
Zr 6 O 8 Node-Catalyzed Butene Hydrogenation and Isomerization in the Metal–Organic Framework NU-1000
Zirconium-based metal–organic frameworks (Zr-MOFs) have been increasingly studied over the past two decades as heterogeneous catalysts due to their synthetic tunability, well-defined nature, and chemical stability. In contrast to traditional zirconia-based heterogeneous catalysts, the community has assumed that Zr-MOFs are inert catalyst supports that do not participate directly in hydrocarbon transformations, such as olefin hydrogenation and isomerization. Here, we report that the Zr-MOF NU-1000 is capable of catalyzing olefin hydrogenation and isomerization, without any postsynthetic modifications, under a hydrogen atmosphere. We probe H 2 activation over the nodes of NU-1000 via spectroscopic and computational techniques revealing that H 2 dissociation can occur heterolytically across coordinatively unsaturated Zr sites and proximal hydroxide and μ3-oxo ligands. These results, along with catalytic experiments, suggest that H 2 activation results in node-supported zirconium hydrides capable of the hydrogenation and isomerization of 1-butene. When examining rate dependence on the partial pressure of H 2 , we observe first-order dependence for hydrogenation and half-order dependence for isomerization. Half-order H 2 rate dependence is consistent with a mechanism where both fragments of cleaved H 2 are active for 1-butene isomerization, suggesting that heterolytic cleavage generates acidic protons resulting in parallel, acid-, and hydride-catalyzed isomerization pathways. In conclusion, this work shows that Zr-MOFs have more diverse reactivity than the current literature may suggest and opens possibilities for ways in which Zr-MOFs can be used as heterogeneous catalysts and supports.
Tuning the Atrazine Binding Sites in an Indium-Based Flexible Metal–Organic Framework
Not provided.
Exploring the Effects of Node Topology, Connectivity, and Metal Identity on the Binding of Nerve Agents and Their Hydrolysis Products in Metal–Organic Frameworks
Recent studies have shown that metal–organic frameworks (MOFs) built from hexanuclear M(IV) oxide cluster nodes are effective catalysts for nerve agent hydrolysis, where the properties of the active sites on the nodes can strongly influence the reaction energetics. The connectivity and metal identity of these M6 nodes can be easily tuned, offering extensive opportunities for computational screening to predict promising new materials. Thus, we used density functional theory (DFT) to examine the effects of node topology, connectivity, and metal identity on the binding energies of multiple nerve agents and their corresponding hydrolysis products. By computing an optimization metric based on the relative binding strengths of key hydrolysis reaction species (water, agent, and bidentate-bound products), we predicted optimal M6 nodes for hydrolyzing specific nerve agent and simulant molecules, where our results are in qualitative agreement with observed experimental trends. This analysis highlighted the notion that no single metal or node topology is optimal for all possible organophosphates, suggesting that MOFs should be selected based on the agent of interest. Using the large amount of data generated from our DFT calculations, we then derived quantitative structure–activity relationship (QSAR) models to help explain the complex trends observed in the binding energies. Through linear regression, we identified the most important descriptors for describing the binding of nerve agents and their hydrolysis products to M6 nodes. These results suggested that both molecular and node properties, including both structural and chemical features, collectively contribute to the binding energetics. By performing a thorough statistical analysis, we showed that our QSAR models are capable of making quantitatively accurate binding energy predictions for nerve agents and their hydrolysis products in a wide variety of M(IV)-MOFs. In conclusion, the insights gained herein can be used to guide future experiments for the synthesis of MOFs with enhanced catalytic activity for organophosphate hydrolysis.
Topological Effects on Separation of Alkane Isomers in Metal-Organic Frameworks.
Abstract not provided.
Comparing GGA, GGA+ U , and meta-GGA functionals for redox-dependent binding at open metal sites in metal–organic frameworks
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 O 2 and N 2 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.
Topological effects on separation of alkane isomers in metal–organic frameworks
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. Here, 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.