Effect of temperature and branching on the nature and stability of alkene cracking intermediates in H-ZSM-5

Detailed analysis of olefins and diolefins in hydrotreated plastic waste pyrolysis oils by GC-VUV

Recently, the hydrotreatment of plastic waste pyrolysis oils has been considered a promising route to remove heteroatoms and unsaturates from plastic pyrolysis oils to be used as steam cracking feedstock. As plastic waste pyrolysis oils differ vastly in composition from conventional crude oils, accurately identifying the different molecular species present in the mixture is essential to acquire further insights into their hydrogenation pathways. This work used gas chromatography coupled to a vacuum ultraviolet detector (GC-VUV) to obtain accurate quantification of (di)olefins compared to more conventional methods. It was found that there is no difference in intrinsic hydrogenation rates between diolefins and olefins, at least for chain lengths above C10. Moreover, an explicit dependency of the carbon chain length is identified on the hydrogenation rate for both olefins and diolefins with optimal hydrogenation rates up to C15. This is a consequence of the competition between hydrogenation and cracking reactions. Shorter unsaturated chains are thus formed because intrapore diffusional limitations prevent the hydrogenation of larger species. In addition, the degree of substitution (non-branched, methyl, dimethyl, trimethyl, etc.) of the unsaturated molecules mattered concerning the hydrogenation rate. Linear olefins were more rapidly hydrogenated than branched components. Apart from hydrogenation, skeletal isomerization reactions are also important, especially at high olefin conversions. These fundamental insights demonstrate the capability of GC-VUV to analyze hydrotreated plastic pyrolysis oils and provide essential knowledge for designing and optimizing the hydrotreatment processes of plastic waste pyrolysis oils.

The Operando Nature of Isobutene Adsorbed in Zeolite H-SSZ-13 Unraveled by Machine Learning Potentials Beyond DFT Accuracy

Unraveling the nature of adsorbed olefins in zeolites is crucial to understand numerous zeolite-catalyzed processes. A well-grounded theoretical description critically depends on both an accurate determination of the potential energy surface (PES) and a reliable account of entropic effects at operating conditions. Herein, we propose a transfer learning approach to perform random phase approximation (RPA) quality enhanced sampling molecular dynamics simulations, thereby approaching chemical accuracy on both the determination and exploration of the PES. The proposed methodology is used to investigate isobutene adsorption in H-SSZ-13 as prototypical system to estimate the relative stability of physisorbed olefins, carbenium ions and surface alkoxide species (SAS) in Brønsted-acidic zeolites. We show that the tert-butyl carbenium ion formation is highly endothermic and no entropic stabilization is observed compared to the physisorbed complex within H-SSZ-13. Hence, its predicted concentration and lifetime are negligible, making a direct experimental observation unlikely. Yet, it remains a shallow minimum on the free energy surface over the whole considered temperature range (273-873 K), being therefore a short-lived reaction intermediate rather than a transition state species.

Predicting the effect of framework and hydrocarbon structure on the zeolite-catalyzed beta-scission

Developing improved zeolites is essential in novel sustainable processes such as the catalytic pyrolysis of plastic waste. This study used density functional theory to investigate how alkyl chain length, unsaturated bonds, and branching affect β-scission kinetics in four zeolite frameworks, a key reaction in hydrocarbon cracking. The activation enthalpy was evaluated for a wide variety of 23 hydrocarbons, with 6 to 12 carbon atoms, in FAU, MFI, MOR, and TON. The consideration of both branched and linear olefin and diolefin reactants for the β-scission indicates how the reactant structure influences the intrinsic cracking kinetics, which is especially relevant for the catalytic cracking of plastic waste feedstocks. Intrinsic chemical effects, such as resonance stabilization, the inductive effect, and pore stabilization were found to provide an essential contribution to the activation enthalpy. Additionally, a predictive group additive model incorporating a novel so-called "pore confinement descriptor" was developed for fast prediction of the β-scission activation barrier of a wide range of molecules in the four zeolites. The obtained model can serve as an input for detailed kinetic models in zeolite-catalyzed cracking reactions. The acquired fundamental insights in the cracking of hydrocarbons, relevant for renewable feedstocks, correspond well with experimental observations and will facilitate an improved rational zeolite design.

