Computational elucidation of the transition state shape selectivity phenomenon

Tailoring parameters for QM/MM simulations: accurate modeling of adsorption and catalysis in zirconium-based metal-organic frameworks

Quantum mechanics/molecular mechanics (QM/MM) simulations offer an efficient way to model reactions occurring in complex environments. This study introduces a specialized set of charge and Lennard-Jones parameters tailored for electrostatically embedded QM/MM calculations, aiming to accurately model both adsorption processes and catalytic reactions in zirconium-based metal-organic frameworks (Zr-MOFs). To validate our approach, we compare adsorption energies derived from QM/MM simulations against experimental results and Monte Carlo simulation outcomes. The developed parameters showcase the ability of QM/MM simulations to represent long-range electrostatic and van der Waals interactions faithfully. This capability is evidenced by the prediction of adsorption energies with a low root mean square error of 1.1 kcal mol-1 across a wide range of adsorbates. The practical applicability of our QM/MM model is further illustrated through the study of glucose isomerization and epimerization reactions catalyzed by two structurally distinct Zr-MOF catalysts, UiO-66 and MOF-808. Our QM/MM calculations closely align with experimental activation energies. Importantly, the parameter set introduced here is compatible with the widely used universal force field (UFF). Moreover, we thoroughly explore how the size of the cluster model and the choice of density functional theory (DFT) methodologies influence the simulation outcomes. This work provides an accurate and computationally efficient framework for modeling complex catalytic reactions within Zr-MOFs, contributing valuable insights into their mechanistic behaviors and facilitating further advancements in this dynamic area of research.

Molecular Views on Mechanisms of Brønsted Acid-Catalyzed Reactions in Zeolites

The Brønsted acidity of proton-exchanged zeolites has historically led to the most impactful applications of these materials in heterogeneous catalysis, mainly in the fields of transformations of hydrocarbons and oxygenates. Unravelling the mechanisms at the atomic scale of these transformations has been the object of tremendous efforts in the last decades. Such investigations have extended our fundamental knowledge about the respective roles of acidity and confinement in the catalytic properties of proton exchanged zeolites. The emerging concepts are of general relevance at the crossroad of heterogeneous catalysis and molecular chemistry. In the present review, emphasis is given to molecular views on the mechanism of generic transformations catalyzed by Brønsted acid sites of zeolites, combining the information gained from advanced kinetic analysis, in situ, and operando spectroscopies, and quantum chemistry calculations. After reviewing the current knowledge on the nature of the Brønsted acid sites themselves, and the key parameters in catalysis by zeolites, a focus is made on reactions undergone by alkenes, alkanes, aromatic molecules, alcohols, and polyhydroxy molecules. Elementary events of C-C, C-H, and C-O bond breaking and formation are at the core of these reactions. Outlooks are given to take up the future challenges in the field, aiming at getting ever more accurate views on these mechanisms, and as the ultimate goal, to provide rational tools for the design of improved zeolite-based Brønsted acid catalysts.

Carbocation chemistry confined in zeolites: spectroscopic and theoretical characterizations

Acid-catalyzed reactions inside zeolites are one type of broadly applied industrial reactions, where carbocations are the most common intermediates of these reaction processes, including methanol to olefins, alkene/aromatic alkylation, and hydrocarbon cracking/isomerization. The fundamental research on these acid-catalyzed reactions is focused on the stability, evolution, and lifetime of carbocations under the zeolite confinement effect, which greatly affects the efficiency, selectivity and deactivation of zeolite catalysts. Therefore, a profound understanding of the carbocations confined in zeolites is not only beneficial to explain the reaction mechanism but also drive the design of new zeolite catalysts with ideal acidity and cages/channels. In this review, we provide both an in-depth understanding of the stabilization of carbocations by the pore confinement effect and summary of the advanced characterization methods to capture carbocations in zeolites, including UV-vis spectroscopy, solid-state NMR, fluorescence microscopy, IR spectroscopy and Raman spectroscopy. Also, we clarify the relationship between the activity and stability of carbocations in zeolite-catalyzed reactions, and further highlight the role of carbocations in various hydrocarbon conversion reactions inside zeolites with diverse frameworks and varying acidic properties.

Design and Synthesis of the Active Site Environment in Zeolite Catalysts for Selectively Manipulating Mechanistic Pathways

