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

Assignment of IR spectra of ethanol at Brønsted sites of H-ZSM-5 to monomer adsorption using a Fermi resonance model

Understanding how alcohol molecules interact with the Brønsted acid sites (BAS) of zeolites is a prerequisite to the design of zeolite catalysts and catalytic processes. Here, we report IR spectra for the adsorption of ethanol on a highly crystalline sample of H-ZSM-5 zeolites exposed to ethanol gas at increasing pressure. We use density functional theory in combination with a FERMI resonance model to assign the measured spectra to a single adsorbed ethanol molecule per BAS. Specifically, we assign the bands at 2450 cm-1 and 1670 cm-1 to a FERMI resonance between the fundamental (Z)O-H stretching band of a single-ethanol-loaded BAS and the first overtone of the (Z)O-H out-of-plane bending. We conclude that adsorbed dimers do not contribute in a noticeable way up to a concentration of almost one ethanol molecule per BAS site. We further show that hybrid functionals (B3LYP) are required to get a close match between the predicted and experimental spectra, whereas commonly used generalized gradient type functionals such as PBE incorrectly describe the potential energy surface. They overestimate the redshift of the OH stretching band on hydrogen bond formation which results in an erroneous assignment of the IR bands.

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.

Understanding deNOx mechanisms in transition metal exchanged zeolites

Transition-metal-containing zeolites have wide-ranging applications in several catalytic processes including the selective catalytic reduction (SCR) of NOx species. To understand how transition metal ions (TMIs) can effect NOx reduction chemistry, both structural and mechanistic aspects at the atomic level are needed. In this review, we discuss the coordination chemistry of TMIs and their mobility within the zeolite framework, the reactivity of active sites, and the mechanisms and intermediates in the NH3-SCR reaction. We emphasise the key relationship between TMI coordination and structure and mechanism and discuss approaches to enhancing catalytic activity via structural modifications.

NMR Relaxation to Probe Zeolites: Mobility of Adsorbed Molecules, Surface Acidity, Pore Size Distribution and Connectivity

Unique structural and chemical properties, such as ion exchange, developed inner surface, etc., as well as the wide possibilities and flexibility of regulating these properties, cause a keen interest in zeolites. They are widely used in industry as molecular sieves, ion exchangers and catalysts. Current trends in the development of zeolite-based catalysts include the adaptation of their cationic composition, acidity and porosity for a specific catalytic process. Recent studies have shown that mesoporosity is beneficial to the rational design of catalysts with controlled product selectivity and an improved catalyst lifetime due to its efficient mass-transport properties. Nuclear magnetic resonance (NMR) has proven to be a reliable method for studying zeolites. Solid-state NMR spectroscopy allows for the quantification of both Lewis and Brønsted acidity in zeolite catalysts and, nowadays, 27Al and 29Si magic angle spinning NMR spectroscopy has become firmly established in the set of approved methods for characterizing zeolites. The use of probe molecules opens up the possibility for the indirect measurement of the characteristics of acid sites. NMR relaxation is less common, although it is especially informative and enlightening for studying the mobility of guest molecules in the porous matrix. Moreover, the NMR relaxation of guest molecules and NMR cryoporometry can quantify pore size distribution on a broader scale (compared to traditional methods), which is especially important for systems with complex pore organization. Over the last few years, there has been a growing interest in the use of 2D NMR relaxation techniques to probe porous catalysts, such as 2D T1-T2 correlation to study the acidity of the surface of catalysts and 2D T2-T2 exchange to study pore connectivity. This contribution provides a comprehensive review of various NMR relaxation techniques for studying porous media and recent results of their applications in probing micro- and mesoporous zeolites, mainly focused on the mobility of adsorbed molecules, the acidity of the zeolite surface and the pore size distribution and connectivity of zeolites with hierarchical porosity.

Unravelling the Crucial of Spatial Al Distribution to Realize Precise Alkali-Treatment for Target Acid-Catalyzed Reactions

Constructing mesoporous structure within zeolites by alkali-treatment is an effective protocol to improve their diffusion properties. However, undesirable changes in Brönsted acid site (BAS) densities always offset this advantage in acid-catalyzed reactions. In this context, the crucial roles of spatial aluminum (Al) distribution were unraveled during alkali-treatment of MFI zeolite and the desirable BAS density was achieved in obtained hierarchical samples for the target reactions. Various characterization methods, particularly the multiple one- and two-dimensional magic-angle spinning (MAS) NMR techniques, were performed to track the alkali-treatment processes. For the sample with a more uniform spatial Al distribution, more tetrahedral Al sites would fall off and migrate around the Si-OH in zeolite as Al(OH)4 -. Those re-deposited Al(OH)4 - sites were easily transformed into NMR-invisible Al sites during the calcination process, which contributed negligibly to both Brönsted and Lewis acidities, thus being referred to "acid-free" Al species. While most tetrahedral Al sites were preserved after the alkali-treatment of sample with non-uniform Al distribution and the BAS density gradually increased with treatment time. According to the requirements of typical acid-catalyzed reactions, such as catalytic cracking of 1,3,5-triisopropylbenzene and methanol-to-olefins, the desired hierarchical zeolite catalysts were developed by matching the amounts of extracted Si and generated "acid-free" Al during the precise alkali-treatment.

Strength of London Dispersion Forces in Organic Structure Directing Agent-Zeolite Assemblies

