Toward the Atomic Scale Simulation of Intricate Acidic Aluminosilicate Catalysts

Water structures on acidic zeolites and their roles in catalysis

The local reaction environment of catalytic active sites can be manipulated to modify the kinetics and thermodynamic properties of heterogeneous catalysis. Because of the unique physical-chemical nature of water, heterogeneously catalyzed reactions involving specific interactions between water molecules and active sites on catalysts exhibit distinct outcomes that are different from those performed in the absence of water. Zeolitic materials are being applied with the presence of water for heterogeneous catalytic reactions in the chemical industry and our transition to sustainable energy. Mechanistic investigation and in-depth understanding about the behaviors and the roles of water are essentially required for zeolite chemistry and catalysis. In this review, we focus on the discussions of the nature and structures of water adsorbed/stabilized on Brønsted and Lewis acidic zeolites based on experimental observations as well as theoretical calculation results. The unveiled functions of water structures in determining the catalytic efficacy of zeolite-catalyzed reactions have been overviewed and the strategies frequently developed for enhancing the stabilization of zeolite catalysts are highlighted. Recent advancement will contribute to the development of innovative catalytic reactions and the rationalization of catalytic performances in terms of activity, selectivity and stability with the presence of water vapor or in condensed aqueous phase.

Are the Brønsted acid sites in amorphous silica-alumina bridging?

Competing models exist to explain the differences in the activity of zeolites and amorphous silica-aluminas. Some postulate that silica-alumina contains dilute zeolitic bridging acid sites, while others favor a pseudo-bridging silanol model. We employed a selective isotope labeling strategy to assess the existence of Si-O(H)-Al bonds using NMR-based distance measurements.

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.

Challenges in modelling dynamic processes in realistic nanostructured materials at operating conditions

The question is addressed in how far current modelling strategies are capable of modelling dynamic phenomena in realistic nanostructured materials at operating conditions. Nanostructured materials used in applications are far from perfect; they possess a broad range of heterogeneities in space and time extending over several orders of magnitude. Spatial heterogeneities from the subnanometre to the micrometre scale in crystal particles with a finite size and specific morphology, impact the material's dynamics. Furthermore, the material's functional behaviour is largely determined by the operating conditions. Currently, there exists a huge length-time scale gap between attainable theoretical length-time scales and experimentally relevant scales. Within this perspective, three key challenges are highlighted within the molecular modelling chain to bridge this length-time scale gap. Methods are needed that enable (i) building structural models for realistic crystal particles having mesoscale dimensions with isolated defects, correlated nanoregions, mesoporosity, internal and external surfaces; (ii) the evaluation of interatomic forces with quantum mechanical accuracy albeit at much lower computational cost than the currently used density functional theory methods and (iii) derivation of the kinetics of phenomena taking place in a multi-length-time scale window to obtain an overall view of the dynamics of the process. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

Site-Averaged Ab Initio Kinetics: Importance Learning for Multistep Reactions on Amorphous Supports

Single-atom centers on amorphous supports include catalysts for polymerization, partial oxidation, metathesis, hydrogenolysis, and more. The disordered environment makes each site different, and the kinetics exponentially magnifies these differences to make ab initio site-averaged kinetics calculations extremely difficult. This work extends the importance learning algorithm for efficient and precise site-averaged kinetics estimates to ab initio calculations and multistep reaction mechanisms. Specifically, we calculate site-averaged proton transfer relaxation rates on an ensemble of cluster models representing Brønsted acid sites on silica-alumina. We include direct and water-assisted proton transfer pathways and simultaneously estimate the water adsorption and activation enthalpies for forward and backward proton transfers. We use density functional theory (DFT) to obtain a site-averaged rate, somewhat like a turnover frequency, for the proton transfer relaxation rate. Finally, we show that importance learning can provide orders-of-magnitude acceleration over standard sampling methods for site-averaged rate calculations in cases where the rate is dominated by a few highly active sites.

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.

A computational investigation of the decomposition of acetic acid in H-SSZ-13 and its role in the initiation of the MTO process

The zeolite-catalyzed reaction of acetic acid is important in the direct utilization of biomass and also plays a role in the reactivity of oxygenates in the methanol-to-olefins (MTO) process.

How Water and Ion Mobility Affect the NMR Fingerprints of the Hydrated JBW Zeolite: A Combined Computational-Experimental Investigation

An important aspect within zeolite synthesis is to make fully tunable framework materials with controlled aluminium distribution. A major challenge in characterising these zeolites at operating conditions is the presence of water. In this work, we investigate the effect of hydration on the 27 Al NMR parameters of the ultracrystalline K,Na-compensated aluminosilicate JBW zeolite using experimental and computational techniques. The JBW framework, with Si/Al ratio of 1, is an ideal benchmark system as a stepping stone towards more complicated zeolites. The presence and mobility of water and extraframework species directly affect NMR fingerprints. Excellent agreement between theoretical and experimental spectra is obtained provided dynamic methods are employed with hydrated structural models. This work shows how NMR is instrumental in characterising aluminium distributions in zeolites at operating conditions.

