Accurate Ab Initio Description of Adsorption on Coordinatively Unsaturated Cu2+ and Fe3+ Sites in MOFs

Survey on Adsorption of Low Molecular Weight Compounds in Cu-BTC Metal-Organic Framework: Experimental Results and Thermodynamic Modeling

This contribution aims at providing a critical overview of experimental results for the sorption of low molecular weight compounds in the Cu-BTC Metal-Organic Framework (MOF) and of their interpretation using available and new, specifically developed, theoretical approaches. First, a literature review of experimental results for the sorption of gases and vapors is presented, with particular focus on the results obtained from vibrational spectroscopy techniques. Then, an overview of theoretical models available in the literature is presented starting from semiempirical theoretical approaches suitable to interpret the adsorption thermodynamics of gases and vapors in Cu-BTC. A more detailed description is provided of a recently proposed Lattice Fluid approach, the Rigid Adsorbent Lattice Fluid (RALF) model. In addition, to deal with the cases where specific self- and cross-interactions (e.g., H-bonding, Lewis acid/Lewis base interactions) play a role, a modification of the RALF model, i.e., the RALFHB model, is introduced here for the first time. An extension of both RALF and RALFHB is also presented to cope with the cases in which the heterogeneity of the rigid adsorbent displaying a different kind of adsorbent cages is of relevance, as it occurs for the adsorption of some low molecular weight substances in Cu-BTC MOF.

Dressing of Cu Atom over Nickel Cluster Stimulating the Poisoning-Free CO Oxidation: An Ab Initio Study

In this work using first-principles calculations based on spin-polarized density functional theory (DFT), the role of the Cu atom in degrading the poisoning of carbon monoxide (CO) over Ni_nCu clusters is unveiled. The search has been initiated with the examination of structural, magnetic, and electronic properties of Ni_n+1 and Ni_nCu clusters (1 ≤ n ≤ 12). X-ray absorption near-edge structure (XANES) spectra of Ni K-edge are computed to extract the information on the oxidation states and coordination environment of metal sites of the clusters. This study is operated with the two forms of dispersion corrections, i.e., D2 and D3, with standard DFT (with LDA and GGA functionals) for the consideration of van der Waals interactions during CO adsorption. The PBE and PBE-D3 approaches are found to be capable of yielding the experimentally observed preferential site for CO adsorption. The effect of spin-polarization on the reactivity of transition metals (TMs) toward CO adsorption is crucially assessed by the electronic reactivity descriptors such as d-band center, d-band width, and fractional filling of d-band using a spin-polarized d-band center model. The effective charge transfer from Cu to Ni atoms makes Ni atoms more efficient of charge and is attributed to the degrading adsorption of CO over Ni_nCu clusters. The Ni₁₂Cu cluster stands out with good CO oxidation activity for the Langmuir-Hinshelwood (L-H) reaction pathway.

Importance of Lattice Constants in QM/MM Calculations on Metal-Organic Frameworks

Metal-organic frameworks (MOFs) are crystalline materials with novel physical and chemical properties. Computational simulations have become powerful complements to experiment for understanding catalysis in MOFs and developing new MOFs and their applications. However, due to their relatively large and complex structures, MOFs can be burdensome for fully quantum mechanical calculations. A combined quantum mechanical and molecular mechanical (QM/MM) method that combines the accuracy of fully quantum mechanical methods and the efficiency of MM methods is therefore attractive. In this study, we employ a QM/MM method for the study of two classes of chemical process in a MOF: the conversion of reaction intermediates in an Ir-containing borylation catalyst supported on MOF UiO-67 and the diffusion of a diborylated methane molecule in the pristine UiO-67 framework. We compare the QM/MM results with full-quantum mechanical results on large systems to validate the accuracy of the applied QM/MM method. In the first case, we consider a model of the entire system by partitioning it into subsystems that interact covalently, and in the second case the subsystem interaction is mainly steric. We observe that the QM/MM results agree with the full-quantum mechanical results within an average of 4 kcal/mol in the first case with strong electronic interactions and within an average of 3 kcal/mol in the case with only noncovalent interactions. An important lesson learned from the present study is that the quantitative results are very sensitive to the lattice constants predicted by the MM method used in the QM/MM calculations.