Herein, we study the London dispersion forces between organic structure directing agents (OSDAs)-here tetraalkyl-ammonium or -phosphonium molecules-and silica zeolite frameworks (FWs). We demonstrate that the interaction energy for these dispersion forces is correlated to the number of H atoms in OSDAs, irrespective of the structures of OSDAs or FWs, and of variations in charges and thermal motions. All calculations considered-DFT-D3 and BOMD undertaken by us, and molecular mechanics from an accessible database-led to the same trend. The mean energy of these dispersion forces is ca. -2 kcal.mol-1 per H for efficient H-O contacts.

Engineered species-selective ion-exchange in tuneable dual-phase zeolite composites

Controllable sorption selectivity in zeolites is crucial for their application in catalysis, gas separation and ion-exchange. Whilst existing approaches to achieving sorption selectivity with natural zeolites typically rely on screening for specific geological deposits, here we develop partial interzeolite transformation as a straightforward and highly tuneable method to achieve sorption selectivity via forming dual-phase composites with simultaneous control of both phase-ratio and morphology. The dual-cation (strontium and caesium) exchange properties of a series of granular mordenite/zeolite P composites formed from a parent natural mordenite material are demonstrated in complex, industrially relevant multi-ion environments pertinent to nuclear waste management. The relative uptake of caesium and strontium is controlled via the extent of transformation: composites exhibit significantly increased ion-exchange affinity for strontium compared to both the parent mordenite and physical mixtures of mordenite/zeolite P phases with similar phase ratios. The composite with a 40 : 60 mordenite : zeolite P ratio composite achieves higher uptake rates than the natural clinoptilolite material currently used to decontaminate nuclear waste streams at the Sellafield site, UK. In situ X-ray image-guided diffraction experiments during caesium exchange demonstrate that the mordenite core retains rapid caesium uptake likely responsible for the unique ion-exchange chemistry achievable through the partial inter-zeolite transformation. These results offer a straightforward and controllable route to optimised zeolite functionality and a strategy to engineer composites from low-grade natural sources at low cost and with formulation advantages for industrial deployment.

A "Liquid-In-Solid" Electrolyte for High-Voltage Anode-Free Rechargeable Sodium Batteries

Developing anode-free batteries is the ultimate goal in pursuit of high energy density and safety. It is more urgent for sodium (Na)-based batteries due to its inherently low energy density and safety hazards induced by highly reactive Na metal anodes. However, there is no electrolyte that can meet the demanding Na plating-stripping Coulomb efficiency (CE) while resisting oxidative decomposition at high voltages for building stable anode-free Na batteries. Here, a "liquid-in-solid" electrolyte design strategy is proposed to integrate target performances of liquid and solid-state electrolytes. Breaking through the Na+ transport channel of Na-containing zeolite molecular sieve by ion-exchange and confining aggregated liquid ether electrolytes in the nanopore and void of zeolites, it achieves excellent high-voltage stability enabled by solid-state zeolite electrolytes, while inheriting the ultra-high CE (99.84%) from liquid ether electrolytes. When applied in a 4.25 V-class anode-free Na battery, an ultra-high energy density of 412 W h kg-1 (based on the active material of both cathodes and anodes) can be reached, which is comparable to the state-of-the-art graphite||LiNi0.8Co0.1Mn0.1O2 lithium-ion batteries. Furthermore, the assembled anode-free pouch cell exhibits excellent cycling stability, and a high capacity retention of 89.2% can be preserved after 370 cycles.

Real-Space Pseudopotential Method for the Calculation of Third-Row Elements X-ray Photoelectron Spectroscopic Signatures

X-ray photoelectron spectroscopy (XPS) is a powerful characterization technique that unveils subtle chemical environment differences via core-electron binding energy (CEBE) analysis. We extend the development of real-space pseudopotential methods to calculating 1s, 2s, and 2p3/2 CEBEs of third-row elements (S, P, and Si) within the framework of Kohn-Sham density-functional theory (KS-DFT). The new approach systematically prevents variational collapse and simplifies core-excited orbital selection within dense energy level distributions. However, careful error cancellation analysis is required to achieve accuracy comparable to all-electron methods and experiments. Combined with real-space KS-DFT implementation, this development enables large-scale simulations with both Dirichlet boundary conditions and periodic boundary conditions.

Synergy between Brønsted Acid Sites and Carbonaceous Deposits during Skeletal 1-Butene Isomerization over Ferrierite

During skeletal 1-butene isomerization over ferrierite carbonaceous deposits block 98% of the micropores within 24 h, rendering them effectively inaccessible to reactants, while the catalytic activity improves continuously for 100 h on stream. Ex-situ pyridine adsorption shows that the concentration of conventional Brønsted acid sites in the 10-R channels decreases below the detection threshold of infrared spectroscopy within 2 h. However, the operando addition of the base triethyl amine to the feed quenches the reaction, showing that mediated acidity is necessary. The larger base 2,2,6,6-tetramethyl piperidine only deactivates catalytic activity after several hours because it cannot directly bind to active sites at the sterically restricted pore mouths. The communication of internal Brønsted acid sites to the external reactants via a concerted mechanism involving protonated monoaromatic deposits trapped in the pore mouths explains the promoting effects of coke species in zeolite-catalyzed skeletal butene isomerization. This work presents a consolidated explanation of the synergy of solid acidity, structural confinement, and carbonaceous deposits in zeolites.

