ZSM-5 Zeolite: Complete Al Bond Connectivity and Implications on Structure Formation from Solid-State NMR and Quantum Chemistry Calculations

Advanced solid-state NMR spectroscopy and its applications in zeolite chemistry

Solid-state NMR spectroscopy (ssNMR) can provide details about the structure, host-guest/guest-guest interactions and dynamic behavior of materials at atomic length scales. A crucial use of ssNMR is for the characterization of zeolite catalysts that are extensively employed in industrial catalytic processes. This review aims to spotlight the recent advancements in ssNMR spectroscopy and its application to zeolite chemistry. We first review the current ssNMR methods and techniques that are relevant to characterize zeolite catalysts, including advanced multinuclear and multidimensional experiments, in situ NMR techniques and hyperpolarization methods. Of these, the methodology development on half-integer quadrupolar nuclei is emphasized, which represent about two-thirds of stable NMR-active nuclei and are widely present in catalytic materials. Subsequently, we introduce the recent progress in understanding zeolite chemistry with the aid of these ssNMR methods and techniques, with a specific focus on the investigation of zeolite framework structures, zeolite crystallization mechanisms, surface active/acidic sites, host-guest/guest-guest interactions, and catalytic reaction mechanisms.

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.

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.

Regulation of the Si/Al ratios and Al distributions of zeolites and their impact on properties

Zeolites are typically a class of crystalline microporous aluminosilicates that are constructed by SiO₄ and AlO₄ tetrahedra. Because of their unique porous structures, strong Brönsted acidity, molecular-level shape selectivity, exchangeable cations, and high thermal/hydrothermal stability, zeolites are widely used as catalysts, adsorbents, and ion-exchangers in industry. The activity, selectivity, and stability/durability of zeolites in applications are closely related to their Si/Al ratios and Al distributions in the framework. In this review, we discussed the basic principles and the state-of-the-art methodologies for regulating the Si/Al ratios and Al distributions of zeolites, including seed-assisted recipe modification, interzeolite transformation, fluoride media, and usage of organic structure-directing agents (OSDAs), etc. The conventional and newly developed characterization methods for determining the Si/Al ratios and Al distributions were summarized, which include X-ray fluorescence spectroscopy (XRF), solid state ²⁹Si/²⁷Al magic-angle-spinning nuclear magnetic resonance spectroscopy (²⁹Si/²⁷Al MAS NMR), Fourier-transform infrared spectroscopy (FT-IR), etc. The impact of Si/Al ratios and Al distributions on the catalysis, adsorption/separation, and ion-exchange performance of zeolites were subsequently demonstrated. Finally, we presented a perspective on the precise control of the Si/Al ratios and Al distributions of zeolites and the corresponding challenges.

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.

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.

Recent advances in solid-state NMR of zeolite catalysts

Zeolites are important inorganic crystalline microporous materials with a broad range of applications in the areas of catalysis, ion exchange, and adsorption/separations. Solid-state nuclear magnetic resonance (NMR) spectroscopy has proven to be a powerful tool in the study of zeolites and relevant catalytic reactions because of its advantage in providing atomic-level insights into molecular structure and dynamic behavior. In this review, we provide a brief discussion on the recent progress in exploring framework structures, catalytically active sites and intermolecular interactions in zeolites and metal-containing ones by using various solid-state NMR methods. Advances in the mechanistic understanding of zeolite-catalysed reactions including methanol and ethanol conversions are presented as selected examples. Finally, we discuss the prospect of the solid-state NMR technique for its application in zeolites.

Regression Machine Learning Models Used to Predict DFT-Computed NMR Parameters of Zeolites

