Abstract:This study introduces a novel kinetic Monte Carlo (KMC) simulation package which models H-ZSM-5 crystals across experimentally relevant time and length scales to understand the role of transport during arene interconversion reactions (~100 reactions). This small subset of the methanol-to-hydrocarbon (MTH) network was previously modeled using periodic, dispersion-corrected density functional theory (DFT) to determine activation barriers and reaction energies for these KMC methods. Transport of arene molecules through the straight and sinusoidal channels of MFI was modeled as site-hopping and the DFT-calculated barriers are incorporated into the KMC model to account for mass-transport limitations. Barriers of different arene molecules trend well with their effective radii, and species with a smaller effective radii diffusive more readily. A previously published maximum rate analysis of arene interconversion pathways—previously validated by experimental data—is compared to a diffusion-free KMC model to confirm the accuracy of this KMC package. The temperature and pressure dependencies of rates obtained from KMC agree well with those of maximum rate analysis on the diffusion-free model, demonstrating that KMC effectively predicts rates as well as maximum rate analysis methods commonly used in kinetic applications of DFT. Arene interconversion pathways were also analyzed on KMC models incorporating diffusion to and from interior crystal sites. These simulations suggest that large species, such as hexamethylbenzene, become trapped at 10–20% of sites, thus causing site deactivation by limiting diffusion through MFI channels and lowering overall rates of product formation. Benzene diffusion barriers are artificially varied from 20–200 kJ mol −1 and rates of benzene methylation decrease by 4-fold with diffusion barriers greater than 80 kJ mol −1 ; this suggests that species with diffusion barriers greater than 80 kJ mol −1 (such as penta- and hexamethylbenzene) will likely become trapped at interior sites and ultimately cause catalyst deactivation. This study serves as a proof-of-concept for a novel KMC package that expedites kinetic analysis of complex reaction pathways and introduces mass-transport limitations which are not commonly accounted for in kinetic DFT studies. This KMC package can predict the behavior of diffusion-limited species, such as penta- and hexamethylbenzene, and the mechanisms by which they are formed and eventually lead to catalyst deactivation.

Diffusion behaviors of binary mixtures of alkanes and aromatics through ZSM‐5 zeolite: A kinetic Monte Carlo study

Diffusion of Binary Mixtures in Zeolites: Kinetic Monte Carlo versus Molecular Dynamics Simulations

Kinetic Modeling of C2–C7 Olefins Interconversion over ZSM-5 Catalyst

Solvents Molecular Mobility in Coked Catalyst ZSM-5 Studied by NMR Relaxation and Pulsed Field Gradient Techniques

Diffusion of a hydrocarbon mixture in a one-dimensional zeolite channel: An exclusion model approach

Prediction of C6–C12 Interconversion Rates Using DFT and Zeolite-specific kinetic Monte Carlo Simulation Methods

Multiscale exploration of hydrocarbon adsorption and hopping through ZSM-5 channels - from Monte Carlo modelling to experiment

Molecular Simulation of Adsorption and Separation of Binary Mixtures for Short Linear Alkanes in Zeolite

The Adsorption and Localization of Mixtures of C4-C7 Alkane Isomers in Zeolites by Computer Simulation.

Mixtures of Ethane with CO$_2$ and Water Simulated in ZSM-22: The Role of Polarity and Hydrogen Bonding

Methanol diffusion and dynamics in zeolite H-ZSM-5 probed by quasi-elastic neutron scattering and classical molecular dynamics simulations

Adsorption and Separation of Ternary and Quaternary Mixtures of Short Linear Alkanes in Zeolites by Molecular Simulation

Molecular Insights into Adsorption and Diffusion Mechanism of N-Hexane in MFI Zeolites with Different Si-to-Al Ratios and Counterions

Experimental proof for resonant diffusion of normal alkanes in LTL and ZSM-12 zeolites

Untangling Framework Confinements: A Dynamical Study on Bulky Aromatic Molecules in MFI Zeolites

Modulation of B-Axis Thickness Within MFI Zeolite: Correlation with Variation of Product Diffusion and Coke Distribution in the Methanol-to-hydrocarbons Conversion

Insight into the effect of Si/Al ratio on the adsorption and diffusion behavior of butane isomers in HZSM-5 zeolites

Dynamics of Water, CO2, Ethane and Their Mixtures in ZSM-22 Zeolite: The Role of Polarity and Hydrogen Bonding

Probing Monomer and Dimer Adsorption Trends in the MFI Framework

Unraveling Molecular Diffusion Behaviors on the Acidic Sites of H‐ZSM‐5 Zeolite using Time‐Resolved in‐situ FT‐IR Technique

Diffusion/adsorption/dissolution coupling process of methane gas in kerogen nano-slit: Molecular simulation and a chemical potential-based mass transfer model