Construction of diverse hollow MFI zeolites through regulating the micropore filling agents

Constructing hollow structure into microporous zeolites can improve the accessibility of acid sites located at the inner part and the diffusion property. Hence, the development of an efficient synthesis strategy to acquire zeolites with tunable hollow structures and acidity has attracted much attention. In this work, an innovative tandem synthesis route was proposed to prepare MFI zeolites with diverse hollow structure while maintaining solid yields exceeding 90 %. The substitution of ethanol molecules, which previously occupied the micropores, with tetrapropylammonium cations was proved to be the key factor to construct hollow structure. And a crystallization-driven particle dissolution mechanism was proposed. The dimension of the hollow cavity, particle size, and Si/Al ratio can be flexibly regulated. Interestingly, hollow MFI samples featuring the common cavity structure, "eye-like" cavity structure, or double-cavity structure can be directly synthesized by controlling the dissolution of core parts. In the 1-butene catalytic cracking reactions, a much higher conversion of 67.2 % was acquired over hollow ZSM-5 compared with that over conventional ZSM-5 (35.8 %) after 64 h of reaction. This improvement can be attributed to the eightfold increase of diffusivity in hollow ZSM-5. This facile and efficient synthesis method endows accurate regulation of the hollow structure, which is meaningful for both fundamental research and industrial applications of hollow zeolites.

A reactive neural network framework for water-loaded acidic zeolites

Under operating conditions, the dynamics of water and ions confined within protonic aluminosilicate zeolite micropores are responsible for many of their properties, including hydrothermal stability, acidity and catalytic activity. However, due to high computational cost, operando studies of acidic zeolites are currently rare and limited to specific cases and simplified models. In this work, we have developed a reactive neural network potential (NNP) attempting to cover the entire class of acidic zeolites, including the full range of experimentally relevant water concentrations and Si/Al ratios. This NNP has the potential to dramatically improve sampling, retaining the (meta)GGA DFT level accuracy, with the capacity for discovery of new chemistry, such as collective defect formation mechanisms at the zeolite surface. Furthermore, we exemplify how the NNP can be used as a basis for further extensions/improvements which include data-efficient adoption of higher-level (hybrid) references via Δ-learning and the acceleration of rare event sampling via automatic construction of collective variables. These developments represent a significant step towards accurate simulations of realistic catalysts under operando conditions.

Unraveling Spatially Dependent Hydrophilicity and Reactivity of Confined Carbocation Intermediates during Methanol Conversion over ZSM-5 Zeolite

Carbocations play a pivotal role as reactive intermediates in zeolite-catalyzed methanol-to-hydrocarbon (MTH) transformations. However, the interaction between carbocations and water vapor and its subsequent effects on catalytic performance remain poorly understood. Using micro-magnetic resonance imaging (μMRI) and solid-state NMR techniques, this work investigates the hydrophilic behavior of cyclopentenyl cations within ZSM-5 pores under vapor conditions. We show that the polar cationic center of cyclopentenyl cations readily initiates water nucleus formation through water molecule capture. This leads to an inhomogeneous water adsorption gradient along the axial positions of zeolite, correlating with the spatial distribution of carbocation concentrations. The adsorbed water promotes deprotonation and aromatization of cyclopentenyl cations, significantly enhancing the aromatic product selectivity in MTH catalysis. These results reveal the important influence of adsorbed water in modulating the carbocation reactivity within confined zeolite pores.

Catalytic performance and mechanism study of the isomerization of 2,5-dichlorotoluene to 2,4-dichlorotoluene

This work investigates the influence of catalyst HZSM-5 on the isomerization of 2,5-dichlorotoluene (2,5-DCT) to produce 2,4-dichlorotoluene (2,4-DCT). We observe that hydrothermal treatment leads to a decrease in total acidity and Brønsted/Lewis ratio of HZSM-5 while generating new secondary pores. These characteristics result in excellent selectivity for post-hydrothermal modified HZSM-5 in the isomerization reaction from 2,5-DCT to 2,4-DCT. Under atmospheric pressure at 350 °C, unmodified HZSM-5 achieves a selectivity of 66.4% for producing 2,4-DCT, however after hydrothermal modification the selectivity increases to 78.7%. Density Functional Theory (DFT) calculations explore the thermodynamic aspects of adsorption between the HZSM-5 surface and 2,4-DCT. The kinetic perspective investigates the mechanism involving proton attack on the methyl group of 2,5-DCT followed by rearrangement leading to formation of 2,4-DCT during isomerization. The consistency between simulation and experimental results provides evidence for the feasibility of isomerizing 2,5-DCT to 2,4-DCT. This work fills the gap in the low value-added product 2,5-DCT isomer conversion, indicating its significant practical application potential and provides a valuable reference and guidelines for industrial research in this field.

Unveiling the mechanisms of carboxylic acid esterification on acid zeolites for biomass-to-energy: A review of the catalytic process through experimental and computational studies