Construction of an Anion-Pillared MOF Database and the Screening of MOFs Suitable for Xe/Kr Separation

The separation of xenon/krypton (Xe/Kr) mixtures is a challenging process. Many porous materials allow the adsorption of both Xe and Kr but only with low selectivity. Anion-pillared metal-organic frameworks (MOFs), featuring the anion groups as structural pillars, show potential in gas separations, but only a limited number of them have been synthesized. Here, we describe a collection of 936 anion-pillared MOFs based on 22 experimentally available structures. We performed density functional theory (DFT) optimization and then assigned density-derived electrostatic and chemical (DDEC) charges for each MOF to make them well suited to molecular simulations. The structural properties of the MOFs vary more strongly with the choice of the organic ligand than with other aspects like fluorine groups and metal centers. We then screened the entire collection of MOFs in the context of Xe/Kr separation at room temperature. Compared with previously reported MOFs, the interpenetrated MOF SIFSIX-6-Cd-i is predicted to perform better for Xe/Kr separations, with a good balance between working capacity (1.62 mmol/g) and separation selectivity (16.4) at 298 K and 100 kPa. We also found that the heterogeneity of fluorine groups within a MOF can help to enhance Xe working capacity without reducing the Xe/Kr selectivity, suggesting that synthesis of anion-pillared MOFs with mixed fluorine groups may lead to improved Xe/Kr separation performance.

Porous Metal-Organic Polyhedra: Morphology, Porosity, and Guest Binding

Designing porous materials which can selectively adsorb CO₂ or CH₄ is an important environmental and industrial goal which requires an understanding of the host–guest interactions involved at the atomic scale. Metal–organic polyhedra (MOPs) showing permanent porosity upon desolvation are rarely observed. We report a family of MOPs (Cu-1a, Cu-1b, Cu-2), which derive their permanent porosity from cavities between packed cages rather than from within the polyhedra. Thus, for Cu-1a, the void fraction outside the cages totals 56% with only 2% within. The relative stabilities of these MOP structures are rationalized by considering their weak nondirectional packing interactions using Hirshfeld surface analyses. The exceptional stability of Cu-1a enables a detailed structural investigation into the adsorption of CO₂ and CH₄ using in situ X-ray and neutron diffraction, coupled with DFT calculations. The primary binding sites for adsorbed CO₂ and CH₄ in Cu-1a are found to be the open metal sites and pockets defined by the faces of phenyl rings. More importantly, the structural analysis of a hydrated sample of Cu-1a reveals a strong hydrogen bond between the adsorbed CO₂ molecule and the Cu(II)-bound water molecule, shedding light on previous empirical and theoretical observations that partial hydration of metal−organic framework (MOF) materials containing open metal sites increases their uptake of CO₂. The results of the crystallographic study on MOP–gas binding have been rationalized using DFT calculations, yielding individual binding energies for the various pore environments of Cu-1a.

We report a family of metal−organic polyhedra (MOP), which derive their permanent porosity from cavities between packed cages rather than from within the polyhedra. The relative stabilities of these MOP structures are rationalized by considering their weak nondirectional packing interactions using Hirshfeld surface analysis. A detailed structural investigation into the adsorption of CO₂ and CH₄ is reported using in situ X-ray and neutron diffraction, coupled with DFT calculations.

MOF-Based Adsorbents for Atmospheric Emission Control: A Review

