Experimental and Theoretical Evidence for the Promotional Effect of Acid Sites on the Diffusion of Alkenes through Small-Pore Zeolites.

The diffusion of saturated and unsaturated hydrocarbons is of fundamental importance for many zeolite-catalyzed processes. Transport of small alkenes in the confined zeolite pores can become hindered, resulting in a significant impact on the ultimate product selectivity and separation. Herein, intracrystalline light olefin/paraffin diffusion through the 8-ring windows of zeolite SAPO-34 is characterized by a complementary set of first-principle molecular dynamics simulations, PFG-NMR experiments, and pulse-response temporal analysis of products measurements, yielding information at different length and time scales. Our results clearly show a promotional effect of the presence of Brønsted acid sites on the diffusion rate of ethene and propene, whereas transport of alkanes is found to be insensitive to the presence of acid sites. The enhanced diffusivity of unsaturated hydrocarbons is ascribed to the formation of favorable π-H interactions with acid protons, as confirmed by IR spectroscopy measurements. The acid site distribution is proven to be an important design parameter for optimizing product distributions and separations.

Light Olefin Diffusion during the MTO Process on H-SAPO-34: A Complex Interplay of Molecular Factors.

The methanol-to-olefins process over H-SAPO-34 is characterized by its high shape selectivity toward light olefins. The catalyst is a supramolecular system consisting of nanometer-sized inorganic cages, decorated by Brønsted acid sites, in which organic compounds, mostly methylated benzene species, are trapped. These hydrocarbon pool species are essential to catalyze the methanol conversion but may also clog the pores. As such, diffusion of ethene and propene plays an essential role in determining the ultimate product selectivity. Enhanced sampling molecular dynamics simulations based on either force fields or density functional theory are used to determine how molecular factors influence the diffusion of light olefins through the 8-ring windows of H-SAPO-34. Our simulations show that diffusion through the 8-ring in general is a hindered process, corresponding to a hopping event of the diffusing molecule between neighboring cages. The loading of different methanol, alkene, and aromatic species in the cages may substantially slow down or facilitate the diffusion process. The presence of Brønsted acid sites in the 8-ring enhances the diffusion process due to the formation of a favorable π-complex host-guest interaction. Aromatic hydrocarbon pool species severely hinder the diffusion and their spatial distribution in the zeolite crystal may have a significant effect on the product selectivity. Herein, we unveil how molecular factors influence the diffusion of light olefins in a complex environment with confined hydrocarbon pool species, high olefin loadings, and the presence of acid sites by means of enhanced molecular dynamics simulations under operating conditions.

Operando Modeling of Zeolite-Catalyzed Reactions Using First-Principles Molecular Dynamics Simulations.

Within this Perspective, we critically reflect on the role of first-principles molecular dynamics (MD) simulations in unraveling the catalytic function within zeolites under operating conditions. First-principles MD simulations refer to methods where the dynamics of the nuclei is followed in time by integrating the Newtonian equations of motion on a potential energy surface that is determined by solving the quantum-mechanical many-body problem for the electrons. Catalytic solids used in industrial applications show an intriguing high degree of complexity, with phenomena taking place at a broad range of length and time scales. Additionally, the state and function of a catalyst critically depend on the operating conditions, such as temperature, moisture, presence of water, etc. Herein we show by means of a series of exemplary cases how first-principles MD simulations are instrumental to unravel the catalyst complexity at the molecular scale. Examples show how the nature of reactive species at higher catalytic temperatures may drastically change compared to species at lower temperatures and how the nature of active sites may dynamically change upon exposure to water. To simulate rare events, first-principles MD simulations need to be used in combination with enhanced sampling techniques to efficiently sample low-probability regions of phase space. Using these techniques, it is shown how competitive pathways at operating conditions can be discovered and how broad transition state regions can be explored. Interestingly, such simulations can also be used to study hindered diffusion under operating conditions. The cases shown clearly illustrate how first-principles MD simulations reveal insights into the catalytic function at operating conditions, which could not be discovered using static or local approaches where only a few points are considered on the potential energy surface (PES). Despite these advantages, some major hurdles still exist to fully integrate first-principles MD methods in a standard computational catalytic workflow or to use the output of MD simulations as input for multiple length/time scale methods that aim to bridge to the reactor scale. First of all, methods are needed that allow us to evaluate the interatomic forces with quantum-mechanical accuracy, albeit at a much lower computational cost compared to currently used density functional theory (DFT) methods. The use of DFT limits the currently attainable length/time scales to hundreds of picoseconds and a few nanometers, which are much smaller than realistic catalyst particle dimensions and time scales encountered in the catalysis process. One solution could be to construct machine learning potentials (MLPs), where a numerical potential is derived from underlying quantum-mechanical data, which could be used in subsequent MD simulations. As such, much longer length and time scales could be reached; however, quite some research is still necessary to construct MLPs for the complex systems encountered in industrially used catalysts. Second, most currently used enhanced sampling techniques in catalysis make use of collective variables (CVs), which are mostly determined based on chemical intuition. To explore complex reactive networks with MD simulations, methods are needed that allow the automatic discovery of CVs or methods that do not rely on a priori definition of CVs. Recently, various data-driven methods have been proposed, which could be explored for complex catalytic systems. Lastly, first-principles MD methods are currently mostly used to investigate local reactive events. We hope that with the rise of data-driven methods and more efficient methods to describe the PES, first-principles MD methods will in the future also be able to describe longer length/time scale processes in catalysis. This might lead to a consistent dynamic description of all steps-diffusion, adsorption, and reaction-as they take place at the catalyst particle level.

