Computational Screening of Bimetal-Functionalized Zr6O8 MOF Nodes for Methane C–H Bond Activation

Theoretical study of fructose adsorption and conversion to trioses on metal-organic frameworks

The conversion of chemically modified biomass into more valuable chemicals has recently gained significant attention from industry. In this study, we investigate the adsorption of fructose and its conversion into two trioses, glyceraldehyde (GLA) and dihydroxyacetone (DHA), on metal-organic frameworks using density functional theory calculations. The reaction mechanism proceeds through two main steps: first, the opening of the fructose ring; second, the retro-aldol fragmentation, which is favored over intramolecular hydrogen shifts. The substitution of a tetravalent metal in the metal-organic framework leads to different adsorption strengths in the order Hf-NU-1000 > Zr-NU-1000 > Ti-NU-1000. The catalytic activities of Hf-NU-1000 and Zr-NU-1000 are found to be similar. Both are more active than Ti-NU1000, corresponding to their relative Lewis acidity. It was found that functionalization of the organic linkers of the Hf-NU-1000 MOF does not improve its catalytic activity. The catalytic activity follows the order Hf-MOF-808 > Hf-NU-1000 > Hf-UIO-66 when based on either the overall activation energy or the turnover frequency (TOF).

Design and Advanced Manufacturing of NU-1000 Metal-Organic Frameworks with Future Perspectives for Environmental and Renewable Energy Applications

Metal-organic frameworks (MOFs) represent a relatively new family of materials that attract lots of attention thanks to their unique features such as hierarchical porosity, active metal centers, versatility of linkers/metal nodes, and large surface area. Among the extended list of MOFs, Zr-based-MOFs demonstrate comparably superior chemical and thermal stabilities, making them ideal candidates for energy and environmental applications. As a Zr-MOF, NU-1000 is first synthesized at Northwestern University. A comprehensive review of various approaches to the synthesis of NU-1000 MOFs for obtaining unique surface properties (e.g., diverse surface morphologies, large surface area, and particular pore size distribution) and their applications in the catalysis (electro-, and photo-catalysis), CO2 reduction, batteries, hydrogen storage, gas storage/separation, and other environmental fields are presented. The review further outlines the current challenges in the development of NU-1000 MOFs and their derivatives in practical applications, revealing areas for future investigation.

Screening Cu-Zeolites for Methane Activation Using Curriculum-Based Training

Machine learning (ML), when used synergistically with atomistic simulations, has recently emerged as a powerful tool for accelerated catalyst discovery. However, the application of these techniques has been limited by the lack of interpretable and transferable ML models. In this work, we propose a curriculum-based training (CBT) philosophy to systematically develop reactive machine learning potentials (rMLPs) for high-throughput screening of zeolite catalysts. Our CBT approach combines several different types of calculations to gradually teach the ML model about the relevant regions of the reactive potential energy surface. The resulting rMLPs are accurate, transferable, and interpretable. We further demonstrate the effectiveness of this approach by exhaustively screening thousands of [CuOCu]2+ sites across hundreds of Cu-zeolites for the industrially relevant methane activation reaction. Specifically, this large-scale analysis of the entire International Zeolite Association (IZA) database identifies a set of previously unexplored zeolites (i.e., MEI, ATN, EWO, and CAS) that show the highest ensemble-averaged rates for [CuOCu]2+-catalyzed methane activation. We believe that this CBT philosophy can be generally applied to other zeolite-catalyzed reactions and, subsequently, to other types of heterogeneous catalysts. Thus, this represents an important step toward overcoming the long-standing barriers within the computational heterogeneous catalysis community.

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.

Full Spectroscopic Characterization of the Molecular Oxygen-Based Methane to Methanol Conversion over Open Fe(II) Sites in a Metal-Organic Framework

Iron-based enzymes efficiently activate molecular oxygen to perform the oxidation of methane to methanol (MTM), a reaction central to the contemporary chemical industry. Conversely, a very limited number of artificial catalysts have been devised to mimic this process. Herein, we employ the MIL-100(Fe) metal-organic framework (MOF), a material that exhibits isolated Fe sites, to accomplish the MTM conversion using O2 as the oxidant under mild conditions. We apply a diverse set of advanced operando X-ray techniques to unveil how MIL-100(Fe) can act as a catalyst for direct MTM conversion. Single-phase crystallinity and stability of the MOF under reaction conditions (200 or 100 °C, CH4 + O2) are confirmed by X-ray diffraction measurements. X-ray absorption, emission, and resonant inelastic scattering measurements show that thermal treatment above 200 °C generates Fe(II) sites that interact with O2 and CH4 to produce methanol. Experimental evidence-driven density functional theory (DFT) calculations illustrate that the MTM reaction involves the oxidation of the Fe(II) sites to Fe(III) via a high-spin Fe(IV)═O intermediate. Catalyst deactivation is proposed to be caused by the escape of CH3• radicals from the relatively large MOF pore cages, ultimately resulting in the formation of hydroxylated triiron units, as proven by valence-to-core X-ray emission spectroscopy. The O2-based MTM catalytic activity of MIL-100(Fe) in the investigated conditions is demonstrated for two consecutive reaction cycles, proving the MOF potential toward active site regeneration. These findings will desirably lay the groundwork for the design of improved MOF catalysts for the MTM conversion.