Confinement-Driven Dimethyl Ether Carbonylation in Mordenite Zeolite as an Ultramicroscopic Reactor

ConspectusThe conversion of C1 molecules to methyl acetate through the carbonylation of dimethyl ether in mordenite zeolite is an appealing reaction and a crucial step in the industrial coal-to-ethanol process. Mordenite zeolite has large 12-membered-ring (12MR) channels (7.0 × 6.5 Å2) and small 8MR channels (5.7 × 2.6 Å2) connected by a side pocket (4.8 × 3.4 Å2), and this unique pore architecture supplies its high catalytic activity to the key step of carbonylation. However, the reaction mechanism of carbonylation in mordenite zeolite is not thoroughly established in that it is able to explain all experimental phenomena and improve its industrial applications, and the classical potential energy surface exerted by static density function theory calculations cannot reflect the reaction kinetics under realistic conditions because the diffusion kinetics of bulk DME (kinetic dimeter: 4.5 Å) and methyl acetate (MA, kinetic dimeter: 5.5 Å) were not well considered and their restricted diffusion in the narrow side pocket and 8MR channels may greatly alter the integrated kinetics of DME carbonylation in mordenite zeolite. Moreover, the precise illustration of the dynamic behaviors of the ketene intermediate and its derivatives (surface acetate and acylium ion) confined within various voids in mordenite has not been effectively portrayed.Advanced ab initio molecular dynamics (AIMD) simulations with or without the acceleration of enhanced sampling methods provide tremendous opportunities for operando modeling of both reaction and diffusion processes and further identify the geometrical structure and chemical properties of the reactants, intermediates, and products in the different confined voids of mordenite under realistic reaction conditions, which enables high consistency between computations and experiments.In this Account, the carbonylation process in mordenite is comprehensively described by the results of decades of continuous research and newly acquired knowledge from both multiscale simulations and in-(ex-)situ spectroscopic experiments. Three primary steps (DME demethylation to surface methoxy species (SMS), carbon-carbon bond coupling between SMS and CO to acetyl species, and methyl acetate formation by acetyl species and methanol/DME) have been respectively studied with a careful consideration of different molecular factors (reactant distribution, concentration, and attack mode). By utilizing the free-energy surface of diffusion and reaction obtained from AIMD simulations, a comprehensive reaction/diffusion kinetic model was formulated for the first time, illustrating the entire zeolite catalytic process. In this context, a comprehensive and informative analysis of the reaction kinetics of carbonylation in mordenite, including the function of the 12MR channels, 8MR channels, and side pockets in the adsorption, diffusion, and reaction of DME carbonylation, was performed. The different channels of mordenite play different roles in all ordered reaction steps, illustrating a highly organized ultramicroscopic reactor that is encompassed.

Selectivity descriptors of the catalytic n-hexane cracking process over 10-membered ring zeolites

Zeolite-mediated catalytic cracking of alkanes is pivotal in the petrochemical and refining industry, breaking down heavier hydrocarbon feedstocks into fuels and chemicals. Its relevance also extends to emerging technologies such as biomass and plastic valorization. Zeolite catalysts, with shape selectivity and selective adsorption capabilities, enhance efficiency and sustainability due to their well-defined network of pores, dimensionality, cages/cavities, and channels. This study focuses on the alkane cracking over 10-membered ring (10-MR) zeolites under industrially relevant conditions. Through a series of characterizations, including operando UV-vis spectroscopy and solid-state NMR spectroscopy, we intend to address mechanistic debates about the alkane cracking mechanism, aiming to understand the dependence of product selectivity on zeolite topologies. The findings highlight topology-dependent mechanisms, particularly the role of intersectional void spaces in zeolite ZSM-5, influencing aromatic-based product selectivity. This work provides a unique understanding of zeolite-catalyzed hydrocarbon conversion, linking alkane activation steps to the traditional hydrocarbon pool mechanism, contributing to the fundamental knowledge of this crucial industrial process.

