Understanding the diversity of the metal-organic framework ecosystem.

Millions of distinct metal-organic frameworks (MOFs) can be made by combining metal nodes and organic linkers. At present, over 90,000 MOFs have been synthesized and over 500,000 predicted. This raises the question whether a new experimental or predicted structure adds new information. For MOF chemists, the chemical design space is a combination of pore geometry, metal nodes, organic linkers, and functional groups, but at present we do not have a formalism to quantify optimal coverage of chemical design space. In this work, we develop a machine learning method to quantify similarities of MOFs to analyse their chemical diversity. This diversity analysis identifies biases in the databases, and we show that such bias can lead to incorrect conclusions. The developed formalism in this study provides a simple and practical guideline to see whether new structures will have the potential for new insights, or constitute a relatively small variation of existing structures.

Data-driven design of metal-organic frameworks for wet flue gas CO2 capture.

Limiting the increase of CO2 in the atmosphere is one of the largest challenges of our generation1. Because carbon capture and storage is one of the few viable technologies that can mitigate current CO2 emissions2, much effort is focused on developing solid adsorbents that can efficiently capture CO2 from flue gases emitted from anthropogenic sources3. One class of materials that has attracted considerable interest in this context is metal-organic frameworks (MOFs), in which the careful combination of organic ligands with metal-ion nodes can, in principle, give rise to innumerable structurally and chemically distinct nanoporous MOFs. However, many MOFs that are optimized for the separation of CO2 from nitrogen4-7 do not perform well when using realistic flue gas that contains water, because water competes with CO2 for the same adsorption sites and thereby causes the materials to lose their selectivity. Although flue gases can be dried, this renders the capture process prohibitively expensive8,9. Here we show that data mining of a computational screening library of over 300,000 MOFs can identify different classes of strong CO2-binding sites-which we term 'adsorbaphores'-that endow MOFs with CO2/N2 selectivity that persists in wet flue gases. We subsequently synthesized two water-stable MOFs containing the most hydrophobic adsorbaphore, and found that their carbon-capture performance is not affected by water and outperforms that of some commercial materials. Testing the performance of these MOFs in an industrial setting and consideration of the full capture process-including the targeted CO2 sink, such as geological storage or serving as a carbon source for the chemical industry-will be necessary to identify the optimal separation material.

Nucleobase pairing and photodimerization in a biologically derived metal-organic framework nanoreactor.

Biologically derived metal-organic frameworks (bio-MOFs) are of great importance as they can be used as models for bio-mimicking and in catalysis, allowing us to gain insights into how large biological molecules function. Through rational design, here we report the synthesis of a novel bio-MOF featuring unobstructed Watson-Crick faces of adenine (Ade) pointing towards the MOF cavities. We show, through a combined experimental and computational approach, that thymine (Thy) molecules diffuse through the pores of the MOF and become base-paired with Ade. The Ade-Thy pair binding at 40-45% loading reveals that Thy molecules are packed within the channels in a way that fulfill both the Woodward-Hoffmann and Schmidt rules, and upon UV irradiation, Thy molecules dimerize into Thy<>Thy. This study highlights the utility of accessible functional groups within the pores of MOFs, and their ability to 'lock' molecules in specific positions that can be subsequently dimerized upon light irradiation, extending the use of MOFs as nanoreactors for the synthesis of molecules that are otherwise challenging to isolate.

Evaluating Charge Equilibration Methods To Generate Electrostatic Fields in Nanoporous Materials.

Charge equilibration (Qeq) methods can estimate the electrostatic potential of molecules and periodic frameworks by assigning point charges to each atom, using only a small fraction of the resources needed to compute density functional (DFT)-derived charges. This makes possible, for example, the computational screening of thousands of microporous structures to assess their performance for the adsorption of polar molecules. Recently, different variants of the original Qeq scheme were proposed to improve the quality of the computed point charges. One focus of this research was to improve the gas adsorption predictions in metal-organic frameworks (MOFs), for which many different structures are available. In this work, we review the evolution of the method from the original Qeq scheme, understanding the role of the different modifications on the final output. We evaluated the result of combining different protocols and set of parameters, by comparing the Qeq charges with high quality DFT-derived DDEC charges for 2338 MOF structures. We focused on the systematic errors that are attributable to specific atom types to quantify the final precision that one can expect from Qeq methods in the context of gas adsorption where the electrostatic potential plays a significant role, namely, CO2 and H2S adsorption. In conclusion, both the type of algorithm and the input parameters have a large impact on the resulting charges, and we draw some guidelines to help the user to choose the proper combination of the two for obtaining a meaningful set of charges. We show that, considering this set of MOFs, the accuracy of the original Qeq scheme is often still comparable with the most recent variants, even if it clearly fails in the presence of certain atom types, such as alkali metals.

