Mechanistic Characterization of Zeolite-Catalyzed Aromatic Electrophilic Substitution at Realistic Operating Conditions

Use of In Situ X-ray Absorption to Probe Reactivity: A Catalysis Golden Rule

The decomposition of ozone on supported manganese oxide catalysts, studied here, exemplifies reactions involving electron transfer. In situ extended X-ray absorption fine-structure spectra (Mn K-edge) on in situ treated samples show that the supported phase in MnOx/SiO2 resembles Mn3O4 while that in MnOx/Al2O3 samples resembles MnO2. In situ Raman spectroscopy shows the involvement of a common peroxide surface species. Kinetic data indicate a nonuniform surface and a rate-determining step (rds) involving electron transfer from the peroxide intermediate. The activation energy for all the catalysts is the same, indicating that the pre-exponential factor controls the rate. This can be associated with the electronic partition function (unoccupied density of states (DOS)), and X-ray absorption near-edge spectroscopy duly shows that the areas of the pre-edge peaks due to 1s to 3d transitions track the reactivity trends. The relation to the unoccupied DOS is analogous to Fermi's Golden Rule for electronic transitions, and we denote the finding here, applicable to reactions involving electron transfer, as a Catalysis Golden Rule.

Unraveling Spatially Dependent Hydrophilicity and Reactivity of Confined Carbocation Intermediates during Methanol Conversion over ZSM-5 Zeolite

Carbocations play a pivotal role as reactive intermediates in zeolite-catalyzed methanol-to-hydrocarbon (MTH) transformations. However, the interaction between carbocations and water vapor and its subsequent effects on catalytic performance remain poorly understood. Using micro-magnetic resonance imaging (μMRI) and solid-state NMR techniques, this work investigates the hydrophilic behavior of cyclopentenyl cations within ZSM-5 pores under vapor conditions. We show that the polar cationic center of cyclopentenyl cations readily initiates water nucleus formation through water molecule capture. This leads to an inhomogeneous water adsorption gradient along the axial positions of zeolite, correlating with the spatial distribution of carbocation concentrations. The adsorbed water promotes deprotonation and aromatization of cyclopentenyl cations, significantly enhancing the aromatic product selectivity in MTH catalysis. These results reveal the important influence of adsorbed water in modulating the carbocation reactivity within confined zeolite pores.

Water structures on acidic zeolites and their roles in catalysis

The local reaction environment of catalytic active sites can be manipulated to modify the kinetics and thermodynamic properties of heterogeneous catalysis. Because of the unique physical-chemical nature of water, heterogeneously catalyzed reactions involving specific interactions between water molecules and active sites on catalysts exhibit distinct outcomes that are different from those performed in the absence of water. Zeolitic materials are being applied with the presence of water for heterogeneous catalytic reactions in the chemical industry and our transition to sustainable energy. Mechanistic investigation and in-depth understanding about the behaviors and the roles of water are essentially required for zeolite chemistry and catalysis. In this review, we focus on the discussions of the nature and structures of water adsorbed/stabilized on Brønsted and Lewis acidic zeolites based on experimental observations as well as theoretical calculation results. The unveiled functions of water structures in determining the catalytic efficacy of zeolite-catalyzed reactions have been overviewed and the strategies frequently developed for enhancing the stabilization of zeolite catalysts are highlighted. Recent advancement will contribute to the development of innovative catalytic reactions and the rationalization of catalytic performances in terms of activity, selectivity and stability with the presence of water vapor or in condensed aqueous phase.

Reference-Quality Free Energy Barriers in Catalysis from Machine Learning Thermodynamic Perturbation Theory

For the first time, we report calculations of the free energies of activation of cracking and isomerization reactions of alkenes that combine several different electronic structure methods with molecular dynamics simulations. We demonstrate that the use of a high level of theory (here Random Phase Approximation-RPA) is necessary to bridge the gap between experimental and computed values. These transformations, catalyzed by zeolites and proceeding via cationic intermediates and transition states, are building blocks of many chemical transformations for valorization of long chain paraffins originating, e.g., from plastic waste, vegetable oils, Fischer-Tropsch waxes or crude oils. Compared with the free energy barriers computed at the PBE+D2 production level of theory via constrained ab initio molecular dynamics, the barriers computed at the RPA level by the application of Machine Learning thermodynamic Perturbation Theory (MLPT) show a significant decrease for isomerization reaction and an increase of a similar magnitude for cracking, yielding an unprecedented agreement with the results obtained by experiments and kinetic modeling.