Machine learning approaches can drastically decrease the computational time for the predictions of spectroscopic properties in materials, while preserving the quality of the computational approaches. We studied the performance of kernel-ridge regression (KRR) and gradient boosting regressor (GBR) models trained on the isotropic shielding values, computed with density-functional theory (DFT), in a series of different known zeolites containing out-of-frame metal cations or fluorine anion and organic structure-directing cations. The smooth overlap of atomic position descriptors were computed from the DFT-optimised Cartesian coordinates of each atoms in the zeolite crystal cells. The use of these descriptors as inputs in both machine learning regression methods led to the prediction of the DFT isotropic shielding values with mean errors within 0.6 ppm. The results showed that the GBR model scales better than the KRR model.

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.

Hydroxyl environments in zeolites probed by deuterium solid-state MAS NMR combined with IR spectroscopy

Deuterium (2D) solid-state MAS NMR spectroscopy is used for the first time for identifying both silanols and Brønsted acid sites in zeolites. The environments of hydroxyl groups in faujasite-type zeolites...

A combination of Proton Spin Diffusion NMR and molecular simulations to probe supramolecular assemblies of organic molecules in nanoporous materials

In this work we show the use of high-resolution 1H MAS NMR to distinguish between two kinds of aggregation states of (1R,2S)-ephedrine, a chiral organic structure directing agent, occluded within...

Experimental and Theoretical Solid-State 29Si NMR Studies on Defect Structures in La9.33+x(SiO4)6O2+1.5x Apatite Oxide Ion Conductors

Oxide ion conductors can be used as electrolytes in solid oxide fuel cells, a promising energy-conversion technology. Local structures around the defects in oxide ion conductors are key for understanding the defect stabilization and migration mechanisms. As the defect contents are generally low, it is rather difficult to characterize the defect structure and therefore elucidate how oxide ions migrate. Solid-state nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for probing the local structures. However, the interpretation of NMR signals mainly based on the empirical knowledge could lead to unprecise local structures. There is still controversy regarding the defect structures in the apatite-type interstitial oxide ion conductors containing isolated tetrahedral units. Here, we combine the experimental solid-state ²⁹Si NMR spectroscopy with theoretical density functional theory calculations to investigate the defect structures in La_9.33+x(SiO₄)₆O_2+1.5x apatites. The results indicate that the ²⁹Si resonance signals on the high field side of the main peak corresponding to the Si atoms in the bulk structure are related to La vacancies and there is no steady-state SiO₅ in the defect structures. This finding provides new atomic-level understanding to the stabilization and migration of interstitial oxide ions in silicate apatites, which could guide the design and discovery of new solid oxide fuel cell electrolyte materials.

Silanol defect engineering and healing in zeolites: opportunities to fine-tune their properties and performances

Zeolites have been game-changing materials in oil refining and petrochemistry over the last 60 years and have the potential to play the same role in the emerging processes of the energy and environmental transition. Although zeolites are crystalline inorganic solids, their structures are not perfect and the presence of defect sites - mainly Brønsted acid sites and silanols - influences their thermal and chemical resistance as well as their performances in key areas such as catalysis, gas and liquid separations and ion-exchange. In this paper, we review the type of defects in zeolites and the characterization techniques used for their identification and quantification with the focus on diffraction, spectroscopic and modeling approaches. More specifically, throughout the review, we will focus on silanol (Si-OH) defects located within the micropore structure and/or on the external surface of zeolites. The main approaches applied to engineer and heal defects and their consequences on the properties and applications of zeolites in catalysis and separation processes are highlighted. Finally, the challenges and opportunities of silanol defect engineering in tuning the properties of zeolites to meet the requirements for specific applications are presented.

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.

The accuracy challenge of the DFT-based molecular assignment of 13C MAS NMR characterization of surface intermediates in zeolite catalysis