Biporous Metal-Organic Framework with Tunable CO2/CH4 Separation Performance Facilitated by Intrinsic Flexibility.

In this work, we report the synthesis of SION-8, a novel metal-organic framework (MOF) based on Ca(II) and a tetracarboxylate ligand TBAPy4- endowed with two chemically distinct types of pores characterized by their hydrophobic and hydrophilic properties. By altering the activation conditions, we gained access to two bulk materials: the fully activated SION-8F and the partially activated SION-8P with exclusively the hydrophobic pores activated. SION-8P shows high affinity for both CO2 ( Qst = 28.4 kJ/mol) and CH4 ( Qst = 21.4 kJ/mol), while upon full activation, the difference in affinity for CO2 ( Qst = 23.4 kJ/mol) and CH4 ( Qst = 16.0 kJ/mol) is more pronounced. The intrinsic flexibility of both materials results in complex adsorption behavior and greater adsorption of gas molecules than if the materials were rigid. Their CO2/CH4 separation performance was tested in fixed-bed breakthrough experiments using binary gas mixtures of different compositions and rationalized in terms of molecular interactions. SION-8F showed a 40-160% increase (depending on the temperature and the gas mixture composition probed) of the CO2/CH4 dynamic breakthrough selectivity compared to SION-8P, demonstrating the possibility to rationally tune the separation performance of a single MOF by manipulating the stepwise activation made possible by the MOF's biporous nature.

Improving the Mechanical Stability of Metal-Organic Frameworks Using Chemical Caryatids.

Metal-organic frameworks (MOFs) have emerged as versatile materials for applications ranging from gas separation and storage, catalysis, and sensing. The attractive feature of MOFs is that, by changing the ligand and/or metal, they can be chemically tuned to perform optimally for a given application. In most, if not all, of these applications one also needs a material that has a sufficient mechanical stability, but our understanding of how changes in the chemical structure influence mechanical stability is limited. In this work, we rationalize how the mechanical properties of MOFs are related to framework bonding topology and ligand structure. We illustrate that the functional groups on the organic ligands can either enhance the mechanical stability through formation of a secondary network of nonbonded interactions or soften the material by destabilizing the bonded network of a MOF. In addition, we show that synergistic effect of the bonding network of the material and the secondary network is required to achieve optimal mechanical stability of a MOF. The developed molecular insights in this work can be used for systematic improvement of the mechanical stability of the materials by careful selection of the functional groups.

Shedding Light on the Protonation States and Location of Protonated N Atoms of Adenine in Metal-Organic Frameworks.

We report the syntheses and structures of five metal-organic frameworks (MOFs) based on transition metals (NiII, CuII, and ZnII), adenine, and di-, tri-, and tetra-carboxylate ligands. Adenine, with multiple N donor sites, was found to coordinate to the metal centers in different binding modes including bidentate (through N7 and N9, or N3 and N9) and tridentate (through N3, N7, and N9). Systematic investigations of the protonation states of adenine in each MOF structure via X-ray photoelectron spectroscopy revealed that adenine can be selectively protonated through N1, N3, or N7. The positions of H atoms connected to the N atoms were found from the electron density maps, and further supported by the study of C-N-C bond angles compared to the literature reports. DFT calculations were performed to geometrically optimize and energetically assess the structures simulated with different protonation modes. The present study highlights the rich coordination chemistry of adenine and provides a method for the determination of its protonation states and the location of protonated N atoms of adenine within MOFs, a task that would be challenging in complicated adenine-based MOF structures.

Text Mining Metal-Organic Framework Papers.

We have developed a simple text mining algorithm that allows us to identify surface area and pore volumes of metal-organic frameworks (MOFs) using manuscript html files as inputs. The algorithm searches for common units (e.g., m2/g, cm3/g) associated with these two quantities to facilitate the search. From the sample set data of over 200 MOFs, the algorithm managed to identify 90% and 88.8% of the correct surface area and pore volume values. Further application to a test set of randomly chosen MOF html files yielded 73.2% and 85.1% accuracies for the two respective quantities. Most of the errors stem from unorthodox sentence structures that made it difficult to identify the correct data as well as bolded notations of MOFs (e.g., 1a) that made it difficult identify its real name. These types of tools will become useful when it comes to discovering structure-property relationships among MOFs as well as collecting a large set of data for references.

Accurate Characterization of the Pore Volume in Microporous Crystalline Materials.

Pore volume is one of the main properties for the characterization of microporous crystals. It is experimentally measurable, and it can also be obtained from the refined unit cell by a number of computational techniques. In this work, we assess the accuracy and the discrepancies between the different computational methods which are commonly used for this purpose, i.e, geometric, helium, and probe center pore volumes, by studying a database of more than 5000 frameworks. We developed a new technique to fully characterize the internal void of a microporous material and to compute the probe-accessible and -occupiable pore volume. We show that, unlike the other definitions of pore volume, the occupiable pore volume can be directly related to the experimentally measured pore volumes from nitrogen isotherms.

