Alternative Pathway of CO2 Hydrogenation by Lewis-Pair-Functionalized UiO-66 MOF Revealed by Metadynamics Simulations

Ab Initio Molecular Dynamics Investigation of Water and Butanone Adsorption on UiO-66 with Defects

Volatile organic compounds (VOCs) are harmful chemicals that are found in minute quantities in the atmosphere and are emitted from a variety of industrial and biological processes. They can be harmful to breathe or serve as biomarkers for disease detection. Therefore, capture and detection of VOCs is important. Here, we have examined if the Zr-based UiO-66 metal-organic framework (MOF) can be used to capture butanone─a well-known VOC. Toward that end, we have performed Born-Oppenheimer ab initio molecular dynamics (AIMD) at 300 and 500 K to probe the energetics and molecular interactions between butanone [CH3C(O)CH2CH3] and open-cage Zr-UiO-66. Such interactions were systematically interrogated using three MOF structures: defective MOF with a missing 1,4-benzene-dicarboxylate linker and two H2O; pristine MOF with two H2O; and pristine dry MOF. These structures were loaded with one and four molecules of butanone to examine the effect of concentration as well. One-molecule loading interacted favorably with the defective structure at 300 K, only. In comparison, interactions with four-molecule loading were energetically favorable for all conditions. Persistent hydrogen bonds between the O atom of butanone, H2O, and μ3-OH hydroxyl attachments at Zr nodes substantially contributed to the intermolecular interactions. At higher loadings, butanone also showed a pronounced tendency to diffuse into the adjoining cages of Zr-UiO-66. The effect of such movement on interaction energies was rationalized using simple statistical mechanics-based models of interacting and noninteracting gases. Broadly, we learn that the presence of prior moisture within the interstitial cages of Zr-UiO-66 significantly impacts the adsorption behavior of butanone.

Experimental and computational aspects of molecular frustrated Lewis pairs for CO2 hydrogenation: en route for heterogeneous systems?

Catalysis plays a crucial role in advancing sustainability. The unique reactivity of frustrated Lewis pairs (FLPs) is driving an ever-growing interest in the transition metal-free transformation of small molecules like CO2 into valuable products. In this area, there is a recent growing incentive to heterogenize molecular FLPs into porous solids, merging the benefits of homogeneous and heterogeneous catalysis - high activity, selectivity, and recyclability. Despite the progress, challenges remain in preventing deactivation, poisoning, and simplifying catalyst-product separation. This review explores the expanding field of FLPs in catalysis, covering existing molecular FLPs for CO2 hydrogenation and recent efforts to design heterogeneous porous systems from both experimental and theoretical perspectives. Section 2 discusses experimental examples of CO2 hydrogenation by molecular FLPs, starting with stoichiometric reactions and advancing to catalytic ones. It then examines attempts to immobilize FLPs in porous matrices, including siliceous solids, metal-organic frameworks (MOFs), covalent organic frameworks, and disordered polymers, highlighting current limitations and challenges. Section 3 then reviews computational studies on the mechanistic details of CO2 hydrogenation, focusing on H2 splitting and hydride/proton transfer steps, summarizing efforts to establish structure-activity relationships. It also covers the computational aspects on grafting FLPs inside MOFs. Finally, Section 4 summarizes the main design principles established so far, while addressing the complexities of translating computational approaches into the experimental realm, particularly in heterogeneous systems. This section underscores the need to strengthen the dialogue between theoretical and experimental approaches in this field.

Enhancing CO2 Hydrogenation Using a Heterogeneous Bimetal NiAl-Deposited Metal-Organic Framework NU-1000: Insights from First-Principles Calculations

