The conversion of methanol to olefins (MTO) over a heterogeneous nanoporous catalyst material is a highly complex process involving a cascade of elementary reactions. The elucidation of the reaction mechanisms leading to either the desired production of ethene and/or propene or undesired deactivation has challenged researchers for many decades. Clearly, catalyst choice, in particular topology and acidity, as well as the specific process conditions determine the overall MTO activity and selectivity; however, the subtle balances between these factors remain not fully understood. In this review, an overview of proposed reaction mechanisms for the MTO process is given, focusing on the archetypal MTO catalysts, H-ZSM-5 and H-SAPO-34. The presence of organic species, that is, the so-called hydrocarbon pool, in the inorganic framework forms the starting point for the majority of the mechanistic routes. The combination of theory and experiment enables a detailed description of reaction mechanisms and corresponding reaction intermediates. The identification of such intermediates occurs by different spectroscopic techniques, for which theory and experiment also complement each other. Depending on the catalyst topology, reaction mechanisms proposed thus far involve aromatic or aliphatic intermediates. Ab initio simulations taking into account the zeolitic environment can nowadays be used to obtain reliable reaction barriers and chemical kinetics of individual reactions. As a result, computational chemistry and by extension computational spectroscopy have matured to the level at which reliable theoretical data can be obtained, supplying information that is very hard to acquire experimentally. Special emphasis is given to theoretical developments that open new perspectives and possibilities that aid to unravel a process as complex as methanol conversion over an acidic porous material.

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.

Methylations of ethene, propene, and butene by methanol over the acidic microporous H-ZSM-5 catalyst are studied by means of state of the art computational techniques, to derive Arrhenius plots and rate constants from first principles that can directly be compared with the experimental data. For these key elementary reactions in the methanol to hydrocarbons (MTH) process, direct kinetic data became available only recently [J. Catal.2005, 224, 115-123; J. Catal.2005, 234, 385-400]. At 350 °C, apparent activation energies of 103, 69, and 45 kJ/mol and rate constants of 2.6 × 10(-4), 4.5 × 10(-3), and 1.3 × 10(-2) mol/(g h mbar) for ethene, propene, and butene were derived, giving following relative ratios for methylation k(ethene)/k(propene)/k(butene) = 1:17:50. In this work, rate constants including pre-exponential factors are calculated which give very good agreement with the experimental data: apparent activation energies of 94, 62, and 37 kJ/mol for ethene, propene, and butene are found, and relative ratios of methylation k(ethene)/k(propene)/k(butene) = 1:23:763. The entropies of gas phase alkenes are underestimated in the harmonic oscillator approximation due to the occurrence of internal rotations. These low vibrational modes were substituted by manually constructed partition functions. Overall, the absolute reaction rates can be calculated with near chemical accuracy, and qualitative trends are very well reproduced. In addition, the proposed scheme is computationally very efficient and constitutes significant progress in kinetic modeling of reactions in heterogeneous catalysis.

Ga/H-MFI was prepared by vapor-phase reaction of GaCl 3 with Brønsted acid O−H groups in dehydrated H-MFI zeolite. The resulting [GaCl 2 ] + cations in the as-exchanged zeolite are treated in H 2 at 823 K to stoichiometrically remove Cl ligands and form [GaH 2 ] + cations. Subsequent oxidation in O 2 and characterization by IR spectroscopy and NH 3 -temperature-programmed desorption (TPD) suggests that, for Ga/Al ratios ≤0.3, Ga 3+ exists predominantly as [Ga(OH) 2 ] + −H + cation pairs and to a lesser degree as [Ga(OH)] 2+ cations at low Ga/Al ratios (∼0.1); while both species are associated with proximate cationexchange sites, calculated free energies of formation suggest that [Ga(OH)] 2+ cations are more stable on cation-exchange sites associated with NNN (next-nearest neighbor) framework Al atoms than on those associated with NNNN (next-next-nearest neighbor) framework Al atoms. Ga K-edge X-ray Absorption Near Edge Spectroscopy (XANES) measurements indicate that, under oxidizing conditions and for all Ga/Al ratios, all Ga species are in the +3 oxidation state and are tetrahedrally coordinated to 4 O atoms. Fourier analysis of Ga K-edge Extended Xray Absorption Fine Structure (EXAFS) data supports the conclusion that Ga 3+ is present predominantly as [Ga(OH) 2 ] + cations (or [Ga(OH) 2 ] + −H + cation pairs). For Ga/Al ratios ≤0.3, wavelet analysis of EXAFS data provide evidence for backscattering from nearest neighboring O atoms and from next-nearest neighboring framework Al atoms. For Ga/Al > 0.3, backscattering from next-nearest neighboring Ga atoms is also evident, characteristic of GaO x species. Upon reduction in H 2 , the oxidized Ga 3+ species produce [Ga(OH)H] + −H + cation pairs, [GaH 2 ] + −H + cation pairs, and [GaH] 2+ cations. Computed phase diagrams indicate that the thermodynamic stability of the reduced Ga 3+ species depends sensitively on temperature, Al−Al interatomic distance, and H 2 and H 2 O partial pressures. For Ga/Al ratios ≤0.2, it is concluded that [GaH 2 ] + −H + cation pairs and [GaH] 2+cations are the predominant species present in Ga/H-MFI reduced above 673 K in 10 5 Pa H 2 and in the absence of water vapor.

Catalytic cracking of alkenes takes place at elevated temperatures in the order of 773–833 K. In this work, the nature of the reactive intermediates at typical reaction conditions is studied in H-ZSM-5 using a complementary set of modeling tools. Ab initio static and molecular dynamics simulations are performed on different C4single bond C5 alkene cracking intermediates to identify the reactive species in terms of temperature. At 323 K, the prevalent intermediates are linear alkoxides, alkene π-complexes and tertiary carbenium ions. At a typical cracking temperature of 773 K, however, both secondary and tertiary alkoxides are unlikely to exist in the zeolite channels. Instead, more stable carbenium ion intermediates are found. Branched tertiary carbenium ions are very stable, while linear carbenium ions are predicted to be metastable at high temperature. Our findings confirm that carbenium ions, rather than alkoxides, are reactive intermediates in catalytic alkene cracking at 773 K