In recent years, there has been significant interest from industrial and academic areas in the esterification of carboxylic acids catalyzed by acidic zeolites, as it represents a sustainable and economically viable approach to producing a wide range of high-value-added products. However, there is a lack of comprehensive reviews that address the intricate reaction mechanisms occurring at the catalyst interface at both the experimental and atomistic levels. Therefore, in this review, we provide an overview of the esterification reaction on acidic zeolites based on experimental and theoretical studies. The combination of infrared spectroscopy with atomistic calculations and experimental strategies using modulation excitation spectroscopy techniques combined with phase-sensitive detection is presented as an approach to detecting short-lived intermediates at the interface of zeolitic frameworks under realistic reaction conditions. To achieve this goal, this review has been divided into four sections: The first is a brief introduction highlighting the distinctive features of this review. The second addresses questions about the topology and activity of different zeolitic systems, since these properties are closely correlated in the esterification process. The third section deals with the mechanisms proposed in the literature. The fourth section presents advances in IR techniques and theoretical calculations that can be applied to gain new insights into reaction mechanisms. Finally, this review concludes with a subtle approach, highlighting the main aspects and perspectives of combining experimental and theoretical techniques to elucidate different reaction mechanisms in zeolitic systems.

Rare Earth Elements Recovery and Waste Management of Municipal Solid Waste Incineration Ash

The advancements in high-tech products and pursuit of renewable energy demand a massive and continuously growing supply of rare earth elements (REE). However, REE production from mining is heavily restricted by technoeconomic limitations and global geopolitical tensions. Municipal solid waste incineration ash (MSWIA) has been recently recognized as a potential alternative for REE recovery. This study applies and optimizes a green modular treatment system using organic ligands for effective REE recovery and concentration from MSWIA with minimal generation of secondary wastes. Citrate extracted >80% of total REE at pH 2.0 and ∼60% at pH 4.0. A subsequent oxalate precipitation step selectively concentrated >98% of extracted REE by ∼7-12 times compared to raw MSWIA. Waste byproducts were upcycled to synthesize zeolites, resulting in an overall solid waste volume reduction of ∼80% and heavy metal immobilization efficiency of ∼75% with negligible leaching, bringing the dual benefits of REE recovery and waste management. This work serves as a pioneer study in REE recovery from an emerging source and provides system level insights on the practicality of a simple three-step treatment system. Compared to existing literature, this system features a low chemical/energy input and a light environmental footprint.

Variation of the Orientations of Organic Structure-Directing Agents inside the Channels of SCM-14 and SCM-15 Germanosilicates Obtained by Ab Initio Molecular Dynamic Simulations

We report ab initio molecular dynamic simulations of the organic structure-directing agent (OSDA) in the channels of SCM-14 and SCM-15 germanosilicates for models with different germanium distribution. Since OSDA was free to move inside the channels, independent of its initial orientation after the simulations in all structures the OSDA, protonated 4-pyrrolidinopyridine, is positioned almost perpendicular to the large channels of SCM-14. The structures obtained from the dynamic simulation are more stable by 157 to 331 kJ/mol than the structures obtained by initial geometry optimization. After simulations, the average distance between the N atom of the pyridine moiety of the OSDA and O from Ge-O-Ge is shorter by 0.2 Å than the same distance obtained from initial optimization. The stretching N-H frequencies in the IR spectra of the OSDA and other calculated vibrational frequencies are not characteristic of the orientation of the molecule and cannot be used to detect it.

Revealing the evolution of local structures in the formation process of alkaline earth metal cation-containing zeolites from glasses

Alkaline earth metal cations are ubiquitously present in natural zeolites but less exploited in synthetic zeolites due to their low solubility in water, and hence it remains elusive how they contribute to zeolite formation. Herein, harmotome, a PHI-type zeolite with Ba2+, is readily synthesized from a Ba-containing aluminosilicate glass. This glass-to-zeolite transformation process, in particular the structure-regulating role of Ba2+, is investigated by anomalous X-ray scattering and high-energy X-ray total scattering techniques. The results demonstrate that the steady Ba2+-aluminosilicate interactions not only help prevent the precipitation of barium species under alkaline synthetic conditions, but also dictate the local structures with distinct interatomic distances between the Ba2+ and the surrounding aluminosilicate species throughout the transformation process, which lead to the successful formation of harmotome without detectable impurities. This study highlights the usefulness of the comprehensive X-ray scattering techniques in revealing the formation scheme of the zeolites containing specific metal species. In addition, a promising alternative approach to design and synthesize zeolites with unique compositions and topologies by using well-crafted glasses with suitable metal cation dopants is demonstrated.

Ab Initio Screening of Divalent Cations for CH4, CO2, H2, and N2 Separations in Chabazite Zeolite

The efficient separation and adsorption of critical gases are, more than ever, a major focus point in important energy processes, such as CH4 enrichment of biogas or natural gas, CO2 separation and capture, and H2 purification and storage. Thanks to its physicochemical properties, cation-exchanged chabazite is a potent zeolite for such applications. Previous computational screening investigations have mostly examined chabazites exchanged with monovalent cations. Therefore, in this contribution, periodic density functional theory (DFT) calculations in combination with dispersion corrections have been used for a systematic screening of divalent cation exchanged chabazite zeolites. The work focuses on cheap and readily available divalent cations, Ca(II), Mg(II), and Zn(II), Fe(II), Sn(II), and Cu(II) and investigates the effect of the cation nature and location within the framework on the adsorption selectivity of chabazite for specific gas separations, namely, CO2/CH4, N2/CH4, and N2/H2. All the cationic adsorption sites were explored to describe the diversity of sites in a typical experimental chabazite with a Si/Al ratio close to 2 or 3. The results revealed that Mg-CHA is the most promising cation for the selective adsorption of CO2. These predictions were further supported by ab initio molecular dynamics simulations performed at 300 K, which demonstrated that the presence of CH4 has a negligible impact on the adsorption of CO2 on Mg-CHA. Ca(II) was found to be the most favorable cation for the selective adsorption of H2 and CO2. Finally, none of the investigated cations were suitable for the preferential capture of N2 and H2 in the purification of CH4 rich mixtures. These findings provide valuable insights into the factors influencing the adsorption behavior of N2, H2, CH4, and CO2 and highlight the crucial role played by theoretical calculations and simulations for the optimal design of efficient adsorbents.