The hydrogenation of CO2 to high-value-added liquid fuels is crucial for greenhouse gas emission reduction and optimal utilization of carbon resources. Developing supported heterogeneous catalysts is a key strategy in this context, as they offer well-defined active sites for in-depth mechanistic studies and improved catalyst design. Here, we conducted extensive first-principles calculations to systematically explore the reaction mechanisms for CO2 hydrogenation on a heterogeneous bimetal NiAl-deposited metal-organic framework (MOF) NU-1000 and its catalytic performance as atomically dispersed catalysts for CO2 hydrogenation to formic acid (HCOOH), formaldehyde (H2CO), and methanol (CH3OH). The present results reveal that the presence of the NiAl-oxo cluster deposited on NU-1000 efficiently activates H2, and the facile heterolysis of H2 on Ni and adjacent O sites serves as a precursor to the hydrogenation of CO2 into various C1 products HCOOH, H2CO, and CH3OH. Generally, H2 activation is the rate-determining step in the entire CO2 hydrogenation process, the corresponding relatively low free energy barriers range from 14.5 to 15.9 kcal/mol, and the desorption of products on NiAl-deposited NU-1000 is relatively facile. Although the Al atom does not directly participate in the reaction, its presence provides exposed oxygen sites that facilitate the heterolytic cleavage of H2 and the hydrogenation of C1 intermediates, which plays an important role in enhancing the catalytic activity of the Ni site. The present study demonstrates that the catalytic performance of NU-1000 can be finely tuned by depositing heterometal-oxo clusters, and the porous MOF should be an attractive platform for the construction of atomically dispersed catalysts.

Catalytic Ring-Opening Polymerisation of Cyclic Ethylene Carbonate: Importance of Elementary Steps for Determining Polymer Properties Revealed via DFT-MTD Simulations Validated Using Kinetic Measurements

The production of CO2-containing polymers is still very demanding in terms of controlling the synthesis of products with pre-defined CO2 content and molecular weight. An elegant way of synthesising these polymers is via CO2-containing building blocks, such as cyclic ethylene carbonate (cEC), via catalytic ring-opening polymerisation. However, to date, the mechanism of this reaction and control parameters have not been elucidated. In this work, using DFT-metadynamics simulations for exploiting the potential of the polymerisation process, we aim to shed more light on the mechanisms of the interaction between catalysts (in particular, the catalysts K3VO4, K3PO4, and Na2SnO3) and the cEC monomer in the propagation step of the polymeric chain and the occurring CO2 release. Confirming the simulation results via subsequent kinetics measurements indicates that, depending on the catalyst's characteristics, it can be attached reversibly to the polymeric chain during polymerisation, resulting in a defined lifetime of the activated polymer chain. The second anionic oxygen of the catalyst can promote the catalyst's transfer to another electrophilic cEC monomer, terminating the growth of the first chain and initiating the propagation of the new polymer chain. This transfer reaction is an essential step in controlling the molecular weight of the products.

Disclosing gate-opening/closing events inside a flexible metal-organic framework loaded with CO2 by reactive and essential dynamics

We have combined reactive molecular dynamics simulations with principal component analysis to provide a clearer view of the interactions and motion of the CO2 molecules inside a metal-organic framework and the movements of the MOF components that regulate storage, adsorption, and diffusion of the guest species. The tens-of-nanometer size of the simulated model, the capability of the reactive force field tuned to reproduce the inorganic-organic material confidently, and the unconventional use of essential dynamics have effectively disclosed the gate-opening/closing phenomenon, possible coordinations of CO2 at the metal centers, all the diffusion steps inside the MOF channels, the primary motions of the linkers, and the effects of their concerted rearrangements on local CO2 relocations.

A DFT-metadynamics study disclosing key properties of ring-opening polymerization catalysts to produce polyethercarbonate polyols from cyclic ethylene carbonate as part of an emerging CCU technology

The ring opening polymerization of cyclic carbonates made from epoxide and CO₂ to CO₂-containing polymers constitutes an emerging technology of particular industrial interest. Considering the reaction of ring-opening polymerization of cyclic ethylene carbonate to produce polyethercarbonate polyols, several types of catalysts were tested experimentally and mechanistic pathways were proposed, but a detailed analysis of structure property relationship including the CO₂-liberation pathways is still lacking. This contribution is using computational methods to investigate reported benchmark catalysts with the lead structure A_xM_yO_z (A: alkali metal or alkyl, M: main group element or transition metal) that are particularly approved as effiecient catalysts for industrial purpose. Employing DFT-metadynamics simulations, free energy surfaces (FESs) for the key-steps in the catalytic polymerization of cyclic ethylene carbonate (cEC) are generated. Important structural criteria and characteristics of the catalysts that influence the catalytic performance and (side)reaction pathways are determined. It turns out that less nucleophilicity of the catalyst anion and more labile cations remain major criteria for prohibiting CO₂ liberation during polymerization. The key learnings of this contribution currently serve as a basis to develop the next generation of catalysts to bring this emerging carbon capture and use (CCU) technology into industrial application.

