Consideration of the Aluminum Distribution in Zeolites in Theoretical and Experimental Catalysis Research

Enhancing the Acidity Window of Zeolites by Low-Temperature Template Oxidation with Ozone

Revisiting the impact of the first and often deemed trivial postsynthetic step, i.e., a high-temperature oxidative calcination to remove organic templates, increases our understanding of thermal acid site evolution and Al distributions. An unprecedented degree of control over the acidity of high-silica zeolites (SSZ-13) was achieved by using a low-temperature ozonation approach. Fourier transform infrared spectroscopy of adsorbed probe molecules and solid-state NMR spectroscopy reveal the complexity of the thermal evolution of acid sites. Low-temperature activated (ozonated) zeolites maintain the original Brønsted acidity content and high defect content and have virtually no Lewis acidity. They also preserve the "as-made" Al distribution after crystallization and show a clear link between synthesis conditions and divalent cation capacity, as measured with aqueous cobalt ion uptake. The synthesis protocol is found to be the main contributor to Al proximity, yielding record high exchange capacity when ozonated. After conventional calcination at 500-600 °C, however, the presence of water leads to the gradual depletion of Brønsted acid sites, in particular, in small crystals. This work indicates that low-temperature ozonation followed by thermal activation at different temperatures can be used as a novel tool for tuning the amount and nature of acid sites, providing insights into the activity of zeolites in acid-catalyzed reactions, such as CO2 hydrogenation to dimethyl ether, and thereby expanding the possibilities of rational acidity tuning.

Advances and challenges in designing active site environments in zeolites for Brønsted acid catalysis

Zeolites contain proton active sites in diverse void environments that stabilize the reactive intermediates and transition states formed in converting hydrocarbons and oxygenates to chemicals and energy carriers. The catalytic diversity that exists among active sites in voids of varying sizes and shapes, even within a given zeolite topology, has motivated research efforts to position and quantify active sites within distinct voids (synthesis-structure) and to link active site environment to catalytic behavior (structure-reactivity). This Feature Article describes advances and challenges in controlling the position of framework Al centers and associated protons within distinct voids during zeolite synthesis or post-synthetic modification, in identifying and quantifying distinct active site environments using characterization techniques, and in determining the influence of active site environments on catalysis. During zeolite synthesis, organic structure directing agents (SDAs) influence Al substitution at distinct lattice positions via intermolecular interactions (e.g., electrostatics, hydrogen bonding) that depend on the size, structure, and charge distribution of organic SDAs and their mobility when confined within zeolitic voids. Complementary post-synthetic strategies to alter intrapore active site distributions include the selective removal of protons by differently-sized titrants or unreactive organic residues and the selective exchange of framework heteroatoms of different reactivities, but remain limited to certain zeolite frameworks. The ability to identify and quantify active sites within distinct intrapore environments depends on the resolution with which a given characterization technique can distinguish Al T-site positions or proton environments in a given zeolite framework. For proton sites in external unconfined environments, various (post-)synthetic strategies exist to control their amounts, with quantitative methods to distinguish them from internal sites that largely depend on using stoichiometric or catalytic probes that only interact with external sites. Protons in different environments influence reactivity by preferentially stabilizing larger transition states over smaller precursor states and influence selectivity by preferentially stabilizing or destabilizing competing transition states of varying sizes that share a common precursor state. We highlight opportunities to address challenges encountered in the design of active site environments in zeolites by closely integrating precise (post-)synthetic methods, validated characterization techniques, well-defined kinetic probes, and properly calibrated theoretical models. Further advances in understanding the molecular details that underlie synthesis-structure-reactivity relationships for active site environments in zeolite catalysis can accelerate the predictive design of tailored zeolites for desired catalytic transformations.

Chemically accurate predictions for water adsorption on Brønsted sites of zeolite H-MFI