Computational Screening of Supported Metal Oxide Nanoclusters for Methane Activation: Insights into Homolytic versus Heterolytic C-H Bond Dissociation

Since its discovery in zeolites, the [CuOCu]²⁺ motif has played an important role in our understanding of selective methane activation over supported metal oxide nanoclusters. Although there are two known C-H bond dissociation mechanisms, namely, homolytic and heterolytic cleavage, most computational studies on optimizing metal oxide nanoclusters for improved methane activation reactivity have focused only on the homolytic mechanism. In this work, both mechanisms were examined for a set of 21 mixed metal oxide complexes of the form of [M₁OM₂]²⁺ (M₁ and M₂ = Mn, Fe, Co, Ni, Cu, and Zn). Except for pure copper, heterolytic cleavage was found to be the dominant C-H bond activation pathway for all systems. Furthermore, mixed systems including [CuOMn]²⁺, [CuONi]²⁺, and [CuOZn]²⁺ are predicted to possess methane activation activity similar to pure [CuOCu]²⁺. These results suggest that both homolytic and heterolytic mechanisms should be considered in computing methane activation energies on supported metal oxide nanoclusters.

Modulating the Energetics of C-H Bond Activation in Methane by Utilizing Metalated Porphyrinic Metal-Organic Frameworks

In recent years, much effort has been directed toward utilizing metal-organic frameworks (MOFs) for activating C-H bonds of light alkanes. The energy demanding steps involved in the catalytic pathway are the formation of metal-oxo species and the subsequent cleavage of the C-H bonds of alkanes. With the intention of exploring the tunability of the activation barriers involved in the catalytic pathway of methane hydroxylation, we have employed density functional theory to model metalated porphyrinic MOFs (MOF-525(M)). We find that the heavier congeners down a particular group have high exothermic oxo-formation enthalpies ΔH_O and hence are associated with low N₂O activation barriers. Independent analyses of activation barriers and structure-activity relationship leads to the conclusion that MOF-525(Ru) and MOF-525(Ir) can act as an effective catalysts for methane hydroxylation. Hence, ΔH_O has been found to act as a guide, in the first place, in choosing the optimum catalyst for methane hydroxylation from a large set of available systems.

Partial Oxidation of Methane to Methanol on the M-O-Ag/Graphene (M = Ag, Cu) Composite Catalyst: A DFT Study

Partial oxidation of methane (CH₄) to methanol (CH₃OH) remains a great challenge in the field of catalysis due to its low selectivity and productivity. Herein, Ag-O-Ag/graphene and Cu-O-Ag/graphene composite catalysts are proposed to oxidize methane (CH₄) to methanol (CH₃OH) by using the first-principles calculations. It is shown that reactive oxygen species (μ-O) on both catalysts can activate the C-H bond of CH₄, and in addition to CH₄ activation, the catalytic activity follows the order of Ag-O-Ag/graphene (singlet) > Ag-O-Ag/graphene (triplet) ≈ Cu-O-Ag/graphene (triplet) > Cu-O-Ag/graphene (singlet). For CH₃OH* formation, the catalytic activity follows the order of Cu-O-Ag/graphene (triplet) > Ag-O-Ag/graphene (triplet) > Ag-O-Ag/graphene (singlet) > Cu-O-Ag/graphene (singlet). It can be inferred that the introduction of Cu not only reduces the use of noble metal Ag but also exhibits a catalytic effect comparable to that of the Ag-O-Ag/graphene catalyst. Our findings will provide a new avenue for understanding and designing highly effective catalysts for the direct conversion of CH₄ to CH₃OH.

Electrostatic potential-derived charge: a universal OER performance descriptor for MOFs

