QM/MM study of the stability of dimethyl ether in zeolites H-ZSM-5 and H-Y

The methanol-to-hydrocarbons (MTH) process transforms C1 carbon sources to higher hydrocarbons, but details of the mechanism that leads to the formation of the first carbon–carbon bond remain unclear. Here, we present a computational investigation of how a crucial intermediate, dimethyl ether (DME), interacts with different zeolite catalysts (H-ZSM-5, H-Y) to gain insight into the initial stages in the MTH process. We use QM/MM computational simulations to model the conversion of methanol to DME in H-ZSM-5, which is a well characterised and important reaction intermediate. We analyse and compare the stability of DME on several acid sites in H-ZSM-5 and H-Y, and show that the more acidic and open “intersection sites” in the H-ZSM-5 framework are able to bond strongest with DME, with complete deprotonation of the acid site occurring. The conversion of methanol to DME in H-ZSM-5 is calculated as requiring a higher activation energy than framework methoxylation, which indicates that a stepwise (indirect) mechanism, through a methoxy intermediate, is the most likely route to DME formation during the initiation of the MTH process

QM/MM study of the reactivity of zeolite bound methoxy and carbene groups

The conversion of methanol-to-hydrocarbons (MTH) is known to occur via an autocatalytic process in zeolites, where framework-bound methoxy species play a pivotal role, especially during catalyst induction. Recent NMR and FT-IR experimental studies suggest that methoxylated zeolites are able to produce hydrocarbons by a mechanism involving carbene migration and association. In order to understand these observations, we have performed QM/MM computational investigations on a range of reaction mechanisms for the reaction of zeolite bound methoxy and carbene groups, which are proposed to initiate hydrocarbon formation in the MTH process. Our simulations demonstrate that it is kinetically unfavourable for methyl species to form on the framework away from the zeolite acid site, and both kinetically and thermodynamically unfavourable for methyl groups to migrate through the framework and aggregate around an acid site. Formation of carbene moieties was considered as an alternative pathway to the formation of C–C bonds; however, the reaction energy for conversion of a methyl to a carbene is unfavourable. Metadynamics simulations help confirm further that methyl species at the framework acid sites would be more reactive towards formed C_{2+} species, rather than inter-framework migration, and that the role of carbenes in the formation of the first C–C bond will be via a concerted type of mechanism rather than stepwise

Methanol loading dependent methoxylation in zeolite H-ZSM-5

We evaluate the effect of the number of methanol molecules per acidic site of H-ZSM-5 on the methoxylation reaction at room temperature by applying operando diffuse reflectance infrared Fourier transformed spectroscopy (DRIFTS) and mass spectrometry (MS), which capture the methoxylation reaction by simultaneously probing surface adsorbed species and reaction products, respectively. To this end, the methanol loading in H-ZSM-5 (Si/Al z 25) pores is systematically varied between 32, 16, 8 and 4 molecules per unit cell, which corresponds to 8, 4, 2 and 1 molecules per Brønsted acidic site, respectively. The operando DRIFTS/MS data show that the room temperature methoxylation depends on the methanol loading: the higher the methanol loading, the faster the methoxylation. Accordingly, the reaction is more than an order of magnitude faster with 8 methanol molecules per Brønsted acidic site than that with 2 molecules, as evident from the evolution of the methyl rock band of the methoxy species and of water as a function of time. Significantly, no methoxylation is observed with #1 molecule per Brønsted acidic site. However, hydrogen bonded methanol occurs across all loadings studied, but the structure of hydrogen bonded methanol also depends on the loading. Methanol loading of #1 molecule per acidic site leads to the formation of hydrogen bonded methanol with no proton transfer (i.e. neutral geometry), while loading $2 molecules per acidic site results in a hydrogen bonded methanol with a net positive charge on the adduct (protonated geometry). The infrared vibrational frequencies of methoxy and hydrogen bonded methanol are corroborated by Density Functional Theory (DFT) calculations. Both the experiments and calculations reflect the methoxy bands at around 940, 1180, 2868–2876 and 2980–2973 cm�1 which correspond to n(C–O), r(CH3), ns(C–H) and nas(C–H), respectively

Computational QM/MM investigation of the adsorption of MTH active species in H-Y and H-ZSM-5

The transformation of methanol-to-hydrocarbons (MTH) has significant potential as a route to synthesise low-cost fuels; however, the initial stages of the zeolite catalysed MTH process are not well understood. Here, we use hybrid quantum- and molecular-mechanical (QM/MM) embedded-cluster simulations to develop our understanding of the interaction between methanol and the zeolite catalysts H-ZSM-5, and for comparison, the larger pore H-Y. Energies and structures, calculated using hybrid-level density functional theory (hybrid-DFT) and higher-level correlated methods, are compared with previous experimental and computational results. We show that hydrogen-bonds between methanol adsorbates, formed through polarizable O-H bonds, substantially influence the adsorption energetics, structural parameters and vibrational frequencies. Our observations are extended by considering polar solvent molecules in the environment, with the presence of both water or methanol around the adsorption site leading to barrier-less transfer of the zeolite proton to an adsorbed methanol, which will significantly influence the reactivity of the adsorbed methanol