Reference-Quality Free Energy Barriers in Catalysis from Machine Learning Thermodynamic Perturbation Theory

For the first time, we report calculations of the free energies of activation of cracking and isomerization reactions of alkenes that combine several different electronic structure methods with molecular dynamics simulations. We demonstrate that the use of a high level of theory (here Random Phase Approximation-RPA) is necessary to bridge the gap between experimental and computed values. These transformations, catalyzed by zeolites and proceeding via cationic intermediates and transition states, are building blocks of many chemical transformations for valorization of long chain paraffins originating, e.g., from plastic waste, vegetable oils, Fischer-Tropsch waxes or crude oils. Compared with the free energy barriers computed at the PBE+D2 production level of theory via constrained ab initio molecular dynamics, the barriers computed at the RPA level by the application of Machine Learning thermodynamic Perturbation Theory (MLPT) show a significant decrease for isomerization reaction and an increase of a similar magnitude for cracking, yielding an unprecedented agreement with the results obtained by experiments and kinetic modeling.

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'.

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.

Reaction Pathway of 1-Decene Cracking to Produce Light Olefins over H-ZSM-5 at Ultrahigh Temperature

The effect of reaction temperature and weight hourly space velocity (WHSV) on the reaction of 1-decene cracking to ethylene and propylene over H-ZSM-5 zeolite was investigated. Also, the thermal cracking reaction of 1-decene was studied by cracking over quartz sand as blank. It was observed that 1-decene undergoes a significant thermal cracking reaction above 600 °C over quartz sand. In the range of 500-750 °C, the conversion remained above 99% for 1-decene cracking over H-ZSM-5, and the catalytic cracking dominated even at 750 °C. With the increase in temperature, the yields of ethylene and propylene gradually increased, and the yields of alkanes and aromatics also increased. The low WHSV was favorable for the yield of light olefins. With the increase of the WHSV, the yields of ethylene and propylene decrease. However, at low WHSV, secondary reactions were accelerated, and the yields of alkanes and aromatics increased significantly. In addition, the possible main and side reaction routes of the 1-decene cracking reaction were proposed based on product distribution.

Dual-Site Model for Ab Initio Calculations of Gibbs Free Energies and Enthalpies of Adsorption: Methane in Zeolite Mobile Five (H-MFI)

Quantum chemical hybrid MP2:PBE+D2 calculations in combination with molecular statistics are employed to calculate enthalpies and Gibbs free energies of adsorption for CH₄ at Brønsted acid sites [bridging Si-O(H)-Al groups] and silica wall sites (Si-O-Si) of the proton form of zeolite MFI (H-ZSM-5) and its purely siliceous analogue Silicalite-1. A Langmuir model is adopted to calculate the amounts of CH₄ adsorbed at each type of site from the Gibbs free energies. The combination of these results according to the ratio of silica wall sites and Brønsted acid sites in the sample yields adsorption isotherms for zeolites with different Si/Al ratios. The zero-coverage isosteric heats of adsorption, calculated as thermal averages of the adsorption enthalpies of the individual sites, vary between 20.2 kJ/mol for the pore wall site and 29.2 kJ/mol for the acid site and agree well within ±1 kJ/mol with experimental results.

Toward Computing Accurate Free Energies in Heterogeneous Catalysis: a Case Study for Adsorbed Isobutene in H-ZSM-5

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.

First-Generation Organic Reaction Intermediates in Zeolite Chemistry and Catalysis

