Ab Initio Simulations Reveal that Reaction Dynamics Strongly Affect Product Selectivity for the Cracking of Alkanes over H-MFI

Selectivity descriptors of the catalytic n-hexane cracking process over 10-membered ring zeolites

Zeolite-mediated catalytic cracking of alkanes is pivotal in the petrochemical and refining industry, breaking down heavier hydrocarbon feedstocks into fuels and chemicals. Its relevance also extends to emerging technologies such as biomass and plastic valorization. Zeolite catalysts, with shape selectivity and selective adsorption capabilities, enhance efficiency and sustainability due to their well-defined network of pores, dimensionality, cages/cavities, and channels. This study focuses on the alkane cracking over 10-membered ring (10-MR) zeolites under industrially relevant conditions. Through a series of characterizations, including operando UV-vis spectroscopy and solid-state NMR spectroscopy, we intend to address mechanistic debates about the alkane cracking mechanism, aiming to understand the dependence of product selectivity on zeolite topologies. The findings highlight topology-dependent mechanisms, particularly the role of intersectional void spaces in zeolite ZSM-5, influencing aromatic-based product selectivity. This work provides a unique understanding of zeolite-catalyzed hydrocarbon conversion, linking alkane activation steps to the traditional hydrocarbon pool mechanism, contributing to the fundamental knowledge of this crucial industrial process.

Spatiotemporal Heterogeneity of Temperature and Catalytic Activation within Individual Catalyst Particles

Temperature is a critical parameter in chemical conversion, significantly affecting the reaction kinetics and thermodynamics. Measuring temperature inside catalyst particles of industrial interest (∼micrometers to millimeters), which is crucial for understanding the evolution of chemical dynamics at catalytic active sites during reaction and advancing catalyst designs, however, remains a big challenge. Here, we propose an approach combining two-photon confocal microscopy and state-of-the-art upconversion luminescence (UL) imaging to measure the spatiotemporal-resolved temperature within individual catalyst particles in the industrially significant methanol-to-hydrocarbons reaction. Specifically, catalyst particles containing zeolites and functional nanothermometers were fabricated using microfluidic chips. Our experimental results directly demonstrate that the zeolite density and particle size can alter the temperature distribution within a single catalyst particle. Importantly, the observed temperature heterogeneity plays a decisive role in the activation of the reaction intermediate and the utilization of active sites. We expect that this work opens a venue for unveiling the reaction mechanism and kinetics within industrial catalyst particles by considering temperature heterogeneity.

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.

A Review of Recent Developments in Molecular Dynamics Simulations of the Photoelectrochemical Water Splitting Process

In this review, we provide a short overview of the Molecular Dynamics (MD) method and how it can be used to model the water splitting process in photoelectrochemical hydrogen production. We cover classical non-reactive and reactive MD techniques as well as multiscale extensions combining classical MD with quantum chemical and continuum methods. Selected examples of MD investigations of various aqueous semiconductor interfaces with a special focus on TiO2 are discussed. Finally, we identify gaps in the current state-of-the-art where further developments will be needed for better utilization of MD techniques in the field of water splitting.

Towards operando computational modeling in heterogeneous catalysis

An increased synergy between experimental and theoretical investigations in heterogeneous catalysis has become apparent during the last decade. Experimental work has extended from ultra-high vacuum and low temperature towards operando conditions. These developments have motivated the computational community to move from standard descriptive computational models, based on inspection of the potential energy surface at 0 K and low reactant concentrations (0 K/UHV model), to more realistic conditions. The transition from 0 K/UHV to operando models has been backed by significant developments in computer hardware and software over the past few decades. New methodological developments, designed to overcome part of the gap between 0 K/UHV and operando conditions, include (i) global optimization techniques, (ii) ab initio constrained thermodynamics, (iii) biased molecular dynamics, (iv) microkinetic models of reaction networks and (v) machine learning approaches. The importance of the transition is highlighted by discussing how the molecular level picture of catalytic sites and the associated reaction mechanisms changes when the chemical environment, pressure and temperature effects are correctly accounted for in molecular simulations. It is the purpose of this review to discuss each method on an equal footing, and to draw connections between methods, particularly where they may be applied in combination.

How Chain Length and Branching Influence the Alkene Cracking Reactivity on H-ZSM-5