Reactions of Diethylazo-Dicarboxylate with Frustrated Lewis Pairs

Reactions of PAr3 /B(C6 F5 )3 (Ar=o-Tol, Mes, Ph) FLPs with diethyl azodicarboxylate (DEAD) afford the corresponding FLP addition products 1-3 in which P-N and B-O linkages are formed. In contrast, the reaction of BPh3 , PPh3 and DEAD gave product 4 where P-N and N-B linkages were confirmed. In all cases, other binding modes were computed to be both higher in energy and readily distinguishable by 31 P and 11 B NMR parameters. These data illustrate the influence of steric demands and electronic structures on the nature of the products of FLP reactions with DEAD.

Ab initio molecular dynamics with enhanced sampling in heterogeneous catalysis

Enhanced sampling ab initio simulations enable to study chemical phenomena in catalytic systems including thermal effects & anharmonicity, & collective dynamics describing enthalpic & entropic contributions, which can significantly impact on reaction free energy landscapes.

Diverse Uses of the Reaction of Frustrated Lewis Pair (FLP) with Hydrogen

The articulation of the notion of "frustrated Lewis pairs" (FLPs) emerged from the discovery that H₂ can be reversibly activated by combinations of sterically encumbered main group Lewis acids and bases. This has prompted numerous studies focused on various perturbations of the Lewis acid/base combinations and the applications to organic reductions. This Perspective focuses on the new directions and developments that are emerging from this FLP chemistry involving hydrogen. Three areas are discussed including new applications and approaches to FLP reductions, the reductions of small molecules, and the advances in heterogeneous FLP systems. These foci serve to illustrate that despite having its roots in main group chemistry, this simple concept of FLPs is being applied across the discipline.

Lewis Acidity of Carbon in Activated Carbonyl Group vs. B(C6 F5 )3 for Metal-Free Catalysis of Hydrogenation of Carbonyl Compounds

In this work, using DFT calculations, we investigated Lewis acidities of carbon (in activated carbonyl group) in comparison to the B(C₆ F₅ )₃ in combination with dioxane as the Lewis base (LB) for metal-free catalysis of heterolytic H₂ splitting and hydrogenation of carbonyl compounds. We found that in case of carbon as the Lewis acid (LA) the reaction is controlled by frontier molecular orbital interactions between the H₂ and LA-LB fragments at shorter distances. The steric effects can be reduced by electrophilic substitutions on the carbonyl carbon. Synergic combination between stronger orbital interactions and reduced steric effects can lower the barrier of the H₂ splitting below 10 kcal/mol. With the B(C₆ F₅ )₃ , the H₂ splitting is controlled by electrostatic interactions, which cause to form an early transition state. An advantage of employing Lewis acidity of the activated carbonyl carbon for hydrogenation is that the hydride-type attack and hydrogenation of the C=O bond occur in a single step throughout H₂ splitting. Hence, stronger Lewis acidity of the C(C=O) reinforces hydrogenation without prohibition of the hydride delivery.

Trends and Prospects in UiO-66 Metal-Organic Framework for CO2 Capture, Separation, and Conversion

Among thousands of known metal-organic frameworks (MOFs), the University of Oslo's MOF (UiO-66) exhibits unique structure topology, chemical and thermal stability, and intriguing tunable properties, that have gained incredible research interest. This paper summarizes the structural advancement of UiO-66 and its role in CO₂ capture, separation, and transformation into chemicals. The first part of the review summarizes the fast-growing literature related to the CO₂ capture reported by UiO-66 during the past ten years. The second part provides an overview of various advancements in UiO-66 membranes in CO₂ purification. The third part describes the role of UiO-66 and its composites as catalysts for CO₂ conversion into useful products. Despite many achievements, significant challenges associated with UiO-66 are addressed, and future perspectives are comprehensively presented to forecast how UiO-66 might be used further for CO₂ management.