We investigate the adsorption of water molecules in the zeolite H-MFI at isolated Brønsted acid sites (BAS) for loadings of 1, 2, and 3 H2O/BAS. We consider two approaches to the O3Al-O(H)-Si sites: the Brønsted-type approach of H2O to the acidic proton and the Lewis-type approach to the aluminium atom of the AlO4 tetrahedron. From the twelve crystallographically inequivalent framework sites for Al, a representative set of six active site positions is chosen. For them, we calculate CCSD(T)-quality adsorption energies at MP2-quality adsorption structures for different approaches, 48 in total. The Brønsted-type approach is favoured for most cases but the Lewis-type approach has similar stability for some framework positions. We predict heats of adsorption per molecule ranging from 60 to 76, 56 to 65, and 56 to 64 kJ mol-1 for loadings of 1, 2, and 3 H2O/BAS, respectively. For 1 H2O/BAS, the experimental result (70 kJ mol-1) falls into the range of our predictions, whereas for 2 and 3 H2O/BAS, the measured adsorption heats per molecule (74 and 70 kJ mol-1, respectively) are larger than our predictions. For 2 H2O/BAS, the ion-pair structure generated by proton transfer to the water dimer competes with the neutral adsorption complex. The DFT adsorption energies (PBE+D2) deviate significantly from the CCSD(T)-quality reference energies, by up to 25 kJ mol-1 for 1 H2O/BAS, 25 kJ mol-1 per H2O for 2 H2O/BAS, and 18 kJ mol-1 per H2O for 3 H2O/BAS. Specifically, PBE+D2 overstabilises the ion-pair structure, i.e. in many cases the PBE+D2 error is much larger for ionic than for neutral adsorption structures.

Synthetic Placement of Active Sites in MFI Zeolites for Selective Toluene Methylation to para-Xylene

Brønsted acidic zeolites are ubiquitous catalysts in fuel and chemical production. Broadening the catalytic diversity of a given zeolite requires strategies to manipulate the acid site placement at framework positions within distinct microporous locations. Here, we combine experiment and theory to elucidate how intermolecular interactions between organic structure-directing agents (OSDAs) and framework Al centers influence the placement of H+ sites in distinct void environments of MFI zeolites and demonstrate the catalytic consequences of active site location on kinetically controlled (403 K) toluene methylation to xylene regioisomers. Kinetic measurements, interpreted using mechanism-derived rate expressions and transition state theory, alongside density functional theory (DFT) calculations show that larger intersection environments similarly stabilize all three xylene isomer transition states without altering well-established aromatic substitution patterns (ortho/para/meta ∼ 60%:30%:10%), while smaller channel environments preferentially destabilize transition states that form bulkier ortho- and meta-isomers, thereby resulting in high intrinsic para-xylene selectivity (∼80%). DFT calculations reveal that the flexibility of nonconventional OSDAs (e.g., 1,4-diazabicyclo[2.2.2]octane) to reorient within MFI intersections and their ability to hydrogen-bond to form protonated complexes favor the placement of Al in smaller channel environments compared to conventional quaternary OSDAs (e.g., tetra-n-propylammonium). These molecular-level insights establish a mechanistic link between OSDA structure, active site placement, and transition state stability in MFI zeolites and provide active site design strategies that are orthogonal to crystallite design approaches harnessing complex reaction-diffusion phenomena to enhance regioisomer selectivity in the industrial production of valuable polymer precursors.

Understanding Alkene Interaction with Metal-Modified Zeolites: Thermodynamics and Mechanism of Bonding in the π-Complex

Zeolites modified with metal cations are perspective catalysts for converting light alkenes to valuable chemicals. A crucial step of the transformation is an alkene interaction with zeolite to afford π-complex with metal cations. The mechanism of alkene bonding with cations is still unclear. To address this problem, propene adsorption on H+ (Bro̷nsted acid site), Na+, Ca2+, Zn2+, Co2+, Cu2+, Cu+, and Ag+ cationic sites in ZSM-5 zeolite has been studied by quantum chemical calculations in terms of adsorption enthalpy, νC═C frequency, and natural bond orbital (NBO) analysis together with natural energy decomposition analysis (NEDA). It is revealed that the conventional concept of σ- and π-bonding is only partially applicable to alkene interaction with metal cations in zeolites. The orbital interaction between an alkene molecule and a metal site is more complex. Several different bonding mechanisms have been identified depending on the nature and electron configuration of the metal cation. This finding explains the complex correlations observed for propene π-complex stability and νC═C frequency shift or charge transfer from the alkene molecule. The results provide the basis for further understanding the interactions between alkenes and inorganic solid Bro̷nsted and Lewis acids.

The need for operando modelling of 27Al NMR in zeolites: the effect of temperature, topology and water