Exploring the In Situ Formation Mechanism of Polymeric Aluminum Chloride-Silica Gel Composites under Mechanical Grinding Conditions: As a High-Performance Nanocatalyst for the Synthesis of Xanthene and Pyrimidinone Compounds

The use of mechanical ball milling to facilitate the synthesis of organic compounds has attracted intense interest from organic chemists. Herein, we report a new process for the preparation of xanthene and pyrimidinone compounds by a one-pot method using polymeric aluminum chloride (PAC), silica gel, and reaction raw materials under mechanical grinding conditions. During the grinding process, polymeric aluminum chloride and silica gel were reconstituted in situ to obtain a new composite catalyst (PAC-silica gel). This catalyst has good stability (six cycles) and wide applicability (22 substrates). The Al-O-Si active center formed by in situ grinding recombination was revealed to be the key to the effective catalytic performance of the PAC-silica gel composites by the comprehensive analysis of the catalytic materials before and after use. In addition, the mechanism of action of the catalyst was verified using density functional theory, and the synthetic pathway of the xanthene compound was reasonably speculated with the experimental data. Mechanical ball milling serves two purposes in this process: not only to induce the self-assembly of silica and PAC into new composites but also to act as a driving force for the catalytic reaction to take place. From a practical point of view, this "one-pot" catalytic method eliminates the need for a complex preparation process for catalytic materials. This is a successful example of the application of mechanochemistry in materials and organic synthesis, offering unlimited possibilities for the application of inorganic polymer materials in green synthesis and catalysis promoted by mechanochemistry.

Machine learning potential era of zeolite simulation

Zeolites, owing to their great variety and complexity in structure and wide applications in chemistry, have long been the hot topic in chemical research. This perspective first presents a short retrospect of theoretical investigations on zeolites using the tools from classical force fields to quantum mechanics calculations and to the latest machine learning (ML) potential simulations. ML potentials as the next-generation technique for atomic simulation open new avenues to simulate and interpret zeolite systems and thus hold great promise for finally predicting the structure-functionality relation of zeolites. Recent advances using ML potentials are then summarized from two main aspects: the origin of zeolite stability and the mechanism of zeolite-related catalytic reactions. We also discussed the possible scenarios of ML potential application aiming to provide instantaneous and easy access of zeolite properties. These advanced applications could now be accomplished by combining cloud-computing-based techniques with ML potential-based atomic simulations. The future development of ML potentials for zeolites in the respects of improving the calculation accuracy, expanding the application scope and constructing the zeolite-related datasets is finally outlooked.

Gas Adsorption in Zeolite and Thin Zeolite Layers: Molecular Simulation, Experiment, and Adsorption Potential Theory

Molecular simulations and experiments are used to investigate methane adsorption in bulk and thin layers of MFI zeolite (silicalite-1). After comparing the theoretical adsorption data obtained using Grand Canonical Monte Carlo simulations for bulk MFI at various temperatures against experiments, zeolite layers with different crystalline orientations and levels of surface flexibility are considered. The data obtained for such prototypical systems allow us to rationalize both the qualitative and quantitative impact of external surface in nanoporous solids. In particular, due to strong confinement in zeolite pores, methane is found to adsorb at low pressures in the core of the zeolite while external surface adsorption occurs at pressures where the internal porosity of zeolite is saturated. Using Polanyi's adsorption potential theory, which is derived here from Hill's general scheme for adsorption, we provide a simple thermodynamic formalism to predict consistently adsorption both in the internal porosity and at the external surface of nanoporous solids. While this seminal theory has been already applied for gases in nanoporous solids, its extension to describe both surface and volume adsorption is important to provide a general rational framework for fluid adsorption in finely divided materials. We also discuss the applicability of this formalism for gas adsorption data under supercritical conditions.

Mechanistic Characterization of Zeolite-Catalyzed Aromatic Electrophilic Substitution at Realistic Operating Conditions

Zeolite-catalyzed benzene ethylation is an important industrial reaction, as it is the first step in the production of styrene for polymer manufacturing. Furthermore, it is a prototypical example of aromatic electrophilic substitution, a key reaction in the synthesis of many bulk and fine chemicals. Despite extensive research, the reaction mechanism and the nature of elusive intermediates at realistic operating conditions is not properly understood. More in detail, the existence of the elusive arenium ion (better known as Wheland complex) formed upon electrophilic attack on the aromatic ring is still a matter of debate. Temperature effects and the presence of protic guest molecules such as water are expected to impact the reaction mechanism and lifetime of the reaction intermediates. Herein, we used enhanced sampling ab initio molecular dynamics simulations to investigate the complete mechanism of benzene ethylation with ethene and ethanol in the H-ZSM-5 zeolite. We show that both the stepwise and concerted mechanisms are active at reaction conditions and that the Wheland intermediate spontaneously appears as a shallow minimum in the free energy surface after the electrophilic attack on the benzene ring. Addition of water enhances the protonation kinetics by about 1 order of magnitude at coverages of one water molecule per Brønsted acidic site. In the fully solvated regime, an overstabilization of the BAS as hydronium ion occurs and the rate enhancement disappears. The obtained results give critical atomistic insights in the role of water to selectively tune the kinetics of protonation reactions in zeolites.