By combining kinetics and theoretical calculations, we show here the benefits of going beyond the concept of static localized and defined active sites on solid catalysts, into a system that globally and dynamically considers the active site located in an environment that involves a scaffold structure particularly suited for a target reaction. We demonstrate that such a system is able to direct the reaction through a preferred mechanism when two of them are competing. This is illustrated here for an industrially relevant reaction, the diethylbenzene–benzene transalkylation. The zeolite catalyst (ITQ-27) optimizes location, density, and environment of acid sites to drive the reaction through the preselected and preferred diaryl-mediated mechanism, instead of the alkyl transfer pathway. This is achieved by minimizing the activation energy of the selected pathway through weak interactions, much in the way that it occurs in enzymatic catalysts. We show that ITQ-27 outperforms previously reported zeolites for the DEB-Bz transalkylation and, more specifically, industrially relevant zeolites such as faujasite, beta, and mordenite.

Insights into Reaction Kinetics in Confined Space: Real Time Observation of Water Formation under a Silica Cover

We offer a comprehensive approach to determine how physical confinement can affect the water formation reaction. By using free-standing crystalline SiO₂ bilayer supported on Ru(0001) as a model system, we studied the water formation reaction under confinement in situ and in real time. Low-energy electron microscopy reveals that the reaction proceeds via the formation of reaction fronts propagating across the Ru(0001) surface. The Arrhenius analyses of the front velocity yield apparent activation energies (E_a^app) of 0.32 eV for the confined and 0.59 eV for the nonconfined reaction. DFT simulations indicate that the rate-determining step remains unchanged upon confinement, therefore ruling out the widely accepted transition state effect. Additionally, H₂O accumulation cannot explain the change in E_a^app for the confined cases studied because its concentration remains low. Instead, numerical simulations of the proposed kinetic model suggest that the H₂ adsorption process plays a decisive role in reproducing the Arrhenius plots.

Synergistic construction of bifunctional and stable Pt/HZSM-5-based catalysts for efficient catalytic cracking of n-butane

Efficient conversion of light alkanes is of essential significance for enhancing the utilization efficiency of resources and exploring the activation and evolution regulation of C-C and C-H bonds in stable molecules. The processes are often executed with catalysts under harsh conditions. The olefin yield and metal stability have been the long-standing concerns. Herein, we report a facile strategy of constructing a bifunctional Pt/HZSM-5-based catalyst by two-step atomic layer deposition (ALD) to achieve a high light olefin formation rate of 0.48 mmol gcat-1·min-1 in the catalytic cracking of n-butane at 600 °C, which is ∼2.2 times higher than that of the conventional Pt/HZSM-5 catalyst (0.22 mmol gcat-1·min-1). Moreover, the bifunctional Pt/HZSM-5-based catalyst exhibited outstanding recyclability and excellent metal stability against sintering in comparison with conventional Pt/HZSM-5. Detailed microscopic and spectroscopic characterization studies demonstrate that the metal oxide (TiO2 or Al2O3) coating not only prevents the metal from high-temperature sintering, but also regulates the proportion of coordinately unsaturated platinum surface atoms. Theoretical calculations further confirm the preference of nucleation of TiO2 or Al2O3 on coordinately unsaturated platinum sites, which in turn modulates the bifunctional dehydrogenation-cracking pathway to improve the olefin formation rate.

Engineering of Transition Metal Catalysts Confined in Zeolites

Transition metal-zeolite composites are versatile catalytic materials for a wide range of industrial and lab-scale processes. Significant advances in fabrication and characterization of well-defined metal centers confined in zeolite matrixes have greatly expanded the library of available materials and, accordingly, their catalytic utility. In this review, we summarize recent developments in the field from the perspective of materials chemistry, focusing on synthesis, postsynthesis modification, (operando) spectroscopy characterization, and computational modeling of transition metal-zeolite catalysts.

Intraparticle Diffusional versus Site Effects on Reaction Pathways in Liquid-Phase Cross Aldol Reactions