Zeolite chemistry and catalysis are expected to play a decisive role in the next decade(s) to build a more decentralized renewable feedstock-dependent sustainable society owing to the increased scrutiny over carbon emissions. Therefore, the lack of fundamental and mechanistic understanding of these processes is a critical "technical bottleneck" that must be eliminated to maximize economic value and minimize waste. We have identified, considering this objective, that the chemistry related to the first-generation reaction intermediates (i.e., carbocations, radicals, carbenes, ketenes, and carbanions) in zeolite chemistry and catalysis is highly underdeveloped or undervalued compared to other catalysis streams (e.g., homogeneous catalysis). This limitation can often be attributed to the technological restrictions to detect such "short-lived and highly reactive" intermediates at the interface (gas-solid/solid-liquid); however, the recent rise of sophisticated spectroscopic/analytical techniques (including under in situ/operando conditions) and modern data analysis methods collectively compete to unravel the impact of these organic intermediates. This comprehensive review summarizes the state-of-the-art first-generation organic reaction intermediates in zeolite chemistry and catalysis and evaluates their existing challenges and future prospects, to contribute significantly to the "circular carbon economy" initiatives.

Catalytic Conversion of Alkenes on Acidic Zeolites: Automated Generation of Reaction Mechanisms and Lumping Technique

Acid-catalyzed hydrocarbon transformations are essential for industrial processes, including oligomerization, cracking, alkylation, and aromatization. However, these chemistries are extremely complex, and computational (automatic) reaction network generation is required to capture these intricacies. The approach relies on the concept that underlying mechanisms for the transformations can be described by a limited number of reaction families applied to various species, with both gaseous and protonated intermediate species tracked. Detailed reaction networks can then be tailored to each industrially relevant process for better understanding or for application in kinetic modeling, which is demonstrated here. However, we show that these networks can grow very large (thousands of species) when they are bound by typical carbon number and rank criteria, and lumping strategies are required to decrease computational expense. For acid-catalyzed hydrocarbon transformations, we propose lumping isomers based on carbon number, branch number, and ion position to reach high carbon limits while maintaining the high resolution of species. Two case studies on propene oligomerization verified the lumping technique in matching a fully detailed model as well as experimental data.

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.

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.

Trimethyloxonium ion – a zeolite confined mobile and efficient methyl carrier at low temperatures: a DFT study coupled with microkinetic analysis

The reaction network of ethene methylation over H-ZSM-5, including methanol dehydration, ethene methylation, and C3H7+ conversion, is investigated by employing a multiscale approach combining DFT calculations and microkinetic modeling.

Ab initio molecular dynamics with enhanced sampling in heterogeneous catalysis

Enhanced sampling ab initio simulations enable to study chemical phenomena in catalytic systems including thermal effects & anharmonicity, & collective dynamics describing enthalpic & entropic contributions, which can significantly impact on reaction free energy landscapes.

Chemically Accurate Vibrational Free Energies of Adsorption from Density Functional Theory Molecular Dynamics: Alkanes in Zeolites

We present a methodology to compute, at reduced computational cost, Gibbs free energies, enthalpies, and entropies of adsorption from molecular dynamics. We calculate vibrational partition functions from vibrational energies, which we obtain from the vibrational density of states by projection on the normal modes. The use of a set of well-chosen reference structures along the trajectories accounts for the anharmonicities of the modes. For the adsorption of methane, ethane, and propane in the H-CHA zeolite, we limit our treatment to a set of vibrational modes localized at the adsorption site (zeolitic OH group) and the alkane molecule interacting with it. Only two short trajectories (1–20 ps) are required to reach convergence (<1 kJ/mol) for the thermodynamic functions. The mean absolute deviations from the experimentally measured values are 2.6, 2.8, and 4.7 kJ/mol for the Gibbs free energy, the enthalpy, and the entropy term (−TΔS), respectively. In particular, the entropy terms show a major improvement compared to the harmonic approximation and almost reach the accuracy of the previous use of anharmonic frequencies obtained with curvilinear distortions of individual modes. The thermodynamic functions so obtained follow the trend of the experimental values for methane, ethane, and propane, and the Gibbs free energy of adsorption at experimental conditions is correctly predicted to change from positive for methane (5.9 kJ/mol) to negative for ethane (−4.8 kJ/mol) and propane (−7.1 kJ/mol).

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.

In Situ Observation of Non-Classical 2-Norbornyl Cation in Confined Zeolites at Ambient Temperature