Bridging the Gap between Zeolites and Dense Silica Polymorphs: Formation of All-Silica Zeolite with High Framework Density from Natural Layered Silicate Magadiite

A silica zeolite (RWZ-1) with a very high framework density (FD) was synthesized from highly crystalline natural layered silicate magadiite, bridging the gap between the two research areas of zeolites and dense silica polymorphs. Magadiite was topotactically converted into a 3D framework through two-step heat treatment. The resulting structure had a 1D micropore system of channel-like cavities with an FD of 22.1 Si atoms/1000 Å3 . This value is higher than those of all other silica zeolites reported so far, approaching those of silica polymorphs (tridymite (22.6) and α-quartz (26.5)). RWZ-1 is a slight negative thermal expansion material with thermal properties approaching those of dense silica polymorphs. It contributes to the creation of a new field on microporous high-density silica/silicates. Synergistic interactions are expected between the micropores with molecular sieving properties and the dense layer-like building units with different topologies which provide thermal and mechanical stabilities.

Methanol diffusion and dynamics in zeolite H-ZSM-5 probed by quasi-elastic neutron scattering and classical molecular dynamics simulations

Zeolite ZSM-5 is a key catalyst in commercially relevant processes including the widely studied methanol to hydrocarbon reaction, and molecular diffusion in zeolite pores is known to be a crucial factor in controlling catalytic reactions. Here, we present critical analyses of recent quasi-elastic neutron scattering (QENS) data and complementary molecular dynamics (MD) simulations. The QENS experiments show that the nature of methanol diffusion dynamics in ZSM-5 pores is dependent both on the Si/Al ratio (11, 25, 36, 40 and 140), which determines the Brønsted acid site density of the zeolite, and that the nature of the type of motion observed may vary qualitatively over a relatively small temperature range. At 373 K, on increasing the ratio from 11 to 140, the observed mobile methanol fraction increases and the nature of methanol dynamics changes from rotational (in ZSM-5 with Si/Al of 11) to translational diffusion. The latter is either confined localized diffusion within a pore in zeolites with ratios up to 40 or non-localized, longer-range diffusion in zeolite samples with the ratio of 140. The complementary MD simulations conducted over long time scales (1 ns), which are longer than those measured in the present study by QENS (≈1-440 ps), at 373 K predict the occurrence of long-range translational diffusion of methanol in ZSM-5, independent of the Si/Al ratios (15, 47, 95, 191 and siliceous MFI). The rate of diffusion increases slightly by increasing the ratio from 15 to 95 and thereafter does not depend on zeolite composition. Discrepancies in the observed mobile methanol fraction between the MD simulations (100% methanol mobility in ZSM-5 pores across all Si/Al ratios) and QENS experiments (for example, ≈80% immobile methanol in ZSM-5 with Si/Al of 11) are attributed to the differences in time resolutions of the techniques. This perspective provides comprehensive information on the effect of acid site density on methanol dynamics in ZSM-5 pores and highlights the complementarity of QENS and MD, and their advantages and limitations. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 2)'.

Exploring the length scales, timescales and chemistry of challenging materials (Part 2)

This themed issue explores the different length and timescales that determine the physics and chemistry of a variety of key of materials, explored from the perspective of a wide range of disciplines, including physics, chemistry materials science, Earth science and biochemistry. The topics discussed include catalysis, chemistry under extreme conditions, energy materials, amorphous and liquid structure, hybrid organic materials and biological materials. The issue is in two parts, with this second set of contributions exploring hybrid organic materials, catalysis low-dimensional and graphitic materials, biological materials and naturally occurring, super-hard material as well as dynamic high pressure and new developments in imaging techniques pressure. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 2)'.

Potential of nanostructured carbon materials for iodine detection in realistic environments revealed by first-principles calculations

In the context of effective detection of iodine species (I2, CH3I) formed in nuclear power plants and nuclear fuel reprocessing facilities, we perform a comparative study of the potential sensing performance of four expectedly promising 2D materials (8-Pmmn borophene, BC3, C3N, and BC6N) towards the iodine-containing gases and, with the view of checking selectivity, towards common inhibiting gases in the containment atmosphere (H2O and CO), applying methods of dispersion-corrected density functional theory with periodic boundary conditions. A covalent bond is formed between the CO molecule and boron in BC3 or in 8-Pmmn borophene, compromising the anticipated applicability of these materials for iodine detection. The presence of nitrogen atoms in BC6N-2 prevents the formation of a covalent bond with CO; however, the closeness of adsorption energies for all the four gases studied does not distinguish this material as specifically sensitive to iodine species. Finally, the energies of adsorption on C3N yield a significant and promising discrimination between the adsorption energies of (I2, CH3I) vs. (CO, H2O), revealing possibilities for this material's use as an iodine sensor. The conclusions are supported by simulations at finite temperature; underlying electronic structures are also discussed.

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