Operando Modeling of Zeolite-Catalyzed Reactions Using First-Principles Molecular Dynamics Simulations

Within this Perspective, we critically reflect on the role of first-principles molecular dynamics (MD) simulations in unraveling the catalytic function within zeolites under operating conditions. First-principles MD simulations refer to methods where the dynamics of the nuclei is followed in time by integrating the Newtonian equations of motion on a potential energy surface that is determined by solving the quantum-mechanical many-body problem for the electrons. Catalytic solids used in industrial applications show an intriguing high degree of complexity, with phenomena taking place at a broad range of length and time scales. Additionally, the state and function of a catalyst critically depend on the operating conditions, such as temperature, moisture, presence of water, etc. Herein we show by means of a series of exemplary cases how first-principles MD simulations are instrumental to unravel the catalyst complexity at the molecular scale. Examples show how the nature of reactive species at higher catalytic temperatures may drastically change compared to species at lower temperatures and how the nature of active sites may dynamically change upon exposure to water. To simulate rare events, first-principles MD simulations need to be used in combination with enhanced sampling techniques to efficiently sample low-probability regions of phase space. Using these techniques, it is shown how competitive pathways at operating conditions can be discovered and how broad transition state regions can be explored. Interestingly, such simulations can also be used to study hindered diffusion under operating conditions. The cases shown clearly illustrate how first-principles MD simulations reveal insights into the catalytic function at operating conditions, which could not be discovered using static or local approaches where only a few points are considered on the potential energy surface (PES). Despite these advantages, some major hurdles still exist to fully integrate first-principles MD methods in a standard computational catalytic workflow or to use the output of MD simulations as input for multiple length/time scale methods that aim to bridge to the reactor scale. First of all, methods are needed that allow us to evaluate the interatomic forces with quantum-mechanical accuracy, albeit at a much lower computational cost compared to currently used density functional theory (DFT) methods. The use of DFT limits the currently attainable length/time scales to hundreds of picoseconds and a few nanometers, which are much smaller than realistic catalyst particle dimensions and time scales encountered in the catalysis process. One solution could be to construct machine learning potentials (MLPs), where a numerical potential is derived from underlying quantum-mechanical data, which could be used in subsequent MD simulations. As such, much longer length and time scales could be reached; however, quite some research is still necessary to construct MLPs for the complex systems encountered in industrially used catalysts. Second, most currently used enhanced sampling techniques in catalysis make use of collective variables (CVs), which are mostly determined based on chemical intuition. To explore complex reactive networks with MD simulations, methods are needed that allow the automatic discovery of CVs or methods that do not rely on a priori definition of CVs. Recently, various data-driven methods have been proposed, which could be explored for complex catalytic systems. Lastly, first-principles MD methods are currently mostly used to investigate local reactive events. We hope that with the rise of data-driven methods and more efficient methods to describe the PES, first-principles MD methods will in the future also be able to describe longer length/time scale processes in catalysis. This might lead to a consistent dynamic description of all steps-diffusion, adsorption, and reaction-as they take place at the catalyst particle level.

Challenges in modelling dynamic processes in realistic nanostructured materials at operating conditions

The question is addressed in how far current modelling strategies are capable of modelling dynamic phenomena in realistic nanostructured materials at operating conditions. Nanostructured materials used in applications are far from perfect; they possess a broad range of heterogeneities in space and time extending over several orders of magnitude. Spatial heterogeneities from the subnanometre to the micrometre scale in crystal particles with a finite size and specific morphology, impact the material's dynamics. Furthermore, the material's functional behaviour is largely determined by the operating conditions. Currently, there exists a huge length-time scale gap between attainable theoretical length-time scales and experimentally relevant scales. Within this perspective, three key challenges are highlighted within the molecular modelling chain to bridge this length-time scale gap. Methods are needed that enable (i) building structural models for realistic crystal particles having mesoscale dimensions with isolated defects, correlated nanoregions, mesoporosity, internal and external surfaces; (ii) the evaluation of interatomic forces with quantum mechanical accuracy albeit at much lower computational cost than the currently used density functional theory methods and (iii) derivation of the kinetics of phenomena taking place in a multi-length-time scale window to obtain an overall view of the dynamics of the process. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

