Karen Hemelsoet

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.

TAMkin: A Versatile Package for Vibrational Analysis and Chemical Kinetics

TAMkin is a program for the calculation and analysis of normal modes, thermochemical properties and chemical reaction rates. At present, the output from the frequently applied software programs ADF, CHARMM, CPMD, CP2K, Gaussian, Q-Chem, and VASP can be analyzed. The normal-mode analysis can be performed using a broad variety of advanced models, including the standard full Hessian, the Mobile Block Hessian, the Partial Hessian Vibrational approach, the Vibrational Subsystem Analysis with or without mass matrix correction, the Elastic Network Model, and other combinations. TAMkin is readily extensible because of its modular structure. Chemical kinetics of unimolecular and bimolecular reactions can be analyzed in a straightforward way using conventional transition state theory, including tunneling corrections and internal rotor refinements. A sensitivity analysis can also be performed, providing important insight into the theoretical error margins on the kinetic parameters. Two extensive examples demonstrate the capabilities of TAMkin: the conformational change of the biological system adenylate kinase is studied, as well as the reaction kinetics of the addition of ethene to the ethyl radical. The important feature of batch processing large amounts of data is highlighted by performing an extended level of theory study, which TAMkin can automate significantly.

Determining the storage, availability and reactivity of NH3 within Cu-Chabazite-based Ammonia Selective Catalytic Reduction systems

Three different types of NH3 species can be simultaneously present on Cu(2+)-exchanged CHA-type zeolites, commonly used in Ammonia Selective Catalytic Reduction (NH3-SCR) systems. These include ammonium ions (NH4(+)), formed on the Brønsted acid sites, [Cu(NH3)4](2+) complexes, resulting from NH3 coordination with the Cu(2+) Lewis sites, and NH3 adsorbed on extra-framework Al (EFAl) species, in contrast to the only two reacting NH3 species recently reported on Cu-SSZ-13 zeolite. The NH4(+) ions react very slowly in comparison to NH3 coordinated to Cu(2+) ions and are likely to contribute little to the standard NH3-SCR process, with the Brønsted groups acting primarily as NH3 storage sites. The availability/reactivity of NH4(+) ions can be however, notably improved by submitting the zeolite to repeated exchanges with Cu(2+), accompanied by a remarkable enhancement in the low temperature activity. Moreover, the presence of EFAl species could also have a positive influence on the reaction rate of the available NH4(+) ions. These results have important implications for NH3 storage and availability in Cu-Chabazite-based NH3-SCR systems.

Opposite regiospecific ring opening of 2-(cyanomethyl)aziridines by hydrogen bromide and benzyl bromide : experimental study and theoretical rationalization

Ring opening of 1-arylmethyl-2-(cyanomethyl)aziridines with HBr afforded 3-(arylmethyl)amino-4-bromobutyronitriles via regiospecific ring opening at the unsubstituted aziridine carbon. Previous experimental and theoretical reports show treatment of the same compounds with benzyl bromide to furnish 4-amino-3-bromobutanenitriles through ring opening at the substituted aziridine carbon. To gain insights into the regioselective preference with HBr, reaction paths have been analyzed with computational methods. The effect of solvation was taken into account by the use of explicit solvent molecules. Geometries were determined at the B3LYP/6-31++G(d,p) level of theory, and a Grimme-type correction term was included for long-range dispersion interactions; relative energies were refined with the meta-hybrid MPW1B95 functional. Activation barriers confirm preference for ring opening at the unsubstituted ring carbon for HBr. HBr versus benzyl bromide ring opening was analyzed through comparison of the electronic structure of corresponding aziridinium intermediates. Although the electrostatic picture fails to explain the opposite regiospecific nature of the reaction, frontier molecular orbital analysis of LUMOs and nucleophilic Fukui functions show a clear preference of attack for the substituted aziridine carbon in the benzyl bromide case and for the unsubstituted aziridine carbon in the HBr case, successfully rationalizing the experimentally observed regioselectivity.

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.

Mechanistic Studies on Chabazite‐Type Methanol‐to‐Olefin Catalysts: Insights from Time‐Resolved UV/Vis Microspectroscopy Combined with Theoretical Simulations