Carbonium ions are an important class of reaction intermediates, but their dynamic evolution is difficult to be monitored by in situ techniques under experimental conditions because of their extremely short lifetime. Probably the most famous case is 2-norbornyl cation (2NB⁺ ): its existing form (classical or non-classical) had been debated for decades, until the concrete proof of non-classical geometry was achieved by X-ray crystallographic characterization at ultra-low temperature (40 K) and super acidic environment. However, we lack the understanding about 2NB⁺ at ambient conditions. Herein, by taking advantage of the confinement effect and delocalized acidic environment of zeolites, we successfully stabilized 2NB⁺ and unequivocally confirmed its "non-classical" structure inside the ZSM-5 zeolite by ab initio molecular dynamics simulations and ¹³ C solid-state nuclear magnetic resonance experiments. It is the first time to in situ observe the non-classical 2NB⁺ without the super acidic environment at ambient temperature, which provides a new strategy to expand the carbocation chemistry.

Anharmonic Correction to Adsorption Free Energy from DFT-Based MD Using Thermodynamic Integration

Adsorption processes are often governed by weak interactions for which the estimation of entropy contributions by means of the harmonic approximation is prone to be inaccurate. Thermodynamic integration (TI) from the harmonic to the fully interacting system (λ-path integration) can be used to compute anharmonic corrections. Here, we combine TI with (curvilinear) internal coordinates in periodic systems to make the formalism available in computational studies. Our implementation of ab initio molecular dynamics in VASP is independent of the reaction path and can be thus applied to study adsorption processes relative to the gas phase and does hence provide a useful tool for computational catalysis. We discuss the application of the approach on three model systems for which exact semianalytical solutions exist and illustrate and quantify the importance of anharmonic vibrations, hindered rotations, and hindered translations (dissociation). Eventually, we apply the method to study the adsorption of small adsorbates in a zeolite (H-SSZ-13).

Dynamic Features of Transition States for β-Scission Reactions of Alkenes over Acid Zeolites Revealed by AIMD Simulations

Zeolite-catalyzed alkene cracking is key to optimize the size of hydrocarbons. The nature and stability of intermediates and transition states (TS) are, however, still debated. We combine transition path sampling and blue moon ensemble density functional theory simulations to unravel the behavior of C₇ alkenes in CHA zeolite. Free energy profiles are determined, linking π-complexes, alkoxides and carbenium ions, for B₁ (secondary to tertiary) and B₂ (tertiary to secondary) β-scissions. B₁ is found to be easier than B₂ . The TS for B₁ occurs at the breaking of the C-C bond, while for B₂ it is the proton transfer from propenium to the zeolite. We highlight the dynamic behaviors of the various intermediates along both pathways, which reduce activation energies with respect to those previously evaluated by static approaches. We finally revisit the ranking of isomerization and cracking rate constants, which are crucial for future kinetic studies.

Chlorella to fuel conversion on amphiphilic SO3H-SBA-15 catalysts: Pyrolysis characteristics and kinetics

The aim of this work was to propose a novel process to make Chlorella pyrolyzed and in situ upgraded to fuel over amphiphilic SO₃H-SBA-15 catalysts. This strategy is developed to build a Pickering emulsion system through the w/o (water/decalin) droplets. Chlorella catalytic pyrolysis has been conducted under the different heating rates to get the activation energy 166 kJ/mol (α = 0.5) according to the kinetic-free model. Palmitic acid, as a model compound, was employed for TG and DRIFTS analysis to elucidate the pyrolysis and deoxygenation reaction pathway. n-hexadecane pyrolysis at 3 MPa N₂ illustrated the peak cracking temperature declining from thermally 422 °C to catalytically 413 °C. N₂ physisorption of the fresh and post-reaction catalysts indicated that there is little catalyst decay. With improved thermal stability and hydrophobicity, the SO₃H-SBA-15 catalysts showed enhanced performance for Chlorella pyrolysis, and revealed the promising application for better fuel production in aqueous conversion.

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.