Catalytic alkene cracking on H-ZSM-5 involves a complex reaction network with many possible reaction routes and often elusive intermediates. Herein, advanced molecular dynamics simulations at 773 K, a typical cracking temperature, are performed to clarify the nature of the intermediates and to elucidate dominant cracking pathways at operating conditions. A series of C₄–C₈ alkene intermediates are investigated to evaluate the influence of chain length and degree of branching on their stability. Our simulations reveal that linear, secondary carbenium ions are relatively unstable, although their lifetime increases with carbon number. Tertiary carbenium ions, on the other hand, are shown to be very stable, irrespective of the chain length. Highly branched carbenium ions, though, tend to rapidly rearrange into more stable cationic species, either via cracking or isomerization reactions. Dominant cracking pathways were determined by combining these insights on carbenium ion stability with intrinsic free energy barriers for various octene β-scission reactions, determined via umbrella sampling simulations at operating temperature (773 K). Cracking modes A (3° → 3°) and B₂ (3° → 2°) are expected to be dominant at operating conditions, whereas modes B₁ (2° → 3°), C (2° → 2°), D₂ (2° → 1°), and E₂ (3° → 1°) are expected to be less important. All β-scission modes in which a transition state with primary carbocation character is involved have high intrinsic free energy barriers. Reactions starting from secondary carbenium ions will contribute less as these intermediates are short living at the high cracking temperature. Our results show the importance of simulations at operating conditions to properly evaluate the carbenium ion stability for β-scission reactions and to assess the mobility of all species in the pores of the zeolite.

Impact of Zeolite Structure on Entropic-Enthalpic Contributions to Alkane Monomolecular Cracking: An IR Operando Study

The monomolecular cracking rates of propane and n-butane over MFI, CHA, FER and TON zeolites were determined simultaneously with the coverage of active sites at reaction condition using IR operando spectroscopy. This allowed direct determination of adsorption thermodynamics and intrinsic rate parameters. The results show that the zeolite confinement mediates enthalpy-entropy trade-offs only at the adsorbed state, leaving the true activation energy insensitive to the zeolite or alkane structure while the activation entropy was found to increase with the confinement. Hence, relative cracking rates of alkanes within zeolite pores are mostly governed by activation entropy.

Understanding Brønsted-Acid Catalyzed Monomolecular Reactions of Alkanes in Zeolite Pores by Combining Insights from Experiment and Theory

Acidic zeolites are effective catalysts for the cracking of large hydrocarbon molecules into lower molecular weight products required for transportation fuels. However, the ways in which the zeolite structure affects the catalytic activity at Brønsted protons are not fully understood. One way to characterize the influence of the zeolite structure on the catalysis is to study alkane cracking and dehydrogenation at very low conversion, conditions for which the kinetics are well defined. To understand the effects of zeolite structure on the measured rate coefficient (k_app ), it is necessary to identify the equilibrium constant for adsorption into the reactant state (K_ads-H+ ) and the intrinsic rate coefficient of the reaction (k_int ) at reaction temperatures, since k_app is proportional to the product of K_ads-H+ and k_int . We show that K_ads-H+ cannot be calculated from experimental adsorption data collected near ambient temperature, but can, however, be estimated accurately from configurational-bias Monte Carlo (CBMC) simulations. Using monomolecular cracking and dehydrogenation of C₃ -C₆ alkanes as an example, we review recent efforts aimed at elucidating the influence of the acid site location and the zeolite framework structure on the observed values of k_app and its components, K_ads-H+ and k_int .

The challenge of catalyst prediction

New insights and successful use of computational catalysis is highlighted.

Towards molecular control of elementary reactions in zeolite catalysis by advanced molecular simulations mimicking operating conditions

Ingeniously selecting zeolite topology and acidity, reaction temperature and guest molecule loading enables tuning the reaction mechanism of zeolite-catalyzed methylation reactions.

Simultaneous realization of high catalytic activity and stability for catalytic cracking of n-heptane on highly exposed (010) crystal planes of nanosheet ZSM-5 zeolite

Nanosheet ZSM-5 zeolite exhibits outstanding reactivity and anti-coking stability for catalytic cracking of n-heptane.

A potential with low point charges for pure siliceous zeolites

A modified CHARMM force-field (ZHB potential) with low point charges for silica was previously proposed by Zimmerman et al. (J. Chem. Theory Comput. 2011, 7, 1695). The ZHB potential is advantageous for quantum mechanics/molecular mechanics simulations as it minimizes the electron spill-out problems. In the same spirit, here we propose a modified ZHB potential (MZHB) by reformulating its bonding potential, while retaining the nonbonding potential as in the ZHB force-field. We show that several structural and dynamic properties of silica, like the IR spectrum, distribution functions, mechanical properties, and negative thermal expansion computed using the MZHB potential agree well with experimental data. Further, transferability of MZHB is also tested for reproducing the crystallographic structures of several polymorphs of silica. © 2015 Wiley Periodicals, Inc.

Advances in theory and their application within the field of zeolite chemistry