The Influence of Intrinsic Framework Flexibility on Adsorption in Nanoporous Materials.

For applications of metal-organic frameworks (MOFs) such as gas storage and separation, flexibility is often seen as a parameter that can tune material performance. In this work we aim to determine the optimal flexibility for the shape selective separation of similarly sized molecules (e.g., Xe/Kr mixtures). To obtain systematic insight into how the flexibility impacts this type of separation, we develop a simple analytical model that predicts a material's Henry regime adsorption and selectivity as a function of flexibility. We elucidate the complex dependence of selectivity on a framework's intrinsic flexibility whereby performance is either improved or reduced with increasing flexibility, depending on the material's pore size characteristics. However, the selectivity of a material with the pore size and chemistry that already maximizes selectivity in the rigid approximation is continuously diminished with increasing flexibility, demonstrating that the globally optimal separation exists within an entirely rigid pore. Molecular simulations show that our simple model predicts performance trends that are observed when screening the adsorption behavior of flexible MOFs. These flexible simulations provide better agreement with experimental adsorption data in a high-performance material that is not captured when modeling this framework as rigid, an approximation typically made in high-throughput screening studies. We conclude that, for shape selective adsorption applications, the globally optimal material will have the optimal pore size/chemistry and minimal intrinsic flexibility even though other nonoptimal materials' selectivity can actually be improved by flexibility. Equally important, we find that flexible simulations can be critical for correctly modeling adsorption in these types of systems.

Force-Field Prediction of Materials Properties in Metal-Organic Frameworks.

In this work, MOF bulk properties are evaluated and compared using several force fields on several well-studied MOFs, including IRMOF-1 (MOF-5), IRMOF-10, HKUST-1, and UiO-66. It is found that, surprisingly, UFF and DREIDING provide good values for the bulk modulus and linear thermal expansion coefficients for these materials, excluding those that they are not parametrized for. Force fields developed specifically for MOFs including UFF4MOF, BTW-FF, and the DWES force field are also found to provide accurate values for these materials' properties. While we find that each force field offers a moderately good picture of these properties, noticeable deviations can be observed when looking at properties sensitive to framework vibrational modes. This observation is more pronounced upon the introduction of framework charges.

An unprecedented Co(II) cuboctahedron as the secondary building unit in a Co-based metal-organic framework.

A cubic metal-organic framework with an unprecedented octanuclear secondary building unit (SBU) was isolated. The obtained SBU is composed of 8 Co(II) ions at each vertex, 6 µ4-OH groups at each face, and 12 cpt(-) ligands framing the metal core. The cuboctahedra arrange in a ubt framework topology, eliciting a highly symmetrical MOF structure. Magnetic measurements as well as DFT calculations on this crystalline MOF reveal intramolecular antiferromagnetic coupling between Co(II) ions in the octanuclear SBU.

Rapid and Accurate Machine Learning Recognition of High Performing Metal Organic Frameworks for CO2 Capture.

In this work, we have developed quantitative structure-property relationship (QSPR) models using advanced machine learning algorithms that can rapidly and accurately recognize high-performing metal organic framework (MOF) materials for CO2 capture. More specifically, QSPR classifiers have been developed that can, in a fraction of a section, identify candidate MOFs with enhanced CO2 adsorption capacity (>1 mmol/g at 0.15 bar and >4 mmol/g at 1 bar). The models were tested on a large set of 292 050 MOFs that were not part of the training set. The QSPR classifier could recover 945 of the top 1000 MOFs in the test set while flagging only 10% of the whole library for compute intensive screening. Thus, using the machine learning classifiers as part of a high-throughput screening protocol would result in an order of magnitude reduction in compute time and allow intractably large structure libraries and search spaces to be screened.

Direct observation and quantification of CO2 binding within an amine-functionalized nanoporous solid.

Understanding the molecular details of CO(2)-sorbent interactions is critical for the design of better carbon-capture systems. Here we report crystallographic resolution of CO(2) molecules and their binding domains in a metal-organic framework functionalized with amine groups. Accompanying computational studies that modeled the gas sorption isotherms, high heat of adsorption, and CO(2) lattice positions showed high agreement on all three fronts. The modeling apportioned specific binding interactions for each CO(2) molecule, including substantial cooperative binding effects among the guest molecules. The validation of the capacity of such simulations to accurately model molecular-scale binding bodes well for the theory-aided development of amine-based CO(2) sorbents. The analysis shows that the combination of appropriate pore size, strongly interacting amine functional groups, and the cooperative binding of CO(2) guest molecules is responsible for the low-pressure binding and large uptake of CO(2) in this sorbent material.