Solid state (ss-) 27Al NMR is one of the most valuable tools for the experimental characterization of zeolites, owing to its high sensitivity and the detailed structural information which can be extracted from the spectra. Unfortunately, the interpretation of ss-NMR is complex and the determination of aluminum distributions remains generally unfeasible. As a result, computational modelling of 27Al ss-NMR spectra has grown increasingly popular as a means to support experimental characterization. However, a number of simplifying assumptions are commonly made in NMR modelling, several of which are not fully justified. In this work, we systematically evaluate the effects of various common models on the prediction of 27Al NMR chemical shifts in zeolites CHA and MOR. We demonstrate the necessity of operando modelling; in particular, taking into account the effects of water loading, temperature and the character of the charge-compensating cation. We observe that conclusions drawn from simple, high symmetry model systems such as CHA do not transfer well to more complex zeolites and can lead to qualitatively wrong interpretations of peak positions, Al assignment and even the number of signals. We use machine learning regression to develop a simple yet robust relationship between chemical shift and local structural parameters in Al-zeolites. This work highlights the need for sophisticated models and high-quality sampling in the field of NMR modelling and provides correlations which allow for the accurate prediction of chemical shifts from dynamical simulations.

Analysis of Al site-directing ability of organic structure-directing agents in FER and CHA zeolites: a computational exploration of energetic preferences

This study explores the control of Al location in zeolites by organic structure-directing agents (OSDAs) using atomistic simulations. We examine several zeolite-OSDA complexes to quantify the Al site-directing ability. The results show that OSDAs induce different energetic preferences to direct Al at certain locations. In particular, these effects can be enhanced by OSDAs with N-H moieties. Our findings will be useful for the development of novel OSDAs that can modulate Al site-directing properties.

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.

Impact of Mineralizing Agents on Aluminum Distribution and Acidity of ZSM-5 Zeolites

Intensive research on improving the catalytic properties of zeolites is focused on modulating their acidity and the distribution of associated Al sites. Herein, by studying a series of ZSM-5 zeolites over a broad range of Al content, we demonstrate how the nature of the mineralizing agent (F^- or OH^- ) used in hydrothermal syntheses directly impacts Al sites distribution. The proportions of Al sites, probed by ²⁷ Al NMR, depend on the Si/Al ratio for F^- , but remain identical for OH^- (from Si/Al=30 to 760). This leads to contrasting variations in weak and strong acidities. Such opposite effect of mineralizers is explained by the spatial location of negative charges and the resulting balance between short- and long-range electrostatic interactions. This understanding paves the way for additional and simple opportunities to control zeolites' acidity.

Heteroatom Manipulation of Zeolite Crystallization: Stabilizing Zn-FAU against Interzeolite Transformation

The preparation of metastable zeolites is often restricted to a limited range of synthesis conditions, which is exemplified in commercial syntheses lacking organics to stabilize the crystal structure. In the absence of an organic structure-directing agent, interzeolite transformation is a common phenomenon that can lead to undesirable products or impurities. Many studies have investigated the substitution of Si and Al in zeolite frameworks with alternative elements (heteroatoms) as a means of tailoring the properties of zeolites; however, relatively few studies have systematically explored the impact of heteroatoms on interzeolite transformations and their concomitant effects on zeolite crystallization. In this study, we examine methods to prepare isostructures of faujasite (FAU), which is one of the most commercially relevant zeolites and also a thermodynamically metastable structure. A survey of multivalent elements revealed that zinc is capable of stabilizing FAU at high temperatures and inhibiting its frequent transformation to zeolite gismondine (GIS). Using combined experimental and computational studies, we show that zinc alters the chemical nature of growth mixtures by sequestering silicates. Zinc heteroatoms incorporate in the FAU framework with a loading-dependent coordination. Our collective findings provide an improved understanding of driving forces for the FAU-to-GIS interzeolite transformation where we observe that heteroatoms (e.g., zinc) can stabilize zeolite FAU over a broad range of synthesis conditions. Given the growing interest in heteroatom-substituted zeolites, this approach to preparing zinc-containing FAU may prove applicable to a broader range of zeolite structures.

17O Labeling Reveals Paired Active Sites in Zeolite Catalysts