Dynamic Catalytic Mechanism of the Methanol-to-Hydrocarbons Reaction over Zeolites

ConspectusThe methanol-to-hydrocarbons (MTH) process has provided a new route to obtaining basic chemicals without relying on an oil resource. Acidity and shape selectivity endow the zeolite with a decisive role in MTH catalysis. However, the inherent reaction characteristics of the MTH reaction over zeolites, such as the complexity of catalytic reaction kinetics, the diversity of catalytic reaction modes, and even the limitations of catalytic and diffusive decoupling, have all confused people with respect to obtaining a comprehensive mechanistic understanding. By examining the zeolite-catalyzed MTH reaction from the perspective of chemical bonding, one would realize that this reaction reflects the dynamic assembly process of C-C bonds from C1 components to multicarbon products. The key to understanding the MTH reaction lies in the mechanism by which C-C bonds are formed and rearranged in the confined microenvironment of the channel or cage structures of zeolite catalysts to achieve shape-selective production.The applications of advanced in situ spectroscopy as well as computational chemistry provide tremendous opportunities for capturing and identifying the details of the structure and properties of reactants, intermediates, and products in the confined reaction space of zeolite channels or cages, observing the real-time dynamic evolution of the catalytic surface, and modeling the elementary reaction steps at the molecular and atomic levels.In this Account, the dynamic catalytic mechanism of the zeolite-catalyzed MTH reaction will be outlined based on decades of continuous research and in-depth understanding. The combination of advanced in situ spectroscopy and theoretical methods allowed us to observe and simulate the formation, growth, and aging process on the working catalyst surface and thus map the dynamical evolution of active sites from a Brønsted acid site (BAS) to an organic-inorganic hybrid supramolecule (OIHS) in the MTH reaction. Moreover, the ever-evolving dynamic succession of the OIHS from surface methoxy species (SMS) to active ion-pair complexes (AIPC) to inert complexes (IC) guided the dynamic autocatalytic process from initiation to sustaining and then to termination, resulting in a complex interlaced hypercycle reaction network. The concept of dynamic catalysis will provide deep insight into the complex catalytic mechanisms as well as the structure-activity relationships in MTH chemistry. More importantly, we are now getting closer to the nature of zeolite catalysis beyond the traditional view of BAS catalysis.

Oil Inclusions Found in Skeleton Crystals of Quartz Indicated the Existence of Organic Matter Surrounding Ancient Growth Environments

In nature, minerals record various origins and information for geology and geobiochemistry. Here, we investigated the origin of organic matter and growth mechanism of quartz with oil inclusion revealing fluorescence under short ultraviolet (UV) light, obtained from the clay vein at Shimanto-cho, Kochi, Shikoku Island, Japan. Geological investigation indicated that the oil-quartz was formed in hydrothermal metamorphic veins found in the late Cretaceous interbedded sandstone and mudstone. The obtained oil-quartz crystals are mostly double-terminated. Micro-X-ray computed tomography (microCT) indicated that oil-quartz crystals have various veins originating as skeleton structures along the quartz crystal {111} and {1-11} faces. Spectroscopic and chromatographic studies indicated that aromatic ester and tetraterpene (lycopene) molecules, which revealed fluorescence, were detected. Large molecular weight sterol molecules, such as C40, were also detected in the vein of oil-quartz. This investigation indicated that organic inclusions in mineral crystals would form with ancient microorganism culture environments.

Loading-Driven Diffusion Pathway Selectivity in Zeolites with Continuum Intersecting Channels

The diffusion processes in zeolites are important for heterogeneous catalysis. Herein, we show that unique zeolites with "continuum intersecting channels" (e.g., BEC, POS, and SOV), in which two intersections are proximal, are greatly significant to the diffusion process with spontaneous switching of the diffusion pathway under varied loading. At low loading, the synergy of strong adsorption sites and molecular reorientation in intersections contribute to almost exclusive molecular diffusion in smaller channels. With an increase in molecular loading, the adsorbates are transported preferentially in larger channels mainly due to the lower diffusion barrier inside continuum intersection channels. This work demonstrates the ability to adjust the prior diffusion pathway by controlling the molecular loading, which may be beneficial for the separation of the product and byproduct in heterogeneous catalysis.

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.

Metal-free phosphorus and boron-doped graphitic carbon nitride/zeolite hetero-linked particles for highly efficient green hydrogen production in methanol

Herein, the development of phosphorus and boron-doped graphitic carbon nitride/zeolite (P- and B-doped g-C₃N₄-zeolite) catalyst under three-step heating conditions was performed. The first step is to prepare g-C₃N₄ synthesis from urea at 500 °C. In the second step, the production of a B-doped zeolite-g-C₃N₄ catalyst by calcination of g-C₃N₄ and zeolite was obtained at a ratio of 1:1 with boric acid at 500 °C. In the third step, the obtained B-doped zeolite- g-C₃N₄ catalyst consists of the preparation of B- and P-doped g-C₃N₄-zeolite catalyst as a result of the hydrothermal method with phosphoric acid. Characterization studies of the obtained catalysts were carried out with X-ray powder diffraction (XRD), scanning electron microscopy (SEM), energy dispersive X-ray spectroscopy (EDS), X-ray photoelectron spectroscopy (XPS), and Fourier-transform infrared spectroscopy (FTIR). These obtained catalysts were used as a metal-free catalyst in the production of hydrogen (H₂-P) by sodium borohydride in methanol(NaBH₄-MR) for the first time in the literature. The hydrogen production rate (HGR) value for P- and B-doped g-C₃N₄-zeolite catalysts was 6250 ml min^-1 g^-1.