Metal-organic frameworks (MOFs) provide opportunities for the design of high-efficiency catalysts attributed to their high compositional and structural tunability. Meanwhile, the huge number of MOFs poses a great challenge to experimental-intensive development of high-performance functional applications. By taking the computationally feasible and structurally representative trigonal prismatic secondary building units (SBUs) of MOFs as the entry point, we introduce a descriptor-based approach for designing high-performance MOFs for the oxygen evolution reaction (OER). The electrostatic potential-derived charge (ESPC) is identified as a robust and universal OER performance descriptor of MOFs, showing a distinct linear relationship with the onset potentials of OER elemental steps. Importantly, we establish an ESPC-based physical pattern of active site-intermediate binding strength, which interprets the rationality of ESPC as an OER performance descriptor. We further reveal that the SBUs with Ni/Cu as active site atoms while Mn/Fe/Co/Ni as spectator atoms have excellent OER activity through the variation pattern of ESPC along with metal composition. The universal correlation between ESPC and OER activity provides a rational rule for designing high-performance MOF-based OER electrocatalysts and can be easily extended to design functional MOFs for a rich variety of catalytic applications.

Exploring mechanistic routes for light alkane oxidation with an iron-triazolate metal-organic framework

In this work, we computationally explore the formation and subsequent reactivity of various iron-oxo species in the iron-triazolate framework Fe₂(μ-OH)₂(bbta) (H₂bbta = 1H,5H-benzo(1,2-d:4,5-d')bistriazole) for the catalytic activation of strong C-H bonds. With the direct conversion of methane to methanol as the probe reaction of interest, we use density functional theory (DFT) calculations to evaluate multiple mechanistic pathways in the presence of either N₂O or H₂O₂ oxidants. These calculations reveal that a wide range of transition metal-oxo sites - both terminal and bridging - are plausible in this family of metal-organic frameworks, making it a unique platform for comparing the electronic structure and reactivity of different proposed active site motifs. Based on the DFT calculations, we predict that Fe₂(μ-OH)₂(bbta) would exhibit a relatively low barrier for N₂O activation and energetically favorable formation of an [Fe(O)]²⁺ species that is capable of oxidizing C-H bonds. In contrast, the use of H₂O₂ as the oxidant is predicted to yield an assortment of bridging iron-oxo sites that are less reactive. We also find that abstracting oxo ligands can exhibit a complex mixture of both positive and negative spin density, which may have broader implications for relating the degree of radical character to catalytic activity. In general, we consider the coordinatively unsaturated iron sites to be promising for oxidation catalysis, and we provide several recommendations on how to further tune the catalytic properties of this family of metal-triazolate frameworks.

Light alkane oxidation over well-defined active sites in metal–organic framework materials

We review structure–catalytic property relationships for MOF materials used in the direct oxidation of light alkanes, focusing specifically on the elucidation of active site structures and probes for reaction mechanisms.

Density Functional Investigation of the Conversion of Furfural to Furfuryl Alcohol by Reaction with i-Propanol over UiO-66 Metal-Organic Framework

Carbonyl C═O bond reduction via catalytic transfer hydrogenation (CTH) is one of the essential processes for biomass conversion to valuable chemicals and fuels. Here, we investigate the CTH of furfural to furfuryl alcohol with i-propanol on UiO-66 metal-organic frameworks using density functional theory calculations and linear scaling relations. Initially, the reaction over two defect sites presented on Zr-UiO-66, namely, dehydrated and hydrated sites, have been compared. The hydrated active site is favored over that on the dehydrated active site since the activation free energy of the rate-determining reaction step occurring on the hydrated active site is lower than that occurring on the dehydrated active site (14.9 vs 17.9 kcal/mol). The catalytic effect of exchanged tetravalent metals (Hf and Ti) on Zr-UiO-66 is also considered. We found that Hf-UiO-66 (13.5 kcal/mol) provides a lower activation energy than Zr-UiO-66 (14.9 kcal/mol) and Ti-UiO-66 (19.4 kcal/mol), which corresponds to it having a higher Lewis acidity. The organic linkers of UiO-66 MOFs play a role in stabilizing all of the species on potential energy surfaces. The linear scaling relationship also reveals the significant role of the UiO-66 active site in activating the carbonyl C═O of furfural, and strong relationships are observed between the activation free energy, the charge of the metal at the MOF active sites, and the complexation energies in reaction coordinates.

Pyrene-based metal organic frameworks: from synthesis to applications

Pyrene is one of the most widely investigated aromatic hydrocarbons given to its unique optical and electronic properties. Hence, pyrene-based ligands have been attractive for the synthesis of metal-organic frameworks (MOFs) in the last few years. In this review, we will focus on the most important characteristics of pyrene, in addition to the development and synthesis of pyrene-based molecules as bridging ligands to be used in MOF structures. We will summarize the synthesis attempts, as well as the post-synthetic modifications of pyrene-based MOFs by the incorporation of metals or ligands in the structure. The discussion of promising results of such MOFs in several applications; including luminescence, photocatalysis, adsorption and separation, heterogeneous catalysis, electrochemical applications and bio-medical applications will be highlighted. Finally, some insights and future prospects will be given based on the studies discussed in the review. This review will pave the way for the researchers in the field for the design and development of novel pyrene-based structures and their utilization for different applications.

Insights into Catalytic Gas-Phase Hydrolysis of Organophosphate Chemical Warfare Agents by MOF-Supported Bimetallic Metal-Oxo Clusters

Zirconium-based metal-organic frameworks (Zr-MOFs) have been reported to be efficient catalysts for the hydrolysis of organophosphate chemical warfare agents (CWAs) in buffered solutions. However, for the gas-phase reaction, which is more relevant to the situation in a battlefield gas mask application, the kinetics of Zr-MOF catalysts may be severely hindered by strong product inhibition. To improve the catalytic performance, we computationally screened a series of synthetically accessible Zr-MOF-supported bimetallic metal-oxo clusters in which the metal-oxygen-metal active motif is preserved, aiming to find catalysts that have lower binding affinities to the hydrolysis product. For the promising catalyst Al₂O₂(OH)₂@NU-1000 identified from the screening using density functional theory, we mapped out the full reaction pathway of gas-phase dimethyl p-nitrophenolphosphate (DMNP) hydrolysis and analyzed the free energy profile as well as the turnover frequency (TOF). We found that the catalytic mechanism on the new catalyst is slightly different from the one on NU-1000, which also led to a different TOF-limiting step. Additional factors that can affect the overall catalytic performance in practical application, such as the amount of ambient moisture and the existence of acid gases that may poison the catalyst, have also been evaluated.

DFT Study on the Catalytic Activity of ALD-Grown Diiron Oxide Nanoclusters for Partial Oxidation of Methane to Methanol

Using density functional theory (DFT), we studied the catalytic activity of iron oxide nanoclusters that mimic the structure of the active site in the soluble form of methane monooxygenase (sMMO) for the partial oxidation of methane to methanol. Using N₂O as the oxidant, we consider a radical-rebound mechanism and a concerted mechanism for the oxidation of methane on either a bridging oxygen (O_b) or a terminal oxygen (O_t) active site. We find that the radical-rebound pathway is preferred over the concerted pathway by 40-50 kJ/mol, but the desorption of methanol and the regeneration of the oxygen site are found to be the highest barriers for the direct conversion of methane to methanol with these catalysts. As demonstrated by a population analysis, the O_x (x = b or t) site behaves as an oxygen radical during the H abstraction, and the [Fe⁺-O_x^-] site behaves as a Lewis acid-base pair during the concerted C-H cleavage. Molecular orbital decomposition analysis further demonstrates electron transfer during the oxidation and reduction steps of the reaction. High-level multireference calculations were also carried out to further assess the DFT results. Understanding how these systems behave during the proposed reaction pathways provides new insights into how they can be tuned for methane partial oxidation.

Selective Methane Oxidation to Methanol on Cu-Oxo Dimers Stabilized by Zirconia Nodes of an NU-1000 Metal-Organic Framework

Mononuclear and dinuclear copper species were synthesized at the nodes of an NU-1000 metal-organic framework (MOF) via cation exchange and subsequent oxidation at 200 °C in oxygen. Copper-exchanged MOFs are active for selectively converting methane to methanol at 150-200 °C. At 150 °C and 1 bar methane, approximately a third of the copper centers are involved in converting methane to methanol. Methanol productivity increased by 3-4-fold and selectivity increased from 70% to 90% by increasing the methane pressure from 1 to 40 bar. Density functional theory showed that reaction pathways on various copper sites are able to convert methane to methanol, the copper oxyl sites with much lower free energies of activation. Combining studies of the stoichiometric activity with characterization by in situ X-ray absorption spectroscopy and density functional theory, we conclude that dehydrated dinuclear copper oxyl sites formed after activation at 200 °C are responsible for the activity.

Computational Approach to Molecular Catalysis by 3d Transition Metals: Challenges and Opportunities

Computational chemistry provides a versatile toolbox for studying mechanistic details of catalytic reactions and holds promise to deliver practical strategies to enable the rational in silico catalyst design. The versatile reactivity and nontrivial electronic structure effects, common for systems based on 3d transition metals, introduce additional complexity that may represent a particular challenge to the standard computational strategies. In this review, we discuss the challenges and capabilities of modern electronic structure methods for studying the reaction mechanisms promoted by 3d transition metal molecular catalysts. Particular focus will be placed on the ways of addressing the multiconfigurational problem in electronic structure calculations and the role of expert bias in the practical utilization of the available methods. The development of density functionals designed to address transition metals is also discussed. Special emphasis is placed on the methods that account for solvation effects and the multicomponent nature of practical catalytic systems. This is followed by an overview of recent computational studies addressing the mechanistic complexity of catalytic processes by molecular catalysts based on 3d metals. Cases that involve noninnocent ligands, multicomponent reaction systems, metal-ligand and metal-metal cooperativity, as well as modeling complex catalytic systems such as metal-organic frameworks are presented. Conventionally, computational studies on catalytic mechanisms are heavily dependent on the chemical intuition and expert input of the researcher. Recent developments in advanced automated methods for reaction path analysis hold promise for eliminating such human-bias from computational catalysis studies. A brief overview of these approaches is presented in the final section of the review. The paper is closed with general concluding remarks.

Computational study of the carbonyl-ene reaction between formaldehyde and propylene encapsulated in coordinatively unsaturated metal-organic frameworks M3(btc)2 (M = Fe, Co, Ni, Cu and Zn)

The carbonyl-ene reaction between encapsulated formaldehyde and propylene over the coordinatively unsaturated metal-organic frameworks M3(btc)2 (M = Fe, Co, Ni, Cu and Zn) has been investigated by means of density functional calculations. Zn3(btc)2 adsorbs formaldehyde strongest due to electron delocalization between Zn and the oxygen atom of the reactant molecule. The reaction is proposed to proceed in a single step involving proton transfer and carbon-carbon bond formation. We find the relative catalytic activity to be Zn3(btc)2 > Fe3(btc)2 ≥ Co3(btc)2 > Ni3(btc)2 > Cu3(btc)2, based on activation energy and turnover frequencies (TOF). The low activation energy for Zn3(btc)2 compared to the others can be explained by the delocalization of electron density between the carbonyl bond and the catalyst active sites, leading to a more stable transition state. The five MOFs are used to propose a descriptor for the relationship between activation energy on one side and electronic properties or adsorption energies on the other side in order to allow a quick screening of other catalytic materials for this reaction.

Mixed-metal metal–organic frameworks

Mixed-metal MOFs contain at least 2 different metal ions presenting promising potential in heterogeneous catalysis, gas sorption/separation, luminescence and sensing.

Metal-Organic Framework (MOF)-Based Materials as Heterogeneous Catalysts for C-H Bond Activation

Converting light hydrocarbons such as methane, ethane, propane, and cyclohexane into value-added chemicals and fuel products by means of direct C-H functionalization is an attractive method in the petrochemical industry. As they emerge as a relatively new class of porous solid materials, metal-organic frameworks (MOFs) are appealing as single-site heterogeneous catalysts or catalytic supports for C-H bond activation. In contrast to the traditional microporous and mesoporous materials, MOFs feature high porosity, functional tunability, and molecular-level characterization for the study of structure-property relationships. These virtues make MOFs ideal platforms to develop catalysts for C-H activation with high catalytic activity, selectivity, and recyclability under relatively mild reaction conditions. This review highlights the research aimed at the implementation of MOFs as single-site heterogeneous catalysts for C-H bond activation. It provides insight into the rational design and synthesis of three types of stable MOF catalysts for C-H bond activation, that is, i) metal nodes as catalytic sites, ii) the incorporation of catalytic sites into organic struts, and iii) the incorporation of catalytically active guest species into pores of MOFs. Here, the rational design and synthesis of MOF catalysts that lead to the distinct catalytic property for C-H bond activation are discussed along with the post-synthesis of MOFs, intriguing functions with MOF catalysts, and microenvironments that lead to the distinct catalytic properties of MOF catalysts.

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.

Computational Design of Functionalized Metal-Organic Framework Nodes for Catalysis

Recent progress in the synthesis and characterization of metal-organic frameworks (MOFs) has opened the door to an increasing number of possible catalytic applications. The great versatility of MOFs creates a large chemical space, whose thorough experimental examination becomes practically impossible. Therefore, computational modeling is a key tool to support, rationalize, and guide experimental efforts. In this outlook we survey the main methodologies employed to model MOFs for catalysis, and we review selected recent studies on the functionalization of their nodes. We pay special attention to catalytic applications involving natural gas conversion.

Theoretical study of CO2 hydrogenation into formic acid on Lewis acid zeolites

The hydrogenation of carbon dioxide (CO₂) to formic acid over Lewis acidic zeolites as catalyst has been investigated by means of density functional theory (DFT) with the M06-L functional.