Cracking of Hydrocarbons on Zeolite Catalysts: Density Functional and Hartree−Fock Calculations on the Mechanism of the β-Scission Reaction

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.

Trimethyloxonium ion – a zeolite confined mobile and efficient methyl carrier at low temperatures: a DFT study coupled with microkinetic analysis

The reaction network of ethene methylation over H-ZSM-5, including methanol dehydration, ethene methylation, and C3H7+ conversion, is investigated by employing a multiscale approach combining DFT calculations and microkinetic modeling.

Dynamic co-catalysis of Au single atoms and nanoporous Au for methane pyrolysis

Nanocatalysts and single-atom catalysts are both vital for heterogeneous catalysis. They are recognized as two different categories of catalysts. Nevertheless, recent theoretical works have indicated that Au nanoparticles/clusters release Au single atoms in CO oxidation, and they co-catalyze the oxidation. However, to date, neither experimental evidence for the co-catalysis nor direct observations on any heterogeneous catalysis process of single-atom catalysts are reported. Here, the dynamic process of nanoporous Au to catalyze methane pyrolysis is monitored by in situ transmission electron microscopy with high spatial–temporal resolutions. It demonstrates that nanoporous Au surfaces partially disintegrate, releasing Au single atoms. As demonstrated by DFT calculation, the single atoms could co-catalyze the reaction with nanoporous Au. Moreover, the single atoms dynamically aggregate into nanoparticles, which re-disintegrate back to single atoms. This work manifests that under certain conditions, the heterogeneous catalysis processes of nanocatalysts and single-atom catalysts are not independent, where their dynamic co-catalysis exists.

Nanocatalysts and single‐atom catalysts are generally considered as two categories with distinct performances. Here, in situ TEM study of catalytic methane pyrolysis over nanoporous Au reveals a highly dynamic process where co‐catalysis exists among various catalyst forms.

Theoretical insight into the single-atom catalytic mechanism of CeO2-supported Ag catalysts in CO oxidation

Revealing the accurate active center structure and the functional mechanism of CeO2-supported Ag catalysts during catalysis is extremely important for their accurate synthesis. In this work, a series of AgnCeO2 (n = 1, 2, 3, 4 and 10) model catalysts was constructed, and a DFT investigation of the reaction mechanism of CO oxidation, as a probe reaction on those catalysts, was carried out. It was found that the entire catalytic reaction was completed coordinately by Ag, lattice O and O vacancies, which could be considered as the active centers. Noticeably, the mobility of Ag atoms played an important role in the reaction process, leading to the observation of a single-atom catalytic mechanism, wherein a series of single Ag atomic species was formed during the reaction, which was beneficial to CO oxidation. With the completion of some elementary reactions, the single Ag formed during the migration of CO-Ag could return to the Ag cluster again. As expected, the single-AgCeO2 catalyst exhibited extremely high activity due to the absence of the binding effect of Ag-Ag. Nevertheless, the AgnCeO2 (n > 1) catalysts showed similar catalytic activity, which was slightly worse than that of single AgCeO2, indicating that the size effect of the Ag cluster was not obvious. These results provide the theoretical basis for further understanding the functional mechanism of the AgnCeO2 catalyst and are helpful for designing various catalysts with tailored functionalities.

Identification of the Reaction Sequence of the MTO Initiation Mechanism Using Ab Initio-Based Kinetics

The initiation of the methanol-to-olefins (MTO) process is investigated using a multiscale modeling approach where more than 100 ab initio computed (MP2:DFT) rate constants for H-SSZ-13 are used in a batch reactor model. The investigated reaction network includes the mechanism for initiation (42 steps) and a representative part of the autocatalytic olefin cycle (63 steps). The simulations unravel the dominant initiation pathway for H-SSZ-13: dehydrogenation of methanol to CO is followed by CO-methylation leading to the formation of the first C-C bond in methyl acetate despite high barriers of >200 kJ/mol. Our multiscale approach is able to shed light on the reaction sequence that ultimately leads to olefin formation and strikingly demonstrates that only with a full reactor model that includes autocatalysis with olefins as cocatalysts is one able to understand the initiation mechanism on the atomic scale. Importantly, the model also shows that autocatalysis takes over long before significant amounts of olefins are formed, thus guiding the interpretation of experimental results.

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.

Olefin methylation and cracking reactions in H-SSZ-13 investigated with ab initio and DFT calculations

The olefin cycle of the methanol-to-olefins process is investigated for the zeolite H-SSZ-13 using periodic, van-der-Waals corrected DFT calculations, together with MP2 corrections derived from cluster models, which are essential for accurate barriers.

Catalyzed β scission of a carbenium ion III — Scission observed in ab initio molecular dynamics simulations