Chemo- and regioselectivity in a heterogeneously catalyzed cross aldol reaction were directed by tuning the nature of the sites, textural properties, and reaction conditions. Catalysts included sulfonic acid-functionalized resins or SBA-15 with varying particle size or pore diameter, H-BEA zeolites, and Sn-BEA zeotype; conditions were 25 °C to 170 °C in organic media. Benzaldehyde and 2-butanone yielded branched (reaction at -CH₂ - of butanone) and linear (reaction at -CH₃ ) addition and condensation products; and fission of the branched aldol led to β-methyl styrene and acetic acid. Strong acids promoted the dehydration step, and regioselectivity originated from preferred formation of the branched aldol. Both, resins and functionalized SBA-15 materials yielded predominantly the branched condensation product, unless particle morphology or temperature moved the reaction into the diffusion-limited regime, in which case more fission products were formed, corresponding to Wheeler Type II selectivity. For H-form zeolites, fission of the branched aldol competed with dehydration of the linear aldol, possibly because weaker acidity or steric restrictions prevented dehydration of the branched aldol.

Acidic Properties and Structure-Activity Correlations of Solid Acid Catalysts Revealed by Solid-State NMR Spectroscopy

Solid acid materials with tunable structural and acidic properties are promising heterogeneous catalysts for manipulating and/or emulating the activity and selectivity of industrially important catalytic reactions. On the other hand, the performances of acid-catalyzed reactions are mostly dictated by the acidic features, namely, type (Brønsted vs Lewis acidity), amount, strength, and local environment of acid sites. The latter is relevant to their location (intra- vs extracrystalline), and possible confinement and Brønsted-Lewis acid synergy effects that may strongly affect the host-guest interactions, reaction mechanism, and shape selectivity of the catalytic system. This account aims to highlight some important applications of state-of-the-art solid-state NMR (SSNMR) techniques for exploring the structural and acidic properties of solid acid catalysts as well as their catalytic performances and relevant reaction pathway invoked. In addition, density functional theory (DFT) calculations may be exploited in conjunction with experimental SSNMR studies to verify the structure-activity correlations of the catalytic system at a microscopic scale. We describe in this Account the developments and applications of advanced ex situ and/or in situ SSNMR techniques, such as two-dimensional (2D) double-quantum magic-angle spinning (DQ MAS) homonuclear correlation spectroscopy for structural investigation of solid acids as well as study of their acidic properties. Moreover, the energies and electronic structures of the catalysts and detailed catalytic reaction processes, including the identification of reaction species, elucidation of reaction mechanism, and verification of structure-activity correlations, made available by DFT theoretical calculations were also discussed. Relevant discussions will focus primarily on results obtained from our laboratories in the past decade, including (i) quantitative and qualitative acidity characterization utilizing assorted probe molecules, (ii) probing the spatial proximity and synergy effect of acid sites, and (iii) influence of acid features and pore confinement effect on catalytic activity, transition-state stability, reaction pathway, and product selectivity of solid acid catalysts such as zeolites, metal oxides, and heteropolyacids. It is conclusive that a synergy of acidity (local effect) and pore confinement (environmental effect) tend to strongly dictate the formations of intermediates and transition states, hence, the reaction pathways and catalytic performance of solid acid catalysts. We hope that these information can provide additional insights toward our understanding in heterogeneous catalysis, especially the roles of structural and acidic properties on catalytic performances and reaction mechanism of acid-catalyzed systems, which should be beneficial for rational design of solid acid catalysts.

Quantum chemical ab initio prediction of proton exchange barriers between CH4 and different H-zeolites

A hybrid MP2:DFT (second-order Møller-Plesset perturbation theory-density functional theory) method that combines MP2 calculations for cluster models with DFT calculations for the full periodic structure is used to localize minima and transition structures for proton jumps at different Brønsted sites in different frameworks (chabazite, faujasite, ferrierite, and ZSM-5) and at different crystallographic positions of a given framework. The MP2 limit for the periodic structures is obtained by extrapolating the results of a series of cluster models of increasing size. A coupled-cluster (CCSD(T)) correction to MP2 energies is calculated for cluster models consisting of three tetrahedra. For the adsorption energies, this difference is small, between 0.1 and 0.9 kJ/mol, but for the intrinsic proton exchange barriers, this difference makes a significant (10.85 ± 0.25 kJ/mol) and almost constant contribution across different systems. The total values of the adsorption energies vary between 22 and 34 kJ/mol, whereas the total proton exchange energy barriers fall in the narrow range of 152-156 kJ/mol. After adding nuclear motion contributions (harmonic approximation, 298 K), intrinsic enthalpy barriers between 134 and 141 kJ/mol and apparent energy barriers between 105 and 118 kJ/mol are predicted for the different sites examined for the different frameworks. These predictions are consistent with experimental results available for faujasite, ferrierite, and ZSM-5.