This review focuses on the use of metal–organic frameworks (MOFs) for adsorbing gas species that are known to weaken the thermal self-regulation capacities of Earth’s atmosphere. A large section is dedicated to the adsorption of carbon dioxide, while another section is dedicated to the adsorption of other different gas typologies, whose emissions, for various reasons, represent a “wound” for Earth’s atmosphere. High emphasis is given to MOFs that have moved enough ahead in their development process to be currently considered as potentially usable in “real-world” (i.e., out-of-lab) adsorption processes. As a result, there is strong evidence of a wide gap between laboratory results and the industrial implementation of MOF-based adsorbents. Indeed, when a MOF that performs well in a specific process is commercially available in large quantities, economic observations still make designers tend toward more traditional adsorbents. Moreover, there are cases in which a specific MOF remarkably outperforms the currently employed adsorbents, but it is not industrially produced, thus strongly limiting its possibilities in large-scale use. To overcome such limitations, it is hoped that the chemical industry will be able to provide more and more mass-produced MOFs at increasingly competitive costs in the future.

Coordinatively unsaturated metal sites (open metal sites) in metal–organic frameworks: design and applications

The defined synthesis of OMS in MOFs is the basis for targeted functionalization through grafting, the coordination of weakly binding species and increased (supramolecular) interactions with guest molecules.

Theoretical Study of Methane Storage in Cu24(m-BDC)24

Calculations on the Cu₂₄(m-BDC)₂₄ (m-BDC = 1,3-benzenedicarboxylate) polyoxometalate (POM) cage with 0, 12, 24, and 40 methane molecules inside were made using the M06 exchange/correlation functional. During filling of the cage with 40 CH₄ molecules, the 12 strongest binding CH₄ molecules are those to the coordination unsaturated sites (CUS) to the inwardly directed Cu(+2) centers via agostic interactions. The next 12 CH₄ molecules are less tightly bound followed by the next 16 CH₄ molecules with average binding energies of 8.27, 7.88, and 7.36 kcal/mol per CH₄, respectively. A section of the Cu₂₄(m-BDC)₂₄ cage was taken with the formula Cu₄(m-BDC)(BC)₆ (BC = benezenecarboxylate) in order to estimate zero-point, thermal, and entropy corrections of the larger cage. Estimating free energies at 1 bar, the Cu₂₄(m-BDC)₂₄ POM is predicted to lose 16, 12, and 12 CH₄ molecules at 67, 123, and 171 °C, respectively. The 40CH₄@Cu₂₄(m-BDC)₂₄ cage, which is isostructural to the main cavity of HKUST-1 with 40 CH₄ molecules inside, is predicted to have a loading of 224 cm³(STP) cm^-3 at 1 bar.

Ab Initio Evaluation of Henry Coefficients Using Importance Sampling

We present a new algorithm that allows for an efficient evaluation of the Henry coefficient of a guest molecule inside a porous material, which permits to use ab initio energy calculations. The Widom insertion method, which is currently used to compute these Henry coefficients, typically requires millions of energy evaluations. Our new methodology reduces this number by more than 1 order of magnitude, enabling the use of an ab initio potential energy surface. The methodology we propose is reminiscent of the well-known importance sampling technique which is frequently used in Monte Carlo integrations. First, a conventional Widom insertion simulation is performed using a force field. In the second step, the Widom results are used to select a limited number of configurations and only for these configurations the ab initio evaluation of the energy is required. Finally, by appropriately reweighting the latter energies, an accurate estimation of the ab initio Henry coefficient is possible at a moderate computational cost. We apply our methodology to the adsorption of CO₂ in Mg-MOF-74, a prototypical case where interactions of a polar guest molecule with unsaturated metal sites dominate the adsorption mechanism. In this case generic force fields such as UFF or Dreiding are inappropriate and the use of ab initio methods is indispensable. In a second case study, we compute Henry coefficients of methane in UiO-66 using different levels of theory. We pay particular attention to the influence of the dispersion corrections and the role of many-body effects. For the final example, we qualitatively investigate adsorption features for a series of functionalized UiO-66 frameworks. Overall the cases we present show that accurate computations of Henry coefficients is extremely challenging, as different levels of theory provide strongly varying results. At the same time ab initio calculations have added value compared to force fields, as they provide a physically more sound description of the adsorption mechanism and in some cases clearly improve correspondence with experiment.