The influence of the model and method choice on the DFT predicted ¹³C NMR chemical shifts of zeolite surface methoxide species has been systematically analyzed. Twelve ¹³C NMR chemical shift calculation protocols on full periodic and hybrid periodic-cluster DFT calculations with varied structural relaxation procedures are examined. The primary assessment of the accuracy of the computational protocols has been carried out for the Si-O(CH₃)-Al surface methoxide species in ZSM-5 zeolite with well-defined experimental NMR parameters (chemical shift, δ(¹³C) value) as a reference. Different configurations of these surface intermediates and their location inside the ZSM-5 pores are considered explicitly. The predicted δ value deviates by up to ±0.8 ppm from the experimental value of 59 ppm due to the varied confinement of the methoxide species at different zeolite sites (model accuracy). The choice of the exchange-correlation functional (method accuracy) introduces ±1.5 ppm uncertainty in the computed chemical shifts. The accuracy of the predicted ¹³C NMR chemical shifts for the computational assignment of spectral characteristics of zeolite intermediates has been further analyzed by considering the potential intermediate species formed upon methane activation by Cu/ZSM-5 zeolite. The presence of Cu species in the vicinity of surface methoxide increases the prediction uncertainty to ±2.5 ppm. The full geometry relaxation of the local environment of an active site at an appropriate level of theory is critical to ensure a good agreement between the experimental and computed NMR data. Chemical shifts (δ) calculated via full geometry relaxation of a cluster model of a relevant portion of the zeolite lattice site are in the best agreement with the experimental values. Our analysis indicates that the full geometry optimization of a cluster model at the PBE0-D3/6-311G(d,p) level of theory followed by GIAO/PBE0-D3/aug-cc-pVDZ calculations is the most suitable approach for the calculation of ¹³C chemical shifts of zeolite surface intermediates.

Single-Atom Catalysts across the Periodic Table

Isolated atoms featuring unique reactivity are at the heart of enzymatic and homogeneous catalysts. In contrast, although the concept has long existed, single-atom heterogeneous catalysts (SACs) have only recently gained prominence. Host materials have similar functions to ligands in homogeneous catalysts, determining the stability, local environment, and electronic properties of isolated atoms and thus providing a platform for tailoring heterogeneous catalysts for targeted applications. Within just a decade, we have witnessed many examples of SACs both disrupting diverse fields of heterogeneous catalysis with their distinctive reactivity and substantially enriching our understanding of molecular processes on surfaces. To date, the term SAC mostly refers to late transition metal-based systems, but numerous examples exist in which isolated atoms of other elements play key catalytic roles. This review provides a compositional encyclopedia of SACs, celebrating the 10th anniversary of the introduction of this term. By defining single-atom catalysis in the broadest sense, we explore the full elemental diversity, joining different areas across the whole periodic table, and discussing historical milestones and recent developments. In particular, we examine the coordination structures and associated properties accessed through distinct single-atom-host combinations and relate them to their main applications in thermo-, electro-, and photocatalysis, revealing trends in element-specific evolution, host design, and uses. Finally, we highlight frontiers in the field, including multimetallic SACs, atom proximity control, and possible applications for multistep and cascade reactions, identifying challenges, and propose directions for future development in this flourishing field.

Local anisotropy in single crystals of zeotypes with the MFI framework structure evidenced by polarised Raman spectroscopy

Polarised Raman spectroscopy is used to characterise the local structure in single crystals of zeotypes, namely silicalite-1 and ZSM-5, which share the MFI framework structure. Attributes favourable for applying polarised Raman spectroscopy are the orthogonal axes of these single crystals and their size, i.e. 10 to 30 micrometers in all three directions. We show that the intensity of certain vibrational modes in silicalite-1 depends on the polarisation of the incident light, reflecting the anisotropic character of the molecular bonds contributing to these vibrations. Using these observations, and by estimating the depolarisation ratio (ρ) and the pseudo-order factor (f), we propose a more accurate assignment of the Raman active modes. More precisely, Raman intensities peaked at 294, 360, 383 and 472 cm^-1 are attributed to bending modes in 10-, 6-, 5- and 4-membered rings, respectively. In the region of stretching modes, the vibration at 832 cm^-1 is assigned to Si-O-Si bonds shared between 5-membered rings, which have an orientation parallel to the a-axis of the crystal. By virtue of having a strongly polarised character, the modes at 472 and 832 cm^-1 can be used as orientational indicators. The proposed assignment is supported by the good agreement between experimental and simulated polar plots, where Raman intensities are plotted as a function of the polarisation angle of the incident light. Finally, upon partial substitution of Si atoms by Al, the crystalline structure is maintained and almost no spectroscopic changes are observed. The only significant difference is the increased width of most vibrational modes, which is consistent with the local lower symmetry. This is also seen in the angular dependence of selected vibrational modes that compared to the case of pure silicalite-1 appear less polarised. In the Raman spectrum of ZSM-5 a new feature at 974 cm^-1 is observed, which we attribute to Al-OH stretching. In the high frequency range, the O-H stretching modes are observed which arise from the Si-O(H)-Al Brønsted acid sites. The intensity of the characteristic mode at 3611 cm^-1 reveals an anisotropic character as well, which is in line with previous findings from solid state NMR that Al atoms distribute nonrandomly within the MFI framework structure.