Multiple catalytic sites in MOF-based hybrid catalysts for organic reactions

Hybrid catalysis provides an effective pathway to improve the catalytic efficiency and simplify the synthesis operation, but multiple catalytic sites are required. Catalysts with multiple functions based on/derived from metal-organic frameworks (MOFs) have received growing attention in organic synthesis due to their wide variety and outstanding designability. This review provides an overview of significant advances in the field of organic reactions by MOF-based hybrid catalysts with emphasis on multiple catalytic sites and their synergies, including inherent sites on host frameworks, sites of MOF composites and metal sites in/on MOF-derived hybrid catalysts.

Heterolytic Scission of Hydrogen Within a Crystalline Frustrated Lewis Pair

We report the heterolysis of molecular hydrogen under ambient conditions by the crystalline frustrated Lewis pair (FLP) 1-{2-[bis(pentafluorophenyl)boryl]phenyl}-2,2,6,6-tetramethylpiperidine (KCAT). The gas-solid reaction provides an approach to prepare the solvent-free, polycrystalline ion pair KCATH2 through a single crystal to single crystal transformation. The crystal lattice of KCATH2 increases in size relative to the parent KCAT by approximately 2%. Microscopy was used to follow the transformation of the highly colored red/orange KCAT to the colorless KCATH2 over a period of 2 h at 300 K under a flow of H₂ gas. There is no evidence of crystal decrepitation during hydrogen uptake. Inelastic neutron scattering employed over a temperature range from 4-200 K did not provide evidence for the formation of polarized H₂ in a precursor complex within the crystal at low temperatures and high pressures. However, at 300 K, the INS spectrum of KCAT transformed to the INS spectrum of KCATH2. Calculations suggest that the driving force is more favorable in the solid state compared to the solution or gas phase, but the addition of H₂ into the KCAT crystal is unfavorable. Ab Initio methods were used to calculate the INS spectra of KCAT, KCATH2, and a possible precursor complex of H₂ in the pocket between the B and N of crystalline KCAT. Ex-situ NMR showed that the transformation from KCAT to KCATH2 is quantitative and our results suggest that the hydrogen heterolysis process occurs via H₂ diffusion into the FLP crystal with a rate-limiting movement of H₂ from inactive positions to reactive sites.

Optimizing the Energetics of FLP-Type H2 Activation by Modulating the Electronic and Structural Properties of the Lewis Acids: A DFT Study

The great potential of frustrated Lewis pairs (FLPs) as metal-free catalysts for activation of molecular hydrogen has attracted increasing interest as an alternative to transition-metal catalysts. However, the complexity of FLP systems, involving the simultaneous interaction of three molecules, impedes a detailed understanding of the activation mechanism and the individual roles of the Lewis acid (LA) and Lewis base (LB). In the present work, using density functional theory (DFT) calculations, we examine the reactivity of 75 FLPs for the H₂ splitting reaction, including a series of experimentally investigated LAs combined with conventional phosphine-based (tBu₃P) and oxygen-based (i.e., ethereal solvent) Lewis bases. We find that the catalytic activity of the FLP is the result of a delicate balance of the LA and LB strengths and their bulkiness. The H₂ splitting reaction can be changed from endergonic to exergonic by tuning the electrophilicity of the LA. Also, a more nucleophilic LB results in a more stable ion pair product and a lower barrier for the hydrogen splitting. The bulkiness of the LB leads to an early transition state to reduce steric hindrance and lower the barrier height. The bulkiness of the fragments determines the cavity size in the FLP complex, and a large cavity allows for a larger charge separation in the ion pair configuration. A shorter proton-hydride distance in this product complex correlates with a stronger attraction between the fragments, which forms more reactive ion pairs and facilitates the proton and hydride donations in the subsequent hydrogenation process. These insights may help with rationalizing the experimentally observed reactivities of FLPs and with designing better FLP systems for hydrogenation catalysis and hydrogen storage.

Ab initio metadynamics calculations of dimethylamine for probing pKb variations in bulk vs. surface environments

Free energy landscape obtained from ab initio metadynamics calculations for dimethylamine protonation at the air–water interface.