Nuclear quantum effects on zeolite proton hopping kinetics explored with machine learning potentials and path integral molecular dynamics

Proton hopping is a key reactive process within zeolite catalysis. However, the accurate determination of its kinetics poses major challenges both for theoreticians and experimentalists. Nuclear quantum effects (NQEs) are known to influence the structure and dynamics of protons, but their rigorous inclusion through the path integral molecular dynamics (PIMD) formalism was so far beyond reach for zeolite catalyzed processes due to the excessive computational cost of evaluating all forces and energies at the Density Functional Theory (DFT) level. Herein, we overcome this limitation by training first a reactive machine learning potential (MLP) that can reproduce with high fidelity the DFT potential energy surface of proton hopping around the first Al coordination sphere in the H-CHA zeolite. The MLP offers an immense computational speedup, enabling us to derive accurate reaction kinetics beyond standard transition state theory for the proton hopping reaction. Overall, more than 0.6 μs of simulation time was needed, which is far beyond reach of any standard DFT approach. NQEs are found to significantly impact the proton hopping kinetics up to ~473 K. Moreover, PIMD simulations with deuterium can be performed without any additional training to compute kinetic isotope effects over a broad range of temperatures.

Active Sites in H-Mordenite Catalysts Probed by NMR and FTIR

Mordenites are widely used in catalysis and environmental protection. The catalytic properties of mordenite are largely determined by the composition of its crystal framework, i.e., the SiO2/Al2O3 molar ratio (MR), and the cationic form. In H-mordenites, the most important characteristic becomes the structure and distribution of acid sites, which depends on the number and distribution of Al tetrahedra in the framework. In the present work, the local structure of these centers in H-mordenite catalysts with a nominal MR varied from 9.9 to 19.8 was studied in detail using a combination of magic angle spinning nuclear magnetic resonance (MAS NMR) and Fourier transform infrared spectroscopy (FTIR). 27Al MAS NMR indicates the presence of extra-framework Al in most of the studied samples that results in a higher real MR of the zeolitic framework compared to the nominal value. Concentrations of Lewis and Brønsted acid sites, as well as of silanol groups were estimated by elemental analysis, NMR, and FTIR spectroscopy. The values of site concentrations obtained from band intensities of adsorbed CO and those of OH groups are compared with the amount of framework and extra-framework aluminum. The advantages and restrictions of different methods of active site characterization are discussed.

Adsorption of Carbamazepine in All-Silica Zeolites Studied with Density Functional Theory Calculations

The anticonvulsant drug carbamazepine (-) is an emerging contaminant of considerable concern due to its hazard potential and environmental persistence. Previous experimental studies proposed hydrophobic zeolites as promising adsorbents for the removal of carbamazepine from water, but only a few framework types were considered in those investigations. In the present work, electronic structure calculations based on dispersion-corrected density functional theory (DFT) were used to study the adsorption of CBZ in eleven all-silica zeolites having different pore sizes and connectivities of the pore system (AFI, ATS, BEA, CFI, DON, FAU, IFR, ISV, MOR, SFH, SSF framework types). It was found that some zeolites with one-dimensional channels formed by twelve-membered rings (IFR, AFI) exhibit the highest affinity towards CBZ. A "good fit" of CBZ into the zeolite pores that maximizes dispersion interactions was identified as the dominant factor determining the interaction strength. Further calculations addressed the role of temperature (for selected systems) and of guest-guest interactions between coadsorbed CBZ molecules. In addition to predicting zeolite frameworks of particular interest as materials for selective CBZ removal, the calculations presented here also contribute to the atomic-level understanding of the interaction of functional organic molecules with all-silica zeolites.

Tracking Sub-Nano-Scale Structural Evolution in Zeolite Synthesis by In Situ High-Energy X-ray Total Scattering Measurement with Pair Distribution Function Analysis

The crystallization mechanism of zeolites remains unclarified to date because of lack of effective techniques in characterizing the local structures of amorphous precursors under synthetic conditions. Herein, in situ high-energy X-ray total scattering measurement with pair distribution function analysis is performed throughout the hydrothermal synthesis of SSZ-13 zeolite to investigate the amorphous-to-crystalline transformation at the sub-nano level in real time. Ordered four-membered rings (4Rs) are dominantly formed during the induction period, prior to the significant increase in the number of symmetric six- and eight-membered rings (6Rs and 8Rs) in the crystal growth stage. These preformed ordered 4Rs contribute to the formation of d6r and cha composite building units containing 6Rs and 8Rs with the assistance of the organic structure-directing agent, leading to the construction of embryonic zeolite crystallites, which facilitate the crystal growth through a particle attachment pathway. This work enriches the toolbox for better understanding the crystallization pathway of zeolites.

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.

Experimental and theoretical study on the regioselective bis- or polyalkylation of 6-amino-2-mercapto-3H-pyrimidin-4-one using zeolite nano-gold catalyst and a quantum hybrid computational method