Zeolite Structure Direction: Identification, Strength and Involvement of Weak CH⋅⋅⋅O Hydrogen Bonds

We demonstrate that weak CH⋅⋅⋅O hydrogen bonds (HBs) are important host-guest interactions in zeolite assemblies involving structure directing organocations. This type of HB is identified between alkyl groups of the organic structure directing agent (OSDA) and the silica framework in as-synthesized silicalite-1 of complex topology (MFI) using a combination of experimental and theoretical data obtained at low and room temperatures. The 28 weak CH⋅⋅⋅O HBs, evidenced along dynamics simulation at room temperature, represent 30 % of the energy of the Coulomb electrostatic interaction between OSDA and the zeolite framework. The strongest and most stable HB found here connects the OSDA to the [4¹ 5² 6² ] cage containing F atoms and should contribute to preserve zeolite topology during crystal growth. An inspection of other as-synthesized zeolites of very different framework topology indicates that the directional CH⋅⋅⋅O HBs have to be considered when discussing zeolite structure directing phenomena.

A Supramolecular View on the Cooperative Role of Brønsted and Lewis Acid Sites in Zeolites for Methanol Conversion

A systematic molecular level and spectroscopic investigation is presented to show the cooperative role of Brønsted acid and Lewis acid sites in zeolites for the conversion of methanol. Extra-framework alkaline-earth metal containing species and aluminum species decrease the number of Brønsted acid sites, as protonated metal clusters are formed. A combined experimental and theoretical effort shows that postsynthetically modified ZSM-5 zeolites, by incorporation of extra-framework alkaline-earth metals or by demetalation with dealuminating agents, contain both mononuclear [MOH]⁺ and double protonated binuclear metal clusters [M(μ-OH)₂M]²⁺ (M = Mg, Ca, Sr, Ba, and HOAl). The metal in the extra-framework clusters has a Lewis acid character, which is confirmed experimentally and theoretically by IR spectra of adsorbed pyridine. The strength of the Lewis acid sites (Mg > Ca > Sr > Ba) was characterized by a blue shift of characteristic IR peaks, thus offering a tool to sample Lewis acidity experimentally. The incorporation of extra-framework Lewis acid sites has a substantial influence on the reactivity of propene and benzene methylations. Alkaline-earth Lewis acid sites yield increased benzene methylation barriers and destabilization of typical aromatic intermediates, whereas propene methylation routes are less affected. The effect on the catalytic function is especially induced by the double protonated binuclear species. Overall, the extra-framework metal clusters have a dual effect on the catalytic function. By reducing the number of Brønsted acid sites and suppressing typical catalytic reactions in which aromatics are involved, an optimal propene selectivity and increased lifetime for methanol conversion over zeolites is obtained. The combined experimental and theoretical approach gives a unique insight into the nature of the supramolecular zeolite catalyst for methanol conversion which can be meticulously tuned by subtle interplay of Brønsted and Lewis acid sites.

