Facing the Challenges of Borderline Oxidation State Assignments Using State-of-the-Art Computational Methods.

The oxidation state (OS) of metals and ligands in inorganic complexes may be defined by carefully curated rules, such as from IUPAC, or by computational procedures such as the effective oxidation state (EOS) or localized orbital bonding analysis (LOBA). Such definitions typically agree for systems with simple ionic bonding and innocent ligands but may disagree as the boundary between ionic and covalent bonds is approached, or as the role of ligand noninnocence becomes nontrivial, or high oxidation states of metals are supported by heavy dative bonding, and so on. This work systematically compares IUPAC, EOS, and LOBA across a series of complexes where OS assignment is challenging. These systems include high-valent transition metal oxides, transition metal complexes with noninnocent ligands such as dithiolate and nitrosyl, metal sulfur dioxide adducts, and two transition metal carbene complexes. The differences in OS assignment by the three methods are carefully discussed, demonstrating the synergy between EOS and LOBA. In addition, a clarity index for LOBA OS assignments is introduced that provides an indication of whether or not its predictions are close to the ionic-covalent boundary.

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

This is a response to the paper published by S. A. Kadam, H. Li, R. F. Wormsbacher, A. Travert, Chem. Eur. J. 2018, 24, 5489. Key consistencies between our reported results and those reported in this work are also highlighted.

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 (kapp ), it is necessary to identify the equilibrium constant for adsorption into the reactant state (Kads-H+ ) and the intrinsic rate coefficient of the reaction (kint ) at reaction temperatures, since kapp is proportional to the product of Kads-H+ and kint . We show that Kads-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 C3 -C6 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 kapp and its components, Kads-H+ and kint .

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.

Enthalpy and entropy barriers explain the effects of topology on the kinetics of zeolite-catalyzed reactions.

The methylation of ethene, propene, and trans-2-butene on zeolites H-ZSM-58 (DDR), H-ZSM-22 (TON), and H-ZSM-5 (MFI) is studied to elucidate the particular influence of topology on the kinetics of zeolite-catalyzed reactions. H-ZSM-58 and H-ZSM-22 are found to display overall lower methylation rates compared to H-ZSM-5 and also different trends in methylation rates with increasing alkene size. These variations may be rationalized based on a decomposition of the free-energy barriers into enthalpic and entropic contributions, which reveals that the lower methylation rates on H-ZSM-58 and H-ZSM-22 have virtually opposite reasons. On H-ZSM-58, the lower methylation rates are caused by higher enthalpy barriers, owing to inefficient stabilization of the reaction intermediates in the large cage-like pores. On the other hand, on H-ZSM-22, the methylation rates mostly suffer from higher entropy barriers, because excessive entropy losses are incurred inside the narrow-channel structure. These results show that the kinetics of crucial elementary steps hinge on the balance between proper stabilization of the reaction intermediates inside the zeolite pores and the resulting entropy losses. These fundamental insights into their inner workings are indispensable for ultimately selecting or designing better zeolite catalysts.

Unraveling the reaction mechanisms governing methanol-to-olefins catalysis by theory and experiment.

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.

First principle kinetic studies of zeolite-catalyzed methylation reactions.

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.

Computational Modeling Predicts the Stability of Both Pd+ and Pd2+ Ion-Exchanged into H-CHA.

Passive NOx adsorbers (PNA) using Pd/zeolites have emerged as a promising solution for the reduction of cold-start emissions from vehicle exhaust. However, the nature of the active sites and the mechanisms underlying NOx adsorption in Pd/zeolites remain a subject of ongoing investigation. In this study, we employ quantum chemical simulations to investigate the structure of Pd species in cation-exchange sites at isolated Al and Al pairs in the 6-ring and 8-ring of the CHA framework, before the introduction of NOx. Our calculations show that the speciation of Pd in these exchange sites strongly depends on the precise Al arrangement within the framework, as well as the operating conditions. Ionically dispersed Pd is found to be the most favorable species over a wide range of oxidizing and reducing conditions. Small oligomers of PdO and metallic Pd do not appear to be competitive at either isolated Al or Al pairs. Notably, our calculations show that ion exchange sites other than next-next-nearest neighbor Al pairs in the 6-ring will be preferentially occupied by Pd+ instead of Pd2+. The stability of Pd+ in the zeolite environment is an interesting contrast with its rareness in molecular Pd compounds. Nonetheless, a detailed analysis of the electronic structure shows that predicted Pd oxidation states are consistent with chemical intuition for all complexes investigated in this study. We also discuss the potential ambiguity in Pd characterization provided by typical experimental techniques such as XANES, EXAFS and UV-VIS, and highlight the need for additional EPR spectroscopy studies to further elucidate the initial Pd speciation in zeolites for PNA applications.

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 C4-C8 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 B2 (3° → 2°) are expected to be dominant at operating conditions, whereas modes B1 (2° → 3°), C (2° → 2°), D2 (2° → 1°), and E2 (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.