Advances in theory and their application within the field of zeolite chemistry

Zeolites are versatile and fascinating materials which are vital for a wide range of industries, due to their unique structural and chemical properties, which are the basis of applications in gas separation, ion exchange and catalysis. Given their economic impact, there is a powerful incentive for smart design of new materials with enhanced functionalities to obtain the best material for a given application. Over the last decades, theoretical modeling has matured to a level that model guided design has become within reach. Major hurdles have been overcome to reach this point and almost all contemporary methods in computational materials chemistry are actively used in the field of modeling zeolite chemistry and applications. Integration of complementary modeling approaches is necessary to obtain reliable predictions and rationalizations from theory. A close synergy between experimentalists and theoreticians has led to a deep understanding of the complexity of the system at hand, but also allowed the identification of shortcomings in current theoretical approaches. Inspired by the importance of zeolite characterization which can now be performed at the single atom and single molecule level from experiment, computational spectroscopy has grown in importance in the last decade. In this review most of the currently available modeling tools are introduced and illustrated on the most challenging problems in zeolite science. Directions for future model developments will be given.

Solid-state nuclear magnetic resonance investigations of the nature, property, and activity of acid sites on solid catalysts

Further progress in the field of heterogeneous catalysis depends on our knowledge of the nature and behavior of surface sites on solid catalysts and of the mechanisms of chemical reactions catalyzed by these materials. In the past decades, solid-state NMR spectroscopy has been developed to an important tool for routine characterization of solid catalysts. The present work gives a review on experimental approaches and applications of solid-state NMR spectroscopy for investigating Brønsted and Lewis sites on solid acids. Studies focusing on the generation of surface sites via post-synthesis modification routes of microporous and mesoporous materials support the development of new and the improvement of existing catalyst systems. High-temperature and flow techniques of in situ solid-state NMR spectroscopy allow a deeper insight into the mechanisms of heterogeneously catalyzed reactions and open the way for studying the activity of acidic surface sites. They help to clarify the activation of reactants on Brønsted and Lewis acid sites and improve our understanding of mechanisms affecting the selectivity of acid-catalyzed reactions.

Electronic view on ethene adsorption in Cu(I) exchanged zeolites

Ethene adsorption on isolated Cu(i) sites in two types of zeolites (faujasite and MFI) is investigated by means of the embedded cluster method. Structures, energetic stabilities and C[double bond, length as m-dash]C stretching vibrations in adsorption complexes are discussed. Furthermore, for interpretative purposes, the interaction energies are decomposed, using novel approaches based on so called natural orbitals for chemical valence. Ethene is always symmetrically bound to Cu(i) ion by both C atoms. In some cases two local minima of similar stability on the potential energy surface, differing by Cu(i) site relaxation can be found that may be simultaneously populated in equilibrium. Binding energies usually decrease with the degree of reconstruction of Cu(i) site after adsorption, however, in particular cases, a more distorted structure can be slightly more stable if favorable pi* back donation overwhelms the distortion effects. Calculated values of binding energies for Cu(i)-Y zeolite (about 80 kJ mol(-1)) agree well with microcalorimetric data. We predict that ethene binding in MFI is over two times stronger (to the best of our knowledge no experimental data are available). The C[double bond, length as m-dash]C stretching frequency is not site specific, but depends only on the type of copper connectivity to oxygen nodes. The appearance of two C[double bond, length as m-dash]C bands in IR spectra of Cu(i)-faujasite can be explained as the effect of coexistence of two types of adsorption complexes, with Cu(i) coordinated to one or two framework tetrahedrons, respectively. In Cu(i)-MFI, only one type of adsorption complex with Cu(i) ion coordinated to a single tetrahedron exists, as only a single C[double bond, length as m-dash]C band is present in IR spectra.

ONIOM study of isomerization reactions of aromatic hydrocarbons in acidic mordenite zeolite