Flue gas adsorption on periodic mesoporous phenylene-silica: a DFT approach

Periodic mesoporous organosilicas (PMOs) were suggested as potential adsorbents for CO2/CH4 separation because of their large affinities towards CO2 and low interaction with CH4. Herewith, we present a comprehensive computational study on the binding properties of flue gas species with the pore walls of periodic mesoporous phenylene-silica (Ph-PMO) for understanding the possible impact of other gaseous species in the CO2/CH4 separation. The calculations considered three exchange-correlation functionals (PBE, PBE-D2 and M06-2X) based on the density functional theory and the walls of the periodic mesoporous phenylene-silica were modelled within the cluster model approach. The components of the flue gas considered were the diatomic CO, H2, N2, O2 and NO molecules, the triatomic CO2, H2O, H2S and SO2 species, the tetratomic SO3 and NH3 gases and the pentatomic CH4 molecule. The calculated data demonstrate that the presence of H2O, SO2, NH3, H2S and SO3 is a significant threat to CO2 capture by Ph-PMO and suggest that the Ph-PMO material would present high selectivity for CO2 over CH4, CO, H2 or N2 adsorption. The adsorption behaviour of flue gas components in Ph-PMO can be directly related to the experimental proton affinities, basicities or even the polarizabilities of the gaseous molecules.

Origin of the Strong Interaction between Polar Molecules and Copper(II) Paddle-Wheels in Metal Organic Frameworks

The copper paddle-wheel is the building unit of many metal organic frameworks. Because of the ability of the copper cations to attract polar molecules, copper paddle-wheels are promising for carbon dioxide adsorption and separation. They have therefore been studied extensively, both experimentally and computationally. In this work we investigate the copper-CO2 interaction in HKUST-1 and in two different cluster models of HKUST-1: monocopper Cu(formate)2 and dicopper Cu2(formate)4. We show that density functional theory methods severely underestimate the interaction energy between copper paddle-wheels and CO2, even including corrections for the dispersion forces. In contrast, a multireference wave function followed by perturbation theory to second order using the CASPT2 method correctly describes this interaction. The restricted open-shell Møller-Plesset 2 method (ROS-MP2, equivalent to (2,2) CASPT2) was also found to be adequate in describing the system and used to develop a novel force field. Our parametrization is able to predict the experimental CO2 adsorption isotherms in HKUST-1, and it is shown to be transferable to other copper paddle-wheel systems.

Metal-organic and covalent organic frameworks as single-site catalysts

Heterogeneous single-site catalysts consist of isolated, well-defined, active sites that are spatially separated in a given solid and, ideally, structurally identical. In this review, the potential of metal-organic frameworks (MOFs) and covalent organic frameworks (COFs) as platforms for the development of heterogeneous single-site catalysts is reviewed thoroughly. In the first part of this article, synthetic strategies and progress in the implementation of such sites in these two classes of materials are discussed. Because these solids are excellent playgrounds to allow a better understanding of catalytic functions, we highlight the most important recent advances in the modelling and spectroscopic characterization of single-site catalysts based on these materials. Finally, we discuss the potential of MOFs as materials in which several single-site catalytic functions can be combined within one framework along with their potential as powerful enzyme-mimicking materials. The review is wrapped up with our personal vision on future research directions.

CO2 Capture and Separations Using MOFs: Computational and Experimental Studies

This Review focuses on research oriented toward elucidation of the various aspects that determine adsorption of CO₂ in metal-organic frameworks and its separation from gas mixtures found in industrial processes. It includes theoretical, experimental, and combined approaches able to characterize the materials, investigate the adsorption/desorption/reaction properties of the adsorbates inside such environments, screen and design new materials, and analyze additional factors such as material regenerability, stability, effects of impurities, and cost among several factors that influence the effectiveness of the separations. CO₂ adsorption, separations, and membranes are reviewed followed by an analysis of the effects of stability, impurities, and process operation conditions on practical applications.