The formation and nature of active sites for methanol conversion over solid acid catalyst materials are studied by using a unique combined spectroscopic and theoretical approach. A working catalyst for the methanol‐to‐olefin conversion has a hybrid organic–inorganic nature in which a cocatalytic organic species is trapped in zeolite pores. As a case study, microporous materials with the chabazite topology, namely, H‐SAPO‐34 and H‐SSZ‐13, are considered with trapped (poly)aromatic species. First‐principle rate calculations on methylation reactions and in situ UV/Vis spectroscopy measurements are performed. The theoretical results show that the structure of the organic compound and zeolite composition determine the methylation rates: 1) the rate increases by 6 orders of magnitude if more methyl groups are added on benzenic species, 2) transition state selectivity occurs for organic species with more than one aromatic core and bearing more than three methyl groups, 3) methylation rates for H‐SSZ‐13 are approximately 3 orders of magnitude higher than on H‐SAPO‐34 owing to its higher acidity. The formation of (poly)aromatic cationic compounds can be followed by using in situ UV/Vis spectroscopy because these species yield characteristic absorption bands in the visible region of the spectrum. We have monitored the growth of characteristic peaks and derived activation energies of formation for various sets of (poly)aromatic compounds trapped in the zeolite host. The formation–activation barriers deduced by using UV/Vis microspectroscopy correlate well with the activation energies for the methylation of the benzenic species and the lower methylated naphthalenic species. This study shows that a fundamental insight at the molecular level can be obtained by using a combined in situ spectroscopic and theoretical approach for a complex catalyst of industrial relevance.

First-Principles Study of Antisite Defect Configurations in ZnGa2O4:Cr Persistent Phosphors.

Zinc gallate doped with chromium is a recently developed near-infrared emitting persistent phosphor, which is now extensively studied for in vivo bioimaging and security applications. The precise mechanism of this persistent luminescence relies on defects, in particular, on antisite defects and antisite pairs. A theoretical model combining the solid host, the dopant, and/or antisite defects is constructed to elucidate the mutual interactions in these complex materials. Energies of formation as well as dopant, and defect energies are calculated through density-functional theory simulations of large periodic supercells. The calculations support the chromium substitution on the slightly distorted octahedrally coordinated gallium site, and additional energy levels are introduced in the band gap of the host. Antisite pairs are found to be energetically favored over isolated antisites due to significant charge compensation as shown by calculated Hirshfeld-I charges. Significant structural distortions are found around all antisite defects. The local Cr surrounding is mainly distorted due to a ZnGa antisite. The stability analysis reveals that the distance between both antisites dominates the overall stability picture of the material containing the Cr dopant and an antisite pair. The findings are further rationalized using calculated densities of states and Hirshfeld-I charges.

Polycaprolactone and polycaprolactone/chitosan nanofibres functionalised with the pH-sensitive dye Nitrazine Yellow

Nanofibres functionalised with pH-sensitive dyes could greatly contribute to the development of stimuli-responsive materials. However, the application of biocompatible polymers is vital to allow for their use in (bio)medical applications. Therefore, this paper focuses on the development and characterisation of pH-sensitive polycaprolactone (PCL) nanofibrous structures and PCL/chitosan nanofibrous blends with 20% chitosan. Electrospinning with added Nitrazine Yellow molecules proved to be an excellent method resulting in pH-responsive non-wovens. Unlike the slow and broad response of PCL nanofibres (time lag of more than 3h), the use of blends with chitosan led to an increased sensitivity and significantly reduced response time (time lag of 5 min). These important effects are attributed to the increased hydrophilic nature of the nanofibres containing chitosan. Computational calculations indicated stronger interactions, mainly based on electrostatic interactions, of the dye with chitosan (ΔG of -132.3 kJ/mol) compared to the long-range interactions with PCL (ΔG of -35.6 kJ/mol), thus underpinning our experimental observations. In conclusion, because of the unique characteristics of chitosan, the use of PCL/chitosan blends in pH-sensitive biocompatible nanofibrous sensors is crucial.

Identification of Intermediates in Zeolite-Catalyzed Reactions by In Situ UV/Vis Microspectroscopy and a Complementary Set of Molecular Simulations

The optical absorption properties of (poly)aromatic hydrocarbons occluded in a nanoporous environment were investigated by theoretical and experimental methods. The carbonaceous species are an essential part of a working catalyst for the methanol-to-olefins (MTO) process. In situ UV/Vis microscopy measurements on methanol conversion over the acidic solid catalysts H-SAPO-34 and H-SSZ-13 revealed the growth of various broad absorption bands around 400, 480, and 580 nm. The cationic nature of the involved species was determined by interaction of ammonia with the methanol-treated samples. To determine which organic species contribute to the various bands, a systematic series of aromatics was analyzed by means of time-dependent density functional theory (TDDFT) calculations. Static gas-phase simulations revealed the influence of structurally different hydrocarbons on the absorption spectra, whereas the influence of the zeolitic framework was examined by using supramolecular models within a quantum mechanics/molecular mechanics framework. To fully understand the origin of the main absorption peaks, a molecular dynamics (MD) study on the organic species trapped in the inorganic host was essential. During such simulation the flexibility is fully taken into account and the effect on the UV/Vis spectra is determined by performing TDDFT calculations on various snapshots of the MD run. This procedure allows an energy absorption scale to be provided and the various absorption bands determined from in situ UV/Vis spectra to be assigned to structurally different species.