Mobility and Reactivity of Cu+ Species in Cu-CHA Catalysts under NH3-SCR-NOx Reaction Conditions: Insights from AIMD Simulations.

The mobility of the copper cations acting as active sites for the selective catalytic reduction of nitrogen oxides with ammonia in Cu-CHA catalysts varies with temperature and feed composition. Herein, the migration of [Cu(NH3)2]+ complexes between two adjacent cavities of the chabazite structure, including other reactant molecules (NO, O2, H2O, and NH3), in the initial and final cavities is investigated using ab initio molecular dynamics (AIMD) simulations combined with enhanced sampling techniques to describe hopping events from one cage to the other. We find that such diffusion is only significantly hindered by the presence of excess NH3 or NO in the initial cavity, since both reactants form with [Cu(NH3)2]+ stable intermediates which are too bulky to cross the 8-ring windows connecting the cavities. The presence of O2 modifies strongly the interaction of NO with Cu+. At low temperatures, we observe NO detachment from Cu+ and increased mobility of the [Cu(NH3)2]+ complex, while at high temperatures, NO reacts spontaneously with O2 to form NO2. The present simulations give evidence for recent experimental observations, namely, an NH3 inhibition effect on the SCR reaction at low temperatures, and transport limitations of NO and NH3 at high temperatures. Our first principle simulations mimicking operating conditions support the existence of two different reaction mechanisms operating at low and high temperatures, the former involving dimeric Cu(NH3)2-O2-Cu(NH3)2 species and the latter occurring by direct NO oxidation to NO2 in one single cavity.

Theoretical and Spectroscopic Evidence of the Dynamic Nature of Copper Active Sites in Cu-CHA Catalysts under Selective Catalytic Reduction (NH3-SCR-NOx) Conditions.

The dynamic nature of the copper cations acting as active sites for selective catalytic reduction of nitrogen oxides with ammonia is investigated using a combined theoretical and spectroscopic approach. Ab initio molecular dynamics simulations of Cu-CHA catalysts in contact with reactants and intermediates at realistic operating conditions show that only ammonia is able to release Cu+ and Cu2+ cations from their positions coordinated to the zeolite framework, forming mobile Cu+(NH3)2 and Cu2+(NH3)4 complexes that migrate to the center of the cavity. Herein, we give evidence that such mobilization of copper cations modifies the vibrational fingerprint in the 800-1000 cm-1 region of the IR spectra. Bands associated with the lattice asymmetric T-O-T vibrations are perturbed by the presence of coordinated cations, and allow one to experimentally follow the dynamic reorganization of the active sites at operating conditions.

Mechanistic Insight into the Framework Methylation of H-ZSM-5 for Varying Methanol Loadings and Si/Al Ratios Using First-Principles Molecular Dynamics Simulations.

The methanol-to-hydrocarbon 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 densities 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 the temperature and acid site density. Enhanced sampling molecular dynamics simulations give evidence that the barrier for methanol conversion is significantly affected by the neighborhood of an additional acid site, suggesting that cooperative effects influence methanol clustering and reactivity. The insights obtained are important steps in optimizing the catalyst and engineering the induction period of the methanol-to-hydrocarbon process.

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 C4-C8 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 B2 (3° → 2°) are expected to be dominant at operating conditions, whereas modes B1 (2° → 3°), C (2° → 2°), D2 (2° → 1°), and E2 (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.