The β scission (cracking) of branched carbenium ions have been observed in molecular dynamics simulations, possibly for the first time. Simulations were performed with molecular dynamics based on PW91 density functional theory, and which included three-dimensional periodic boundary replication of the unit cell to mimic long-range bulk effects. A rising-temperature algorithm was used to encourage reaction within the narrow time windows (∼10 ps) of the simulations. Twenty-eight simulations were performed, featuring alkyl ions in three different catalytic systems: the ionic liquid, [(C5H5NH+)5(Al2Cl7−)6]−, the chabazite zeolite, [AlSi23O48]−, and the chabazite zeolite, [Al4Si20O45(OH)3]−. Twenty-four runs began with unbranched sec-n-alkyl ions, but only one exhibited β scission, and only after branching to a tertiary ion and under extreme heating. In contrast, the four simulations that began with branched alkyl ions were all successful in demonstrating β scission at lower temperatures: 2,4,4-trimethyl-2-pentyl ion and 2,4-dimethyl-2-hexyl ion in each of the first two catalysts. The lifetimes of desorbed alkyl ions in the chabazite models were < 5 ps at 1000–1500 K. The β scission results support the classical Weitkamp et al. ( Appl. Catal. 1983, 8, 123 ) mechanism over the nonclassical Sie ( Ind. Eng. Chem. Res. 1992, 31, 1881 ) and the chemisorping Kazansky et al. ( J. Catal. 1989, 119, 108 ) mechanisms.

Palladium-Catalyzed Oxidative Activation of Arylcyclopropanes

Palladium chloride-catalyzed intramolecular activation of electroneutral cyclopropane derivatives results in cleavage of the cyclopropane ring followed by formation of heterocyclic derivatives. Phenols, carboxylic acids, and amide groups were considered as substituents ortho to the cyclopropane ring in this catalytic activation chemistry. The regioselectivity observed in the case of amide-containing substrates was different from that of carboxylic acid-containing substrates, ruling out simple cyclopropane isomerization followed by a Wacker oxidation as the mechanistic pathway. [reaction: see text]

Catalyzed β scission of a carbenium ion II — Variations leading to a general mechanism

The β-scission mechanisms of catalytically chemisorbed carbenium ions are further investigated using density functional theory computations and explicit-contact modelling, but with slightly larger catalyst fragment models than in our previous work. Some variations are seen, including the existence of formal one-step and three-step (rather than two-step) mechanisms. The activation barriers are most affected by the basicity of the catalyst model than by any other characteristics: the stronger the base, the greater the barrier. A general mechanism for β scission is presented, as are the specific mechanisms for all the step variations observed from computations to date.Key words: C–C bond fission, β scission, carbenium ion, catalysis, chloroaluminate, mechanism.

Reactivity of Alkanes on Zeolites: A Computational Study of Propane Conversion Reactions

In this work, quantum chemical methods were used to study propane conversion reactions on zeolites; these reactions included protolytic cracking, primary hydrogen exchange, secondary hydrogen exchange, and dehydrogenation reactions. The reactants, products, and transition-state structures were optimized at the B3LYP/6-31G level and the energies were calculated with CBS-QB3, a complete basis set composite energy method. The computed activation barriers were 62.1 and 62.6 kcal/mol for protolytic cracking through two different transition states, 30.4 kcal/mol for primary hydrogen exchange, 29.8 kcal/mol for secondary hydrogen exchange, and 76.7 kcal/mol for dehydrogenation reactions. The effects of basis set for the geometry optimization and zeolite acidity on the reaction barriers were also investigated. Adding extra polarization and diffuse functions for the geometry optimization did not affect the activation barriers obtained with the composite energy method. The largest difference in calculated activation barriers is within 1 kcal/mol. Reaction activation barriers do change as zeolite acidity changes, however. Linear relationships were found between activation barriers and zeolite deprotonation energies. Analytical expressions for each reaction were proposed so that accurate activation barriers can be obtained when using different zeolites as catalysts, as long as the deprotonation energies are first acquired.

Catalyzed β scission of a carbenium ion — Mechanistic differences from varying catalyst basicity

The β-scission mechanism of physisorbed and chemisorbed pentenium ions, as catalyzed by AlH2(OH)2– and by AlHCl3– anions, was investigated using density functional theory computations and explicit-contact modelling. A thorough search of intermediates was performed for each catalyst. On the aluminum chloride, β scission of an aliphatic, secondary carbenium ion featured chemisorbed and physisorbed ion intermediates, while on the aluminum hydroxide, β scission featured chemisorbed ions but physisorbed neutral species. The importance of this work is its demonstration of a qualitatively different mechanism, with qualitatively different intermediates, due only to the different basicity of the two catalysts.Key words: C—C bond fission, β scission, carbenium ion, catalysis, mechanism.

Cluster and periodic calculations of the ethene protonation reaction catalyzed by theta-1 zeolite: influence of method, model size, and structural constraints

The protonation of ethene by three different acid sites of theta-1 zeolite was theoretically studied to analyze the extent and relevance of the following aspects of heterogeneous catalysis: the local geometry of the Brønsted acid site in a particular zeolite, the size of the cluster used to model the catalyst, the degree of geometry relaxation around the active site, and the effects related to medium- and long-range interactions between the reaction site and its environment. It has been found that while the reaction energy is very sensitive to the local geometry of the site, the activation energy is mainly affected by the methodology used and by electrostatic effects on account of the carbocationic nature of the transition state.