Zeolites are versatile and fascinating materials which are vital for a wide range of industries, due to their unique structural and chemical properties, which are the basis of applications in gas separation, ion exchange and catalysis. Given their economic impact, there is a powerful incentive for smart design of new materials with enhanced functionalities to obtain the best material for a given application. Over the last decades, theoretical modeling has matured to a level that model guided design has become within reach. Major hurdles have been overcome to reach this point and almost all contemporary methods in computational materials chemistry are actively used in the field of modeling zeolite chemistry and applications. Integration of complementary modeling approaches is necessary to obtain reliable predictions and rationalizations from theory. A close synergy between experimentalists and theoreticians has led to a deep understanding of the complexity of the system at hand, but also allowed the identification of shortcomings in current theoretical approaches. Inspired by the importance of zeolite characterization which can now be performed at the single atom and single molecule level from experiment, computational spectroscopy has grown in importance in the last decade. In this review most of the currently available modeling tools are introduced and illustrated on the most challenging problems in zeolite science. Directions for future model developments will be given.

First principle chemical kinetics in zeolites: the methanol-to-olefin process as a case study

To optimally design next generation catalysts a thorough understanding of the chemical phenomena at the molecular scale is a prerequisite. Apart from qualitative knowledge on the reaction mechanism, it is also essential to be able to predict accurate rate constants. Molecular modeling has become a ubiquitous tool within the field of heterogeneous catalysis. Herein, we review current computational procedures to determine chemical kinetics from first principles, thus by using no experimental input and by modeling the catalyst and reacting species at the molecular level. Therefore, we use the methanol-to-olefin (MTO) process as a case study to illustrate the various theoretical concepts. This process is a showcase example where rational design of the catalyst was for a long time performed on the basis of trial and error, due to insufficient knowledge of the mechanism. For theoreticians the MTO process is particularly challenging as the catalyst has an inherent supramolecular nature, for which not only the Brønsted acidic site is important but also organic species, trapped in the zeolite pores, must be essentially present during active catalyst operation. All these aspects give rise to specific challenges for theoretical modeling. It is shown that present computational techniques have matured to a level where accurate enthalpy barriers and rate constants can be predicted for reactions occurring at a single active site. The comparison with experimental data such as apparent kinetic data for well-defined elementary reactions has become feasible as current computational techniques also allow predicting adsorption enthalpies with reasonable accuracy. Real catalysts are truly heterogeneous in a space- and time-like manner. Future theory developments should focus on extending our view towards phenomena occurring at longer length and time scales and integrating information from various scales towards a unified understanding of the catalyst. Within this respect molecular dynamics methods complemented with additional techniques to simulate rare events are now gradually making their entrance within zeolite catalysis. Recent applications have already given a flavor of the benefit of such techniques to simulate chemical reactions in complex molecular environments.

Experimental and Theoretical Methods in Kinetic Studies of Heterogeneously Catalyzed Reactions

This review aims to illustrate the potential of kinetic analysis in general and microkinetic modeling in particular for rational catalyst design. Both ab initio calculations and experiments providing intrinsic kinetic data allow us to assess the effects of catalytic properties and reaction conditions on the activity and selectivity of the targeted reactions. Three complementary approaches for kinetic analysis of complex reaction networks are illustrated, using select examples of acid zeolite–catalyzed reactions from the authors' recent work. Challenges for future research aimed at defining targets for synthesis strategies that enable us to tune zeolite properties are identified.

Conversion of CO2 and C2H6 to propanoic acid over a Au-exchanged MCM-22 zeolite

Finding novel catalysts for the direct conversion of CO2 to fuels and chemicals is a primary goal in energy and environmental research. In this work, density functional theory (DFT) is used to study possible reaction mechanisms for the conversion of CO2 and C2H6 to propanoic acid over a gold-exchanged MCM-22 zeolite catalyst. The reaction begins with the activation of ethane to produce a gold ethyl hydride intermediate. Hydrogen transfers to the framework oxygen leads then to gold ethyl adsorbed on the Brønsted-acid site. The energy barriers for these steps of ethane activation are 9.3 and 16.3 kcal mol(-1), respectively. Two mechanisms of propanoic acid formation are investigated. In the first one, the insertion of CO2 into the Au-H bond of the first intermediate yields gold carboxyl ethyl as subsequent intermediate. This is then converted to propanoic acid by forming the relevant C-C bond. The activation energy of the rate-determining step of this pathway is 48.2 kcal mol(-1). In the second mechanism, CO2 interacts with gold ethyl adsorbed on the Brønsted-acid site. Propanoic acid is formed via protonation of CO2 by the Brønsted acid and the simultaneous formation of a bond between CO2 and the ethyl group. The activation energy there is 44.2 kcal mol(-1), favoring this second pathway at least at low temperatures. Gold-exchanged MCM-22 zeolite can therefore, at least in principle, be used as the catalyst for producing propanoic acid from CO2 and ethane.