A new mechanistic proposal for the aromatic cycle of the MTO process based on a computational investigation for H-SSZ-13

The paring mechanism of the aromatic cycle of the hydrocarbon pool is reinvestigated based on the heptamethylbenzenium cation adsorbed within H-SSZ-13 using quantum chemical calculations.

Revealing Brønsted Acidic Bridging SiOHAl Groups on Amorphous Silica-Alumina by Ultrahigh Field Solid-State NMR

Amorphous silica-aluminas (ASAs) are important acidic catalysts and supports for many industrially essential and sustainable processes. The identification of surface acid sites with their local structures on ASAs is of critical importance for tuning their catalytic properties but still remains a great challenge and is under debate. Here, ultrahigh magnetic field (35.2 T) ²⁷Al-{¹H} D-HMQC (dipolar-mediated heteronuclear multiple-quantum correlation) two-dimensional NMR experiments demonstrate two types of Brønsted acid sites in ASA catalysts. In addition to the known pseudobridging silanol acid sites, the use of ultrahigh field NMR provides the first direct experimental evidence for the existence of bridging silanol (BS: SiOHAl) acid sites in ASAs, which has been hotly debated in the past few decades. This discovery provides new opportunities for scientists and engineers to develop and apply ASAs in various reaction processes due to the significance of BS in chemical and fuel productions based on its strong Brønsted acidity.

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.

Engineering the Distinct Structure Interface of Subnano-alumina Domains on Silica for Acidic Amorphous Silica-Alumina toward Biorefining

Amorphous silica-aluminas (ASAs) are important solid catalysts and supports for many industrially essential and sustainable processes, such as hydrocarbon transformation and biorefining. However, the wide distribution of acid strength on ASAs often results in undesired side reactions, lowering the product selectivity. Here we developed a strategy for the synthesis of a unique class of ASAs with unvarying strength of Brønsted acid sites (BAS) and Lewis acid sites (LAS) using double-flame-spray pyrolysis. Structural characterization using high-resolution transmission electron microscopy (TEM) and solid-state nuclear magnetic resonance (NMR) spectroscopy showed that the uniform acidity is due to a distinct nanostructure, characterized by a uniform interface of silica-alumina and homogeneously dispersed alumina domains. The BAS population density of as-prepared ASAs is up to 6 times higher than that obtained by classical methods. The BAS/LAS ratio, as well as the population densities of BAS and LAS of these ASAs, could be tuned in a broad range. In cyclohexanol dehydration, the uniform Brønsted acid strength provides a high selectivity to cyclohexene and a nearly linear correlation between acid site densities and cyclohexanol conversion. Moreover, the concerted action of these BAS and LAS leads to an excellent bifunctional Brønsted-Lewis acid catalyst for glucose dehydration, affording a superior 5-hydroxymethylfurfural yield.

Zeolites at the Molecular Level: What Can Be Learned from Molecular Modeling

This review puts the development of molecular modeling methods in the context of their applications to zeolitic active sites. We attempt to highlight the utmost necessity of close cooperation between theory and experiment, resulting both in advances in computational methods and in progress in experimental techniques.

Effect of Impurities on the Initiation of the Methanol-to-Olefins Process: Kinetic Modeling Based on Ab Initio Rate Constants

Abstract

The relevance of a selection of organic impurities for the initiation of the MTO process was quantified in a kinetic model comprising 107 elementary steps with ab initio computed reaction barriers (MP2:DFT). This model includes a representative part of the autocatalytic olefin cycle as well as a direct initiation mechanism starting from methanol through CO-mediated direct C–C bond formation. We find that the effect of different impurities on the olefin evolution varies with the type of impurity and their partial pressures. The reactivity of the considered impurities for initiating the olefin cycle increases in the order formaldehyde < di-methoxy methane < CO < methyl acetate < ethanol < ethene < propene. In our kinetic model, already extremely low quantities of impurities such as ethanol lead to faster initiation than through direct C–C bond formation which only matters in complete absence of impurities.

Graphic Abstract

Shape selectivity effects in the hydroconversion of perhydrophenanthrene over bifunctional catalysts

The zeolite pore structure dictates the formation of isomers, which in turn influences the preferred ring opening products and the distribution of cracking products.

Theoretical investigation of the side-chain mechanism of the MTO process over H-SSZ-13 using DFT and ab initio calculations

The side-chain mechanism of the methanol-to-olefins process over the H-SSZ-13 acidic zeolite was investigated using periodic density functional theory with corrections from highly accurate ab intio calculations on large cluster models.