The synthetic utility of 6-amino-2-mercapto-3H-pyrimidin-4-one 3 as building blocks for new poly (pyrimidine) by alkylation using the bis(halo) compounds and zeolite nano-gold as a catalyst was investigated. Furthermore, the experimental findings by the theoretical Density functional theory (DFT) computations at the DFT/B3LYP level of theory, utilizing the 6-311++G (d,p) basis set in the gas phase, were used to investigate the distinct phases for Regio isomer 11a & 12a and 11b & 12b compounds was fair and of good quality. The stability of the 12a and 12b phases is higher than the other Regio isomer 11a and 11b phases, according to DFT modelling. By computing HOMO and LUMO pictures, the electronic parameters: dipole moment of these compounds in the ground state were theoretically investigated. Non-linear optical (NLO) characteristics and quantum chemical parameters were examined using frontier molecular orbital (FMO) analysis. Natural bond orbital analysis was used to characterize the charge transfer of the electron density in the investigated compounds (NBO). The molecular electrostatic potential surfaces (MEPS) plots have been generated, and absorption spectral analysis in different solvents has been theoretically and experimentally examined to better understand the reactivity spots. At the B3LYP/6-311G (d,p) level of theory, thermodynamic properties were also calculated. Finally, DFT calculations were used to connect the structure-activity relationship (SAR) with real antibacterial results for compounds 12a and 12b.

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.

Physicochemical Properties of Water-Based Copolymer and Zeolite Composite Sustained-Release Membrane Materials

A nitrogen fertilizer slow-release membrane was proposed using polyvinyl alcohol (PVA), polyvinylpyrrolidone (PVP), epoxy resin, and zeolite as raw materials. The effects of the water-based copolymer (PVA:PVP) solution ratio A (A1−A4) and zeolite amount B (B1−B4) on the water absorption rate (XS), water permeability (TS), fertilizer permeability (TF), tensile strength (KL), elongation at break (DSL), and viscosity (ND) of the membrane were explored using the swelling method, a self-made device, and a universal testing machine. The optimal combination of the water-based copolymer and zeolite amount was determined by the coefficient-of-variation method. The results show that the effects of the decrease in A on KL and the increase in B on KL and DSL are promoted first and then inhibited. DSL and ND showed a negative response to the A decrease, whereas XS, TS, and TF showed a positive response. The effect of increasing B on ND, TS, and TF showed a zigzag fluctuation. In the condition of A1−A3, XS showed a negative response to the B increase, whereas in the condition of A4, XS was promoted first and then inhibited. Adding PVP and zeolite caused the hydroxyl stretching vibration peak of PVA at 3300 cm−1 to widen; the former caused the vibration peak to move to low frequencies, and the latter caused it to move to high frequencies. The XRD pattern shows that the highest peak of zeolite is located at 2θ = 7.18° and the crystallization peak of the composite membrane increases with the rise in the proportion of zeolite. Adding PVP made the surface of the membrane smooth and flat, and adding a small amount of zeolite improved the mechanical properties of the membrane and exhibited good compatibility with water-based copolymers. In the evaluation model of the physicochemical properties of sustained-release membrane materials, the weight of all indicators was in the following order: TF > ND > TS > KL > XL > DSL. The optimal membrane material for comprehensive performance was determined to be A2B3.

Charge-separation driven mechanism via acylium ion intermediate migration during catalytic carbonylation in mordenite zeolite

By employing ab initio molecular dynamic simulations, solid-state NMR spectroscopy, and two-dimensional correlation analysis of rapid scan Fourier transform infrared spectroscopy data, a new pathway is proposed for the formation of methyl acetate (MA) via the acylium ion (i.e.,CH3 − C ≡ O+) in 12-membered ring (MR) channel of mordenite by an integrated reaction/diffusion kinetics model, and this route is kinetically and thermodynamically more favorable than the traditional viewpoint in 8MR channel. From perspective of the complete catalytic cycle, the separation of these two reaction zones, i.e., the C-C bond coupling in 8MR channel and MA formation in 12MR channel, effectively avoids aggregation of highly active acetyl species or ketene, thereby reducing undesired carbon deposit production. The synergistic effect of different channels appears to account for the high carbonylation activity in mordenite that has thus far not been fully explained, and this paradigm may rationalize the observed catalytic activity of other reactions.

Optimal machine learning feature selection for assessing the mechanical properties of a zeolite framework

In this study, 45 and 249 critical features were discovered among 896 zeolite descriptors generated by the matminer package for estimating the shear and bulk moduli of zeolites, respectively. A database containing the mechanical properties of 873 zeolite structures, calculated using density functional theory, was used to train the machine learning regression model. The results of using these critical features with the LightGBM algorithm were rigorously compared with those from other regressors as well as with different sets of features. The comparison results indicate that the surrogate model with critical features increases the prediction accuracy of the bulk and shear moduli of zeolites by 17.3% and 10.6%, respectively, and reduces the prediction uncertainty by one-third of that achieved using previously available features. The suggested features originating from several physical and chemical groups highlight the unveiled relationships between the features and mechanical properties of zeolites. The robustness of the constructed model with 356 features was confirmed by applying a set of different training-test set ratios. We believe that the suggested critical features of zeolites can help to understand the mechanical behavior of a half million unlabeled hypothetical zeolite structures and accelerate the discovery of novel zeolites with unprecedented mechanical properties.

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.