Current needs for extending zeolite catalysts beyond traditional gas-phase hydrocarbon chemistry demand detailed characterization of active site structures, distributions, and hydrothermal impacts. A broad suite of homonuclear and heteronuclear NMR correlation experiments on dehydrated H-ZSM-5 catalysts with isotopically enriched ¹⁷O frameworks reveals that at least two types of paired active sites exist, the amount of which depends on the population of fully framework-coordinated tetrahedral Al (Al(IV)-1) and partially framework-coordinated tetrahedral Al (Al(IV)-2) sites, both of which can be denoted as (SiO)_4-n-Al(OH)_n. The relative amounts of Al(IV)-1 and Al(IV)-2 sites, and subsequent pairing, cannot be inferred from the catalyst Si/Al ratio, but depend on synthetic and postsynthetic modifications. Correlation experiments demonstrate that, on average, acidic hydroxyl groups from Al(IV)-1/Al(IV)-2 pairs are closer to one another than those from Al(IV)-1/Al(IV)-1 pairs, as supported by computational DFT calculations. Through-bond and through-space polarization transfer experiments exploiting ¹⁷O nuclei reveal a number of different acidic hydroxyl groups in varying Si/Al catalysts, the relative amounts of which change following postsynthetic modifications. Using room-temperature isotopic exchange methods, it was determined that ¹⁷O was homogeneously incorporated into the zeolite framework, while ¹⁷O → ²⁷Al polarization transfer experiments demonstrated that ¹⁷O incorporation does not occur for extra-framework Al_nO_m species. Data from samples exposed to controlled hydrolysis indicates that nearest neighbor Al pairs in the framework are more susceptible to hydrolytic attack. The data reported here suggest that Al(IV)-1/Al(IV)-2 paired sites are synergistic sites leading to increased reactivity in both low- and high-temperature reactions. No evidence was found for paired framework/nonframework sites.

Driving the Active Site Incorporation in Zeolitic Materials via the Organic Structure-Directing Agent Through Development of H-Bonds with Hydroxyl Groups

(1S,2S)-N-methyl-pseudoephedrine (MPS) was used as organic structure-directing agent (OSDA) for the synthesis of Mg-doped nanoporous aluminophosphates. This molecule displays a particular conformational behavior, where the presence of H-bond donor and acceptor groups provide a rigid conformational space with one asymmetric conformation preferentially occurring. MPS drives the crystallization of Mg-containing AFI materials. Characterization of these materials shows that the OSDA incorporate as protonated species, arranged as head-to-tail monomers. Combination of three-dimensional electron diffraction with high-resolution synchrotron powder X-ray diffraction allowed to locate both the Mg and the organic species. Interestingly, results showed that the spatial incorporation of Mg is driven by the hydroxyl groups of the organic cation through the development of H-bonds with negatively-charged MgO4 tetrahedra. This work demonstrates that H-bond forming groups can be used to drive the spatial incorporation of low-valent dopants within zeolitic frameworks, a highly desired aim in order to control their catalytic activity and selectivity.

Interconversion of Atomically Dispersed Platinum Cations and Platinum Clusters in Zeolite ZSM-5 and Formation of Platinum gem-Dicarbonyls

Catalysts composed of platinum dispersed on zeolite supports are widely applied in industry, and coking and sintering of platinum during operation under reactive conditions require their oxidative regeneration, with the platinum cycling between clusters and cations. The intermediate platinum species have remained only incompletely understood. Here, we report an experimental and theoretical investigation of the structure, bonding, and local environment of cationic platinum species in zeolite ZSM-5, which are key intermediates in this cycling. Upon exposure of platinum clusters to O₂ at 700 °C, oxidative fragmentation occurs, and Pt²⁺ ions are stabilized at six-membered rings in the zeolite that contain paired aluminum sites. When exposed to CO under mild conditions, these Pt²⁺ ions form highly uniform platinum gem-dicarbonyls, which can be converted in H₂ to Pt^δ+ monocarbonyls. This conversion, which weakens the platinum-zeolite bonding, is a first step toward platinum migration and aggregation into clusters. X-ray absorption and infrared spectra provide evidence of the reductive and oxidative transformations in various gas environments. The chemistry is general, as shown by the observation of platinum gem-dicarbonyls in several commercially used zeolites (ZSM-5, Beta, mordenite, and Y).

Increasing the Number of Aluminum Atoms in T3 Sites of a Mordenite Zeolite by Low-Pressure SiCl4 Treatment to Catalyze Dimethyl Ether Carbonylation