Two different mechanisms: Shift 1,2-isomerization and isomerization via the disproportionation reaction are investigated for aromatic hydrocarbons over acidic mordenite zeolite by using our own n-layered integrated molecular orbital and molecular mechanics (ONIOM) scheme. The picture shows a schematic energy profile for the isomerization of toluene catalyzed by acidic mordenite.The geometrical structures and energies of isomerization reactions of aromatic hydrocarbons catalyzed by the 128T cluster model of acidic mordenite zeolite are studied by using our own n-layered integrated molecular orbital and molecular mechanics (ONIOM) scheme. The aromatic hydrocarbons investigated are toluene, xylenes, and trimethylbenzenes. Two proposed intramolecular reaction mechanisms, that is, shift 1,2-isomerization and isomerization via the disproportionation reaction, are studied and discussed in detail. Both of the mechanisms firstly proceed by protonation of aromatic hydrocarbons to form the aromatic-based carbenium ions. Attempts to determine the prevailing reaction mechanism of the isomerization are evaluated. In most cases, the isomerization reactivity sequence is: trimethylbenzenes>xylenes>toluene. The conclusions obtained using the ONIOM scheme agree with the available periodic DFT and experimental results on acidic mordenite zeolite.

A link between reactivity and local structure in acid catalysis on zeolites

The extent to which spatial constraints influence rates and pathways in catalysis depends on the structure of intermediates, transition states, and active sites involved. We aim to answer, as we seek insights into catalytic mechanisms and site requirements, persistent questions about the potential for controlling rates and selectivities by rational design of spatial constraints around active sites within inorganic structures useful as catalysts. This Account addresses these matters for the specific case of reactions on zeolites that contain Brønsted acid sites encapsulated within subnanometer channels. We compare and contrast here the effects of local zeolite structure on the dynamics of the carbonylation of surface methyl groups and of the isotopic exchange of CD4 with surface OH groups on zeolites. Methyl and hydroxyl groups are the smallest monovalent cations relevant in catalysis by zeolites. Their small size, taken together with their inability to desorb except via reactions with other species, allowed us to discriminate between stabilization of cationic transition states and stabilization of adsorbed reactants and products by spatial constraints. We show that apparent effects of proton density and of zeolite channel structure on dimethyl ether carbonylation turnover rates reflect instead the remarkable specificity of eight-membered ring zeolite channels in accelerating kinetically relevant steps that form *COCH3 species via CO insertion into methyl groups. This specificity reflects the selective stabilization of cationic transition states via interactions with framework oxygen anions. These findings for carbonylation catalysts contrast sharply the weak effects of channel structure on the rate of exchange of CD4 with OH groups. This latter reaction involves concerted symmetric transition states with much lower charge than that required for CH3 carbonylation. Our Account extends the scope of shape selectivity concepts beyond those reflecting size exclusion and preferential adsorption. Our ability to discriminate among various effects of spatial constraints depends critically on dissecting chemical conversions into elementary steps of kinetic relevance and on eliminating secondary reactions and accounting for the concomitant effects of zeolite structure on the stability of adsorbed reactants and intermediates.

Specificity of Sites within Eight-Membered Ring Zeolite Channels for Carbonylation of Methyls to Acetyls

The acid-catalyzed formation of carbon-carbon bonds from C1 precursors via CO insertion into chemisorbed methyl groups occurs selectively within eight-membered ring (8-MR) zeolite channels. This elementary step controls catalytic carbonylation rates of dimethyl ether (DME) to methyl acetate. The number of O-H groups within 8-MR channels was measured by rigorous deconvolution of the infrared bands for O-H groups in cation-exchanged and acid forms of mordenite (M,H-MOR) and ferrierite (H-FER) after adsorption of basic probe molecules of varying size. DME carbonylation rates are proportional to the number of O-H groups within 8-MR channels. Na+ cations selectively replaced protons within 8-MR channels and led to a disproportionate decrease in carbonylation turnover rates (per total H+). These conclusions are consistent with the low or undetectable rates of carbonylation on zeolites without 8-MR channels (H-BEA, H-FAU, H-MFI). Such specificity of methyl reactivity upon confinement within small channels appears to be unprecedented in catalysis by microporous solids, which typically select reactions by size exclusion of bulkier transition states.