Kristof De Wispelaere

Herein, we propose a novel computational protocol that enables calculating free energies with improved accuracy by combining the best available techniques for enthalpy and entropy calculation. While the entropy is described by enhanced sampling molecular dynamics techniques, the energy is calculated using ab initio methods. We apply the method to assess the stability of isobutene adsorption intermediates in the zeolite H-SSZ-13, a prototypical problem that is computationally extremely challenging in terms of calculating enthalpy and entropy. We find that at typical operating conditions for zeolite catalysis (400 °C), the physisorbed π-complex, and not the tertiary carbenium ion as often reported, is the most stable intermediate. This method paves the way for sampling-based techniques to calculate the accurate free energies in a broad range of chemistry-related disciplines, thus presenting a big step forward toward predictive modeling.

The methanol-to-hydrocarbons reaction refers collectively to a series of important industrial catalytic processes to produce either olefins or gasoline. Mechanistically, methanol conversion proceeds through a ‘pool’ of hydrocarbon species. For the methanol-to-olefins process, these species can be delineated broadly into ‘desired’ lighter olefins and ‘undesired’ heavier fractions that cause deactivation in a matter of hours. The crux in further catalyst optimization is the ability to follow the formation of carbonaceous species during operation. Here, we report the combined results of an operando Kerr-gated Raman spectroscopic study with state-of-the-art operando molecular simulations, which allowed us to follow the formation of hydrocarbon species at various stages of methanol conversion. Polyenes are identified as crucial intermediates towards formation of polycyclic aromatic hydrocarbons, with their fate determined largely by the zeolite topology. Notably, we provide the missing link between active and deactivating species, which allows us to propose potential design rules for future-generation catalysts.

The methanol-to-hydrocarbons process is known to proceed autocatalytically in H-ZSM-5 after an induction period where framework methoxy species are formed. In this work, we provide mechanistic insight into the framework methylation within H-ZSM-5, at high methanol loadings and varying acid site density, by means of first principles molecular dynamics simulations. The molecular dynamics simulations show that stable methanol clusters form in the zeolite pores, and these clusters commonly deprotonate the active site, however the cluster size is dependent on temperature and acid site density. Enhanced sampling molecular dynamics simulations give evidence that the barrier for methanol conversion is significantly affected by the neighbourhood of an additional acid site, suggesting that cooperative effects influence methanol clustering and reactivity. The insights obtained are important steps in optimising the catalyst and engineering the induction period of the methanol-to-hydrocarbons process.

When exposed to steam, zeolite catalysts are irreversibly deactivated by loss of acidity and framework degradation caused by dealumination. Steaming typically occurs at elevated temperatures, making it challenging to investigate the mechanism with most approaches. Herein, we follow the dynamics of zeolite dealumination in situ, in the presence of a realistic loading of water molecules by means of enhanced sampling molecular dynamics simulations. H-SSZ-13 zeolite is chosen as a target system. Monte Carlo simulations predict a loading of more than 3 water molecules per unit cell at representative steaming conditions (450 °C, 1 bar steam). Our results show that a higher water loading lowers the free energy barrier of dealumination, as water molecules cooperate to facilitate hydrolysis of Al–O bonds. We find free energies of activation for dealumination that agree well with the available experimental measurements. Clearly, the use of enhanced sampling molecular dynamics yields a major step forward in the molecular level understanding of the dealumination; insight which is very hard to derive experimentally.

The active site in ethene oligomerization catalyzed by Ni-zeolites is proposed to be a mobile Ni(II) complex, based on density functional theory-based molecular dynamics (DFT-MD) simulations corroborated by continuous-flow experiments on Ni-SSZ-24 zeolite. The results of the simulations at operating conditions show that ethene molecules reversibly mobilize the active site as they exchange with the zeolite as ligands on Ni during reaction. Microkinetic modeling was conducted on the basis of free-energy profiles derived with DFT-MD for oligomerization on these mobile [(ethene)₂-Ni-alkyl]⁺ species. The model reproduces the experimentally observed high selectivity to dimerization and indicates that the mechanism is consistent with the observed second-order rate dependence on ethene pressure.

In the version of this Article originally published, on the right side of Fig. 4b, the ‘Aromatic cycle’ label was erroneously shifted outside of the central circular arrow into a position on part of the reaction cycle. This has been corrected in the online versions of the Article.

The combination of well-defined acid sites, shape-selective properties and outstanding stability places zeolites among the most practically relevant heterogeneous catalysts. The development of structure–performance descriptors for processes that they catalyse has been a matter of intense debate, both in industry and academia, and the direct conversion of methanol to olefins is a prototypical system in which various catalytic functions contribute to the overall performance. Propylene selectivity and resistance to coking are the two most important parameters in developing new methanol-to-olefin catalysts. Here, we present a systematic investigation on the effect of acidity on the performance of the zeolite ‘ZSM-5’ for the production of propylene. Our results demonstrate that the isolation of Brønsted acid sites is key to the selective formation of propylene. Also, the introduction of Lewis acid sites prevents the formation of coke, hence drastically increasing catalyst lifetime.

UiO-66 is a showcase example of an extremely stable metal–organic framework, which maintains its structural integrity during activation processes such as linker exchange and dehydration. The framework can even accommodate a substantial number of defects without compromising its stability. These observations point to an intrinsic dynamic flexibility of the framework, related to changes in the coordination number of the zirconium atoms. Herein we follow the dynamics of the framework in situ, by means of enhanced sampling molecular dynamics simulations such as umbrella sampling, during an activation process, where the coordination number of the bridging hydroxyl groups capped in the inorganic Zr₆(μ₃-O)₄(μ₃-OH)₄ brick is reduced from three to one. Such a reduction in the coordination number occurs during the dehydration process and in other processes where defects are formed. We observe a remarkable fast response of the system upon structural changes of the hydroxyl group. Internal deformation modes are detected, which point to linker decoordination and recoordination. Detached linkers may be stabilized by hydrogen bonds with hydroxyl groups of the inorganic brick, which gives evidence for an intrinsic dynamic acidity even in the absence of protic guest molecules. Our observations yield a major step forward in the understanding on the molecular level of activation processes realized experimentally but that is hard to track on a purely experimental basis.