Controlling the location of aluminum atoms in a zeolite framework is critical for understanding structure-performance relationships of catalytic reaction systems and tailoring catalyst design. Herein, we report a strategy to preferentially relocate mordenite (MOR) framework Al atoms into the desired T₃ sites by low-pressure SiCl₄ treatment (LPST). High-field ²⁷ Al NMR was used to identify the exact location of framework Al for the MOR samples. The results indicate that 73 % of the framework Al atoms were at the T₃ sites after LPST under optimal conditions, which leads to controllably generating and intensifying active sites in MOR zeolite for the dimethyl ether (DME) carbonylation reaction with higher methyl acetate (MA) selectivity and much longer lifetime (25 times). Further research reveals that the Al relocation mechanism involves simultaneous extraction, migration, and reinsertion of Al atoms from and into the parent MOR framework. This unique method is potentially applicable to other zeolites to control Al location.

Designing Sites in Heterogeneous Catalysis: Are We Reaching Selectivities Competitive With Those of Homogeneous Catalysts?

A critical review of different prominent nanotechnologies adapted to catalysis is provided, with focus on how they contribute to the improvement of selectivity in heterogeneous catalysis. Ways to modify catalytic sites range from the use of the reversible or irreversible adsorption of molecular modifiers to the immobilization or tethering of homogeneous catalysts and the development of well-defined catalytic sites on solid surfaces. The latter covers methods for the dispersion of single-atom sites within solid supports as well as the use of complex nanostructures, and it includes the post-modification of materials via processes such as silylation and atomic layer deposition. All these methodologies exhibit both advantages and limitations, but all offer new avenues for the design of catalysts for specific applications. Because of the high cost of most nanotechnologies and the fact that the resulting materials may exhibit limited thermal or chemical stability, they may be best aimed at improving the selective synthesis of high value-added chemicals, to be incorporated in organic synthesis schemes, but other applications are being explored as well to address problems in energy production, for instance, and to design greener chemical processes. The details of each of these approaches are discussed, and representative examples are provided. We conclude with some general remarks on the future of this field.

Long-Range Spatial Distribution of Single Aluminum Sites in Zeolites

How aluminum distributes during synthesis and rearranges after processing within the zeolite framework is a central question in heterogeneous catalysis, as it determines the structure and location of the catalytically active sites of the one of the most important classes of industrial catalysts. Here, exploiting the dipolar interaction between paramagnetic metal ions, we derive the spatial distribution of single aluminum sites within the ZSM-5 zeolite framework in the nanometer range, in polycrystalline samples lacking long-range order. We use a Monte Carlo approach to validate the findings on a pristine ZSM-5 sample and demonstrate that the method is sensitive enough to monitor aluminum redistribution induced in the framework by chemical stress.

Second-Sphere Lattice Effects in Copper and Iron Zeolite Catalysis

Transition-metal-exchanged zeolites perform remarkable chemical reactions from low-temperature methane to methanol oxidation to selective reduction of NOx pollutants. As with metalloenzymes, metallozeolites have impressive reactivities that are controlled in part by interactions outside the immediate coordination sphere. These second-sphere effects include activating a metal site through enforcing an "entatic" state, controlling binding and access to the metal site with pockets and channels, and directing radical rebound vs cage escape. This review explores these effects with emphasis placed on but not limited to the selective oxidation of methane to methanol with a focus on copper and iron active sites, although other transition-metal-ion zeolite reactions are also explored. While the actual active-site geometric and electronic structures are different in the copper and iron metallozeolites compared to the metalloenzymes, their second-sphere interactions with the lattice or the protein environments are found to have strong parallels that contribute to their high activity and selectivity.

Effect of crystallization mechanism on zeolite BEA textural and acidic properties

Two zeolite BEA catalysts have been synthesized via alternative mechanistic pathways involving solution-mediated transportation and solid-solid hydrogel rearrangement. In the synthetic procedures selected, the composition of starting reaction mixtures and...

Anchoring Boron Atom to the Specific Tetrahedral Sites of Borosilicate MFI by Imidazolium-based Molecules

Regulating heteroatom substituted tetrahedral sites (T sites) within a zeolite crystal significantly influences the zeolite properties. Herein, the as-synthesized MFI zeolite with tailored boron location at various T sites is...

Construction of cucurbit[n]uril-based supramolecular frameworks via host-guest inclusion and functional properties thereof