Magnetization transfer from protons to quadrupolar nuclei in solid-state NMR using PRESTO or dipolar-mediated refocused INEPT methods

In solid-state NMR spectroscopy, the through-space transfer of magnetization from protons to quadrupolar nuclei is employed to probe proximities between those isotopes. Furthermore, such transfer, in conjunction with Dynamic Nuclear Polarization (DNP), can enhance the NMR sensitivity of quadrupolar nuclei, as it allows the transfer of DNP-enhanced ¹H polarization to surrounding nuclei. We compare here the performances of two approaches to achieve such transfer: PRESTO (Phase-shifted Recoupling Effects a Smooth Transfer of Order), which is currently the method of choice to achieve the magnetization transfer from protons to quadrupolar nuclei and which has been shown to supersede Cross-Polarization under Magic-Angle Spinning (MAS) for quadrupolar nuclei and D-RINEPT (Dipolar-mediated Refocused Insensitive Nuclei Enhanced by Polarization Transfer) using symmetry-based SR4₁² recoupling, which has already been employed to transfer the magnetization in the reverse way from half-integer quadrupolar spin to protons. We also test the PRESTO sequence with R16₇⁶ recoupling using 270₀90₁₈₀ composite π-pulses as inversion elements. This recoupling scheme, which has previously been proposed to reintroduce ¹H Chemical Shift Anisotropy (CSA) at high MAS frequencies with high robustness to rf-field inhomogeneity, has not so far been employed to reintroduce dipolar couplings with protons. These various techniques to transfer magnetization from protons to quadrupolar nuclei are analyzed using (i) an average Hamiltonian theory, (ii) numerical simulations of spin dynamics, and (iii) experimental ¹H → ²⁷Al and ¹H → ¹⁷O transfers in as-synthesized AlPO₄-14 and ¹⁷O-labelled fumed silica, respectively. The experiments and simulations are done at two magnetic fields (9.4 and 18.8 T) and several spinning speeds (15, 18-24 and 60 kHz). This analysis indicates that owing to its γ-encoded character, PRESTO yields the highest transfer efficiency at low magnetic fields and MAS frequencies, whereas owing to its higher robustness to rf-field inhomogeneity and chemical shifts, D-RINEPT is more sensitive at high fields and MAS frequencies, notably for protons exhibiting large offset or CSA, such as those involved in hydrogen bonds.

Active sites and mechanism of the direct conversion of methane and carbon dioxide to acetic acid over the zinc-modified H-ZSM-5 zeolite

The catalytic activity of the conversion of CH₄ and CO₂ on zinc modified H-ZSM-5 is strongly dependent on the structure of the active sites.

Determining the structures, acidity and adsorption properties of Al substituted HZSM-5

Determining the locations and distributions of Al substitution in zeolite-based catalysts and catalysis is always very challenging. Despite advanced experimental characterization techniques and improved theoretical models, this issue is not reasonably solved and this is because the locations and distributions of Al substitution in zeolites are more kinetically than thermodynamically controlled. In this work, we computed one Al substitution in the orthorhombic form of MFI (HZSM-5) which contains 12 distinct tetrahedral (T) centers on the basis of a periodic slab model containing 96 T centers including van der Waals dispersion correction (GGA-PBE-D3). For all 12 T centers, there are 48 acidic sites and each site can be considered for the adsorption of probe molecules. Thermodynamically, the energy span of the twelve most stable acidic sites is less than 15 kJ mol-1, and such a small energy difference enables all adjustable possibilities for the locations and distributions of Al substitution under suitable conditions. Excellent agreement between experiment and theory in the adsorption enthalpies of pyridine, methylamine, dimethylamine and trimethylamine shows that the location of Al substitution is most likely at T1, T3, T5, T7 and T11, while much less likely at the often used T12 site. These results provide the basis for identifying Al substitution in new synthesized HZSM-5 catalysts and for studying the acidic site-dependent catalytic activity of HZSM-5 in cracking and hydrogenation reactions.