Effect of Surface Oxidation on Oxidative Propane Dehydrogenation over Chromia: An Ab Initio Multiscale Kinetic Study

Machine learning predicts atomistic structures of multielement solid surfaces for heterogeneous catalysts in variable environments

Solid surfaces usually reach thermodynamic equilibrium through particle exchange with their environment under reactive conditions. A prerequisite for understanding their functionalities is detailed knowledge of the surface composition and atomistic geometry under working conditions. Owing to the large number of possible Miller indices and terminations involved in multielement solids, extensive sampling of the compositional and conformational space needed for reliable surface energy estimation is beyond the scope of ab initio calculations. Here, we demonstrate, using the case of iron carbides in environments with varied carbon chemical potentials, that the stable surface composition and geometry of multielement solids under reactive conditions, which involve large compositional and conformational spaces, can be predicted at ab initio accuracy using an approach that combines the bond valence model, Gaussian process regression, and ab initio thermodynamics. Determining the atomistic structure of surfaces under working conditions paves the way toward identifying the true active sites of multielement catalysts in heterogeneous catalysis.

Speeding up the Detection of Adsorbate Lateral Interactions in Graph-Theoretical Kinetic Monte Carlo Simulations

Kinetic Monte Carlo (KMC) has become an indispensable tool in heterogeneous catalyst discovery, but realistic simulations remain computationally demanding on account of the need to capture complex and long-range lateral interactions between adsorbates. The Zacros software package (https://zacros.org) adopts a graph-theoretical cluster expansion (CE) framework that allows such interactions to be computed with a high degree of generality and fidelity. This involves solving a series of subgraph isomorphism problems in order to identify relevant interaction patterns in the lattice. In an effort to reduce the computational burden, we have adapted two well-known subgraph isomorphism algorithms, namely, VF2 and RI, for use in KMC simulations and implemented them in Zacros. To benchmark their performance, we simulate a previously established model of catalytic NO oxidation, treating the O* lateral interactions with a series of progressively larger CEs. For CEs with long-range interactions, VF2 and RI are found to provide impressive speedups relative to simpler algorithms. RI performs best, giving speedups reaching more than 150× when combined with OpenMP parallelization. We also simulate a recently developed methane cracking model, showing that RI offers significant improvements in performance at high surface coverages.

Coupling acid catalysis and selective oxidation over MoO3-Fe2O3 for chemical looping oxidative dehydrogenation of propane

Redox catalysts play a vital role in chemical looping oxidative dehydrogenation processes, which have recently been considered to be a promising prospect for propylene production. This work describes the coupling of surface acid catalysis and selective oxidation from lattice oxygen over MoO₃-Fe₂O₃ redox catalysts for promoted propylene production. Atomically dispersed Mo species over γ-Fe₂O₃ introduce effective acid sites for the promotion of propane conversion. In addition, Mo could also regulate the lattice oxygen activity, which makes the oxygen species from the reduction of γ-Fe₂O₃ to Fe₃O₄ contribute to selectively oxidative dehydrogenation instead of over-oxidation in pristine γ-Fe₂O₃. The enhanced surface acidity, coupled with proper lattice oxygen activity, leads to a higher surface reaction rate and moderate oxygen diffusion rate. Consequently, this coupling strategy achieves a robust performance with 49% of propane conversion and 90% of propylene selectivity for at least 300 redox cycles and ultimately demonstrates a potential design strategy for more advanced redox catalysts.

Correlation between the Properties of Surface Lattice Oxygen on NiO and Its Reactivity and Selectivity towards the Oxidative Dehydrogenation of Propane

Modified NiO catalysts with controllable vacancies and dopants are promising for alkene production from oxidative dehydrogenation (ODH) of light alkanes, and a molecular understanding of the modification on elementary reaction steps would facilitate the design of highly efficient catalysts and catalytic processes. In this study, density functional theory (DFT) calculations was used to map out the complete reaction pathways of propane ODH on the NiO (100) surfaces with different modifiers. The results demonstrated that the presence of vacancies (O and Ni) and dopants (Li and Al) alters the electrophilicity of surface oxygen species, which in turn affects the reactivity towards C-H bond activation and the overall catalytic activity and selectivity. The strongly electrophilic O favors a radical mechanism for the first C-H activation on O followed by the second C-H activation on O-O site, whereas weak electrophilic O favors concerted C-H bond breaking over Ni-O site. The C-H bond activation proceeds through a late transition state, characterized by the almost completion of the O-H bond formation. Consequently, the adsorption energy of H adatom on O rather than p-band center or Bader charge of O has been identified to be an accurate descriptor to predict the activation barrier for C-H breaking (activity) as well as the difference between the activation barriers of propene and CH₃ CCH₃ (selectivity) of ODH.

Quantitative Evidence to Challenge the Traditional Model in Heterogeneous Catalysis: Kinetic Modeling for Ethane Dehydrogenation over Fe/SAPO-34

The production of ethylene from ethane dehydrogenation (EDH) is of great importance in the chemical industry, where zeolites are reported to be promising catalysts and kinetic simulations using the energetics from quantum mechanical calculations might provide an effective approach to speed up the development. However, the kinetic simulations with rigorous considerations of the zeolite environment are not yet advanced. In this work, EDH over Fe/SAPO-34 is investigated using quantum mechanical calculations with kinetic simulations. We show that an excellent agreement between the reaction rates from the self-consistent kinetic simulations using the coverage-dependent kinetic model developed in this work and the experimental ones can be achieved. We demonstrate that the adsorbate-adsorbate interactions are of paramount importance to the accuracy of kinetic calculations for zeolite catalysts. Our self-consistent kinetic calculations illustrate that the CH₃CH₂• radical rather than CH₃CH₂* is a favored intermediate. Perhaps more importantly, we reveal that the traditional model to describe catalytic reactions in heterogeneous catalysis cannot be used for the kinetics of the system and it may not be appropriate for many real catalytic systems. This work not only builds a framework for accurate kinetic simulations in zeolites, but also emphasizes an important concept beyond the traditional model.

In Situ Generated Boron Peroxo as Mild Oxidant in Propane Oxidative Dehydrogenation Revealed by Density Functional Theory Study

Boron-based materials catalyzing oxidative dehydrogenation is emerging as a promising protocol for efficient conversion of light alkanes to olefins, while the origin of its remarkable selectivity remains unclear. By means of density functional theory calculations, this work addresses the crucial role of boron peroxo as the mild oxidant in propane ODH: (1) Surface boron peroxo species can be generated in situ in the presence of peroxo species, preferably at the >B-O-B< sites of the zigzag edge, and show high activity to dehydrogenate propane (ΔG^⧧ = 13.5 kcal/mol, ΔG = 8.9 kcal/mol). (2) The >B-O-O· site shows high discriminability of secondary H over primary H of the propane molecule, leading to significantly higher yield of iso-propyl (CH₃ĊHCH₃) than n-propyl (CH₃CH₂ĊH₂); thus, propene formation is favored over deep oxidation. This provides physical insights into the origin of the remarkable olefin selectivity in the boron-containing ODH catalytic systems.

Light olefin synthesis from a diversity of renewable and fossil feedstocks: state-of the-art and outlook

Light olefins are important feedstocks and platform molecules for the chemical industry. Their synthesis has been a research priority in both academia and industry. There are many different approaches to the synthesis of these compounds, which differ by the choice of raw materials, catalysts and reaction conditions. The goals of this review are to highlight the most recent trends in light olefin synthesis and to perform a comparative analysis of different synthetic routes using several quantitative characteristics: selectivity, productivity, severity of operating conditions, stability, technological maturity and sustainability. Traditionally, on an industrial scale, the cracking of oil fractions has been used to produce light olefins. Methanol-to-olefins, alkane direct or oxidative dehydrogenation technologies have great potential in the short term and have already reached scientific and technological maturities. Major progress should be made in the field of methanol-mediated CO and CO₂ direct hydrogenation to light olefins. The electrocatalytic reduction of CO₂ to light olefins is a very attractive process in the long run due to the low reaction temperature and possible use of sustainable electricity. The application of modern concepts such as electricity-driven process intensification, looping, CO₂ management and nanoscale catalyst design should lead in the near future to more environmentally friendly, energy efficient and selective large-scale technologies for light olefin synthesis.

A critical evaluation of the catalytic role of CO2 in propane dehydrogenation catalyzed by chromium oxide from a DFT-based microkinetic simulation

Propane dehydrogenation under CO₂ is an important catalytic route to obtain propene with a good balance between selectivity and stability. However, a precise description of the catalytic role of CO₂ in propane dehydrogenation is still absent. In this work, we focus on the elucidation of the role of CO₂ by using DFT-based microkinetic simulation. The influence of CO₂ is categorized as direct and indirect effects. It was found that the chemisorbed CO₂ can directly abstract hydrogen from propane and propyl with a comparable barrier to the counterpart at the surface oxygen site. On the other hand, the dissociation of CO₂ yields active surface species of CO* and O* which are actively involved in the removal of surface hydroxyls. It is found that the TOFs of both propane conversion and propene formation are significantly increased with the presence of CO₂, which is explained by the reduced apparent activation energy. The primary hydrogen abstraction is identified to be the most influential step from the DRC analysis. The main effects of CO₂ are concluded to be removing hydrogen and restoring oxygen vacancies from reaction pathway analysis.

Kinetic Monte Carlo simulations for heterogeneous catalysis: Fundamentals, current status, and challenges

Kinetic Monte Carlo (KMC) simulations in combination with first-principles (1p)-based calculations are rapidly becoming the gold-standard computational framework for bridging the gap between the wide range of length scales and time scales over which heterogeneous catalysis unfolds. 1p-KMC simulations provide accurate insights into reactions over surfaces, a vital step toward the rational design of novel catalysts. In this Perspective, we briefly outline basic principles, computational challenges, successful applications, as well as future directions and opportunities of this promising and ever more popular kinetic modeling approach.

A theoretical study on the mechanism of conversion of C3H8 and CO2 to C3H6 and HCOOH by M4–B24N28 catalysis

In this work, we explored the reaction mechanism of metallic copper-doped modified aperiodic (BN)28 nanocages for the catalytic oxidation of propane with carbon dioxide. The DFT calculation is performed on all possible paths during the reaction process.