Molecular Views on Mechanisms of Brønsted Acid-Catalyzed Reactions in Zeolites

The Brønsted acidity of proton-exchanged zeolites has historically led to the most impactful applications of these materials in heterogeneous catalysis, mainly in the fields of transformations of hydrocarbons and oxygenates. Unravelling the mechanisms at the atomic scale of these transformations has been the object of tremendous efforts in the last decades. Such investigations have extended our fundamental knowledge about the respective roles of acidity and confinement in the catalytic properties of proton exchanged zeolites. The emerging concepts are of general relevance at the crossroad of heterogeneous catalysis and molecular chemistry. In the present review, emphasis is given to molecular views on the mechanism of generic transformations catalyzed by Brønsted acid sites of zeolites, combining the information gained from advanced kinetic analysis, in situ, and operando spectroscopies, and quantum chemistry calculations. After reviewing the current knowledge on the nature of the Brønsted acid sites themselves, and the key parameters in catalysis by zeolites, a focus is made on reactions undergone by alkenes, alkanes, aromatic molecules, alcohols, and polyhydroxy molecules. Elementary events of C-C, C-H, and C-O bond breaking and formation are at the core of these reactions. Outlooks are given to take up the future challenges in the field, aiming at getting ever more accurate views on these mechanisms, and as the ultimate goal, to provide rational tools for the design of improved zeolite-based Brønsted acid catalysts.

Simulations in the era of exascale computing

Exascale computers - supercomputers that can perform 1018 floating point operations per second - started coming online in 2022: in the United States, Frontier launched as the first public exascale supercomputer and Aurora is due to open soon; OceanLight and Tianhe-3 are operational in China; and JUPITER is due to launch in 2023 in Europe. Supercomputers offer unprecedented opportunities for modelling complex materials. In this Viewpoint, five researchers working on different types of materials discuss the most promising directions in computational materials science.

First-Generation Organic Reaction Intermediates in Zeolite Chemistry and Catalysis

Zeolite chemistry and catalysis are expected to play a decisive role in the next decade(s) to build a more decentralized renewable feedstock-dependent sustainable society owing to the increased scrutiny over carbon emissions. Therefore, the lack of fundamental and mechanistic understanding of these processes is a critical "technical bottleneck" that must be eliminated to maximize economic value and minimize waste. We have identified, considering this objective, that the chemistry related to the first-generation reaction intermediates (i.e., carbocations, radicals, carbenes, ketenes, and carbanions) in zeolite chemistry and catalysis is highly underdeveloped or undervalued compared to other catalysis streams (e.g., homogeneous catalysis). This limitation can often be attributed to the technological restrictions to detect such "short-lived and highly reactive" intermediates at the interface (gas-solid/solid-liquid); however, the recent rise of sophisticated spectroscopic/analytical techniques (including under in situ/operando conditions) and modern data analysis methods collectively compete to unravel the impact of these organic intermediates. This comprehensive review summarizes the state-of-the-art first-generation organic reaction intermediates in zeolite chemistry and catalysis and evaluates their existing challenges and future prospects, to contribute significantly to the "circular carbon economy" initiatives.

Carbocation chemistry confined in zeolites: spectroscopic and theoretical characterizations

Acid-catalyzed reactions inside zeolites are one type of broadly applied industrial reactions, where carbocations are the most common intermediates of these reaction processes, including methanol to olefins, alkene/aromatic alkylation, and hydrocarbon cracking/isomerization. The fundamental research on these acid-catalyzed reactions is focused on the stability, evolution, and lifetime of carbocations under the zeolite confinement effect, which greatly affects the efficiency, selectivity and deactivation of zeolite catalysts. Therefore, a profound understanding of the carbocations confined in zeolites is not only beneficial to explain the reaction mechanism but also drive the design of new zeolite catalysts with ideal acidity and cages/channels. In this review, we provide both an in-depth understanding of the stabilization of carbocations by the pore confinement effect and summary of the advanced characterization methods to capture carbocations in zeolites, including UV-vis spectroscopy, solid-state NMR, fluorescence microscopy, IR spectroscopy and Raman spectroscopy. Also, we clarify the relationship between the activity and stability of carbocations in zeolite-catalyzed reactions, and further highlight the role of carbocations in various hydrocarbon conversion reactions inside zeolites with diverse frameworks and varying acidic properties.