Frameworks utilizing cucurbit[n]uril-based chemistry build on the rapid developments in the fields of metal-organic frameworks (MOFs), covalent-organic frameworks (COFs), and supramolecular organic frameworks (SOFs), and as porous materials have found...

Designing a Multifunctional Catalyst for the Direct Production of Gasoline-Range Isoparaffins from CO2

The production of carbon-neutral fuels from CO₂ presents an avenue for causing an appreciable effect in terms of volume toward the mitigation of global carbon emissions. To that end, the production of isoparaffin-rich fuels is highly desirable. Here, we demonstrate the potential of a multifunctional catalyst combination, consisting of a methanol producer (InCo) and a Zn-modified zeolite beta, which produces a mostly isoparaffinic hydrocarbon mixture from CO₂ (up to ∼85% isoparaffin selectivity among hydrocarbons) at a CO₂ conversion of >15%. The catalyst combination was thoroughly characterized via an extensive complement of techniques. Specifically, operando X-ray absorption spectroscopy (XAS) reveals that Zn (which plays a crucial role of providing a hydrogenating function, improving the stability of the overall catalyst combination and isomerization performance) is likely present in the form of Zn₆O₆ clusters within the zeolite component, in contrast to previously reported estimations.

Selectivity descriptors for the direct hydrogenation of CO2 to hydrocarbons during zeolite-mediated bifunctional catalysis

Cascade processes are gaining momentum in heterogeneous catalysis. The combination of several catalytic solids within one reactor has shown great promise for the one-step valorization of C1-feedstocks. The combination of metal-based catalysts and zeolites in the gas phase hydrogenation of CO2 leads to a large degree of product selectivity control, defined mainly by zeolites. However, a great deal of mechanistic understanding remains unclear: metal-based catalysts usually lead to complex product compositions that may result in unexpected zeolite reactivity. Here we present an in-depth multivariate analysis of the chemistry involved in eight different zeolite topologies when combined with a highly active Fe-based catalyst in the hydrogenation of CO2 to olefins, aromatics, and paraffins. Solid-state NMR spectroscopy and computational analysis demonstrate that the hybrid nature of the active zeolite catalyst and its preferred CO2-derived reaction intermediates (CO/ester/ketone/hydrocarbons, i.e., inorganic-organic supramolecular reactive centers), along with 10 MR-zeolite topology, act as descriptors governing the ultimate product selectivity.

The Role of Organic and Inorganic Structure-Directing Agents in Selective Al Substitution of Zeolite

Organic and inorganic structure-directing agents (SDAs) impact Al distributions in zeolite, but the insights into how SDAs manipulate Al distribution have not been elucidated yet. Herein, the roles of different SDAs such as cyclohexylamine (CHA), hexamethylenimine (HMI), and Na⁺ in selective Al substitution of MCM-49 zeolite are investigated comprehensively by multinuclear solid-state NMR. The results demonstrate that MCM-49 synthesized with HMI shows relatively more T₆ and T₇ Al, while more T₂ Al is observed using CHA. The formation of T₂ Al in both MCM-49(HMI) and MCM-49(CHA) is derived from Na⁺, while protonated HMIs show bias in incorporation of T₆ and T₇ Al. Most HMIs are occluded in protonated status, and about half of CHAs are occluded in nonprotonated status. The close spatial proximity between nonprotonated CHAs and Na⁺ synergistically promotes the formation of zeolite structure, leading to more Na⁺ ions occluded in the zeolite channel with preferential T₂ Al substitution.

QM/MM study of the reactivity of zeolite bound methoxy and carbene groups

The conversion of methanol-to-hydrocarbons (MTH) is known to occur via an autocatalytic process in zeolites, where framework-bound methoxy species play a pivotal role, especially during catalyst induction. Recent NMR and FT-IR experimental studies suggest that methoxylated zeolites are able to produce hydrocarbons by a mechanism involving carbene migration and association. In order to understand these observations, we have performed QM/MM computational investigations on a range of reaction mechanisms for the reaction of zeolite bound methoxy and carbene groups, which are proposed to initiate hydrocarbon formation in the MTH process. Our simulations demonstrate that it is kinetically unfavourable for methyl species to form on the framework away from the zeolite acid site, and both kinetically and thermodynamically unfavourable for methyl groups to migrate through the framework and aggregate around an acid site. Formation of carbene moieties was considered as an alternative pathway to the formation of C-C bonds; however, the reaction energy for conversion of a methyl to a carbene is unfavourable. Metadynamics simulations help confirm further that methyl species at the framework acid sites would be more reactive towards formed C₂₊ species, rather than inter-framework migration, and that the role of carbenes in the formation of the first C-C bond will be via a concerted type of mechanism rather than stepwise.

Probing Monomer and Dimer Adsorption Trends in the MFI Framework

Porous aluminosilicates such as zeolites are ubiquitous catalysts for the production of high-value and industrially relevant commodity chemicals, including the conversion of hydrocarbons, amines, alcohols, and others. Bimolecular reactions are an important subclass of reactions that can occur on Brønsted acid sites of a zeolite catalyst. Kinetic modeling of these systems at the process scale requires the interaction energetics of reactants and the active sites to be described accurately. It is generally known that adsorption is a coverage-dependent phenomenon, with lower heats of adsorption observed for molecules at higher coverage. However, few studies have systematically investigated the coadsorption of molecules on a single active site, specifically focusing on the strength of interaction of the second adsorbate after the initial adsorption step. In this work, we quantify the unimolecular and bimolecular adsorption energies of varying adsorbates, including paraffins, olefins, alcohols, amines, and noncondensible gases in the acidic and siliceous ZSM-5 frameworks. As a special case, olefin adsorption was examined for physisorption and chemisorption regimes, characterized by π-complex, framework alkoxide and carbenium ion adsorption, respectively. The effects of functional groups and molecular size were quantified, and correlations that relate the adsorption of the second adsorbate identity to that of the first adsorbate are provided.

Role of the ionic environment in enhancing the activity of reacting molecules in zeolite pores

Tailoring the molecular environment around catalytically active sites allows for the enhancement of catalytic reactivity through a hitherto unexplored pathway. In zeolites, the presence of water creates an ionic environment via the formation of hydrated hydronium ions and the negatively charged framework aluminum tetrahedra. The high density of cation-anion pairs determined by the aluminum concentration of a zeolite induces a high local ionic strength that increases the excess chemical potential of sorbed and uncharged organic reactants. Charged transition states (carbocations for example) are stabilized, which reduces the energy barrier and leads to higher reaction rates. Using the intramolecular dehydration of cyclohexanol on H-MFI zeolites in water, we quantitatively show an enhancement of the reaction rate by the presence of high ionic strength as well as show potential limitations of this strategy.

Discovering Relationships between OSDAs and Zeolites through Data Mining and Generative Neural Networks

Organic structure directing agents (OSDAs) play a crucial role in the synthesis of micro- and mesoporous materials especially in the case of zeolites. Despite the wide use of OSDAs, their interaction with zeolite frameworks is poorly understood, with researchers relying on synthesis heuristics or computationally expensive techniques to predict whether an organic molecule can act as an OSDA for a certain zeolite. In this paper, we undertake a data-driven approach to unearth generalized OSDA-zeolite relationships using a comprehensive database comprising of 5,663 synthesis routes for porous materials. To generate this comprehensive database, we use natural language processing and text mining techniques to extract OSDAs, zeolite phases, and gel chemistry from the scientific literature published between 1966 and 2020. Through structural featurization of the OSDAs using weighted holistic invariant molecular (WHIM) descriptors, we relate OSDAs described in the literature to different types of cage-based, small-pore zeolites. Lastly, we adapt a generative neural network capable of suggesting new molecules as potential OSDAs for a given zeolite structure and gel chemistry. We apply this model to CHA and SFW zeolites generating several alternative OSDA candidates to those currently used in practice. These molecules are further vetted with molecular mechanics simulations to show the model generates physically meaningful predictions. Our model can automatically explore the OSDA space, reducing the amount of simulation or experimentation needed to find new OSDA candidates.

Benchmarking binding energy calculations for organic structure-directing agents in pure-silica zeolites

Molecular modeling plays an important role in the discovery of organic structure-directing agents (OSDAs) for zeolites. By quantifying the intensity of host-guest interactions, it is possible to select cost-effective molecules that maximize binding toward a given zeolite framework. Over the last few decades, a variety of methods and levels of theory have been used to calculate these binding energies. Nevertheless, there is no consensus on the best calculation strategy for high-throughput virtual screening undertakings. In this work, we compare binding affinities from density functional theory (DFT) and Dreiding force field calculations for 272 zeolite-OSDA pairs obtained from static and time-averaged simulations. Enabled by automation software, we show that Dreiding binding energies from the frozen pose method correlate best with DFT energies. They are also less sensitive to the choice of initial lattice parameters and optimization algorithms, as well as less computationally expensive than their time-averaged counterparts. Furthermore, we demonstrate that a broader exploration of the conformation space from molecular dynamics simulations does not provide significant improvements in binding energy trends over the frozen pose method despite being orders of magnitude more expensive. The code and benchmark data are open-sourced and provide robust and computationally efficient guidelines to calculating binding energies in zeolite-OSDA pairs.

Distribution of Aluminum Species in Zeolite Catalysts: 27Al NMR of Framework, Partially-Coordinated Framework, and Non-Framework Moieties

The structure of aluminum-containing moieties in and within zeolite H-ZSM-5 catalysts is a complex function of the elemental composition of the catalyst, synthesis conditions, exposure to moisture, and thermal history. 27Al NMR data collected at field strengths ranging from 7.05 to 35.2 T, i.e., 1H Larmor frequencies from 300 to 1500 MHz, reveal that Al primarily exists as framework or partially coordinated framework species in commercially available dehydrated H-ZSM-5 catalysts with Si/Al ranging from 11.5 to 40. Quantitative direct-excitation and sensitivity-enhanced 27Al NMR techniques applied over the wide range of magnetic field strengths used in this study show that prior to significant hydrothermal exposure, detectable amounts of nonframework Al species do not exist. Two-dimensional 27Al multiple-quantum magic-angle spinning (MQMAS) along with 1H-27Al and 29Si-27Al dipolar correlation (D-HMQC) NMR experiments confirm this conclusion and show that generation of nonframework species following varying severities of hydrothermal exposure are clearly resolved from partially coordinated framework sites. The impact of hydration on the appearance and interpretation of conventional direct-excitation 27Al spectra, commonly used to assess framework and nonframework Al, is discussed. Aluminum sites in dehydrated catalysts, which are representative of typical operating conditions, are characterized by large quadrupole interactions and are best assigned by obtaining data at multiple field strengths. On the basis of the results here, an accurate initial assessment of Al sites in high-Al content MFI catalysts prior to any hydrothermal treatment can be used to guide reaction conditions, anticipate potential water impacts, and identify contributions from hydroxyl groups other than those associated with the framework bridging acid site.

Dynamic Interconversion of Metal Active Site Ensembles in Zeolite Catalysis

Catalysis science is founded on understanding the structure, number, and reactivity of active sites. Kinetic models that consider active sites to be static and noninteracting entities are routinely successful in describing the behavior of heterogeneous catalysts. Yet, active site ensembles often restructure in response to their external environment and even during steady-state catalytic turnover, sometimes requiring non-mean-field kinetic treatments to describe distance-dependent interactions among sites. Such behavior is being recognized more frequently in modern catalysis research, with the advent of experimental methods to quantify turnover rates with increasing precision, an expanding arsenal of operando characterization tools, and computational descriptions of atomic structure and motion at chemical potentials and timescales increasingly relevant to reaction conditions. This review focuses on dynamic changes to metal active site ensembles on zeolite supports, which are silica-based crystalline materials substituted with Al that generate binding sites for isolated and low-nuclearity metal site ensembles. Metal sites can become solvated and mobilized during reaction, facilitating interactions among sites that change their nuclearity and function. Such intersite communication can be regulated by the zeolite support, resulting in non-single-site and potentially non-mean-field kinetic behavior arising from mechanisms of catalytic action that combine elements of those canonically associated with homogeneous and heterogeneous catalysis.We discuss recent literature examples that document dynamic active site behavior in metal-zeolites and outline methodologies to identify and interpret such behavior. We conclude with our outlook on future research directions to develop this evolving branch of catalysis science and harness it for practical applications.

To Every Rule There is an Exception: A Rational Extension of Loewenstein's Rule

Loewenstein's rule, which states that Al-O-Al motifs are energetically unstable, is fundamental to the understanding and design of zeolites. Here, using a combination of electronic structure calculations and lattice models, we show under which circumstances this rule becomes invalid and how it can be rationally extended using the chabasite framework for demonstration.