Interaction of water with zeolites: a review

Solvent effect in H-BEA catalyzed cyclohexanol dehydration reaction

The solvent effect on H-BEA catalyzed cyclohexanol dehydration was investigated in water, dioxane, and cyclohexanol. The dynamic evolution of the Brønsted acid site of zeolite and its interaction with reactant molecules in different solvents were explored with ab initio molecular dynamics simulations, providing reliable configuration sampling to obtain configurations at equilibrium. Solvent profoundly changes the adsorption as well as the dehydration reaction of cyclohexanol in H-BEA, where the reaction is determined to follow the E2 mechanism in water and dioxane but the E1 mechanism in cyclohexanol untill saturation uptake. Near saturation uptake, all three solvents significantly reduce the cyclohexanol dehydration rates in H-BEA. Cyclohexanol loading also dramatically affects the kinetics of the dehydration reaction, displaying an overall decreasing trend with a local minimum present at intermediate loading of 6 molecules per unit cell, which is a result of the entropic effect associated with greater freedom of motion of the transition state. Rigorous quantification of enthalpy and entropy contributions to cyclohexanol adsorption and activation shed light on the solvent effect of zeolite-catalyzed alcohol dehydration.

Solvent structure and dynamics over Brønsted acid MWW zeolite nanosheets

In the liquid phase of heterogeneous catalysis, solvent plays an important role and governs the kinetics and thermodynamics of a reaction. Although it is often difficult to quantify the role of the solvent, it becomes particularly challenging when a zeolite is used as the catalyst. This difficulty arises from the complex nature of the liquid/zeolite interface and the different solvation environments around catalytically active sites. Here, we use ab initio molecular dynamics simulations to probe the local solvation structure and dynamics of methanol and water over MWW zeolite nanosheets with varying Brønsted acidity. We find that the zeolite framework and the number and location of the acid sites in the zeolite influence the structure and dynamics of the solvent. In particular, methanol is more likely to be in the vicinity of the aluminum (Al3+) at the T4 site than at T1 due to easy accessibility. The methanol oxygen binds strongly to the Al at the T4 site, weakening the Al-O for the bridging acid site, which results in the formation of the silanol group, significantly reducing the acidity of the site. The behavior of methanol is in direct contrast to that of water, where protons can easily propagate from the zeolite to the solvent molecules regardless of the acid site location. Our work provides molecular-level insights into how solvent interacts with zeolite surfaces, leading to an improved understanding of the catalytic site in the MWW zeolite nanosheet.

Incorporating solvent effects in DFT: insights from cation exchange in faujasites

Zeolites are versatile materials renowned for their extra-framework cation exchange capabilities, with applications spanning diverse fields, including nuclear waste treatment. While detailed experimental characterization offers valuable insight, density functional theory (DFT) proves particularly adept at investigating ion exchange in zeolites, owing to its atomic and electronic resolution. However, the prevalent occurrence of zeolitic ion exchange in aqueous environments poses a challenge to conventional DFT modeling, traditionally conducted in a vacuum. This study seeks to enhance zeolite modeling by systematically evaluating predictive differences across varying degrees of aqueous solvent inclusion. Specifically focusing on monovalent cation exchange in Na-X zeolites, we explore diverse modeling approaches. These range from simple dehydrated systems (representing bare reference states in vacuum) to more sophisticated models that incorporate aqueous solvent effects through explicit water molecules and/or a dielectric medium. Through comparative analysis of DFT and semi-empirical DFT approaches, along with their validation against experimental results, our findings underscore the necessity to concurrently consider explicit and implicit solvent effects for accurate prediction of zeolitic ionic exchange.

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.

Promoting the circular economy: Valorization of a residue from industrial char to activated carbon with potential environmental applications as adsorbents

Pyrolysis of residues enriched with carbon, such as in agroforestry or industrial activities, has been postulated as an emerging technology to promote the production of biofuels, contributing to the circular economy and minimizing waste. However, during the pyrolysis processes a solid fraction residue is generated. This work aims to study the viability of these chars to develop porous carbonaceous materials that can be used for environmental applications. Diverse chars discharged by an industrial pyrolysis factory have been activated with KOH. Concretely, the char residues came from the pyrolysis of olive stone, pine, and acacia splinters, spent residues fuel, and cellulose artificial casings. The changes in the textural, structural, and composition characteristics after the activation process were studied by N2 adsorption-desorption isotherms, scanning electron microscopy, FTIR, elemental analysis, and XPS. A great porosity was developed, SBET within 776-1186 m2 g-1 and pore volume of 0.37-0.59 cm3 g-1 with 70-90% of micropores contribution. The activated chars were used for the adsorption of CO2, leading to CO2 maximum uptakes of 90-130 mg g-1. There was a good correlation between the CO2 uptake with microporosity and oxygenated surface groups of the activated chars. Moreover, their ability to adsorption of contaminants in aqueous solution was also evaluated. Concretely, there was studied the adsorption of aqueous heavy metals, i.e., Cd, Cu, Ni, Pb, and Zn, and organic pollutants of emerging concern such as caffeine, diclofenac, and acetaminophen.

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.

Water structures on acidic zeolites and their roles in catalysis

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

Hydrogen Solubility in Confined Water

Confined water has demonstrated distinct structural and dynamic properties compared to bulk water. Although many studies have explored the water structure within simple geometries using materials such as carbon and silica, studies on gas solubility in confined water and the underlying physics of water structure-solubility remain limited. Recent research has illuminated the concept of "oversolubility", wherein gases display increased solubility within liquids confined in small pores compared to their bulk form. This study focuses on zeolites, naturally abundant materials with versatile applications, to study the hydrogen solubility within confined water through careful experimentation. Our findings underscore the relationship between the pore dimension and gas solubility enhancement within confined water. Hydrogen solubility is closely associated with the rearrangement of water molecules within the porous framework of the zeolite. Our research shows that a 2 nm pore size results in the greatest increase in hydrogen solubility in the water trapped inside the zeolite framework. The double donor-double acceptor (DDAA) bonds play a critical role in hydrogen solubility. Our research provides fundamental insight into the role of the molecular bonding type on hydrogen solubility in water, paving the way for potential applications in hydrogen storage and utilization.

Water-Induced Micro-Hydrophobic Effect Regulates Benzene Methylation in Zeolite

Water is a ubiquitous component in heterogeneous catalysis over zeolites and can significantly influence the catalyst performance. However, the detailed mechanism insights into zeolite-catalyzed reactions under microscale aqueous environment remain elusive. Here, using multiple dimensional solid-state NMR experiments coupled with ultrahigh magic angle spinning technique and theoretical simulations, we establish a fundamental understanding of the role of water in benzene methylation over ZSM-5 zeolite under water vapor conditions. We show that water competes with benzene for the active sites of zeolite and facilitates the bimolecular reaction mechanism. The growth of water clusters induces a micro-hydrophobic effect in zeolite pores, which reorients benzene molecules and drives their interactions with surface methoxy species (SMS) on zeolite. We identify the formation and evolution of active SMS-Benzene complexes in a microscale aqueous environment and demonstrate that their accumulation in zeolite pores boosts benzene conversion and methylation.

Impact of Low-Temperature Water Exposure and Removal on Zeolite HY

Aqueous-phase postsynthetic modifications of the industrially important Y-type zeolite are commonly used to change overall acid site concentrations, introduce stabilizing rare-earth cations, impart bifunctional character through metal cation exchange, and tailor the distribution of Brønsted and Lewis acid sites. Zeolite Y is known to undergo framework degradation in the presence of both vapor- and liquid-phase water at temperatures exceeding 100 °C, and rare-earth exchanged and stabilized HY catalysts are commonly used for fluidized catalytic cracking due to their increased hydrothermal resilience. Here, using detailed spectroscopy, crystallography, and flow-reactor experiments, we reveal unexpected decreases in Brønsted acid site (BAS) density for zeolite HY following exposure even to room-temperature liquid water. These data indicate that aqueous-phase ion-exchange procedures commonly used to modify zeolite Y are impacted by the liquid water and its removal, even when fractional heating rates and inert conditions much less severe than standard practice are used for catalyst dehydration. X-ray diffraction, thermogravimetric, and spectroscopic analyses reveal that the majority of framework degradation occurs during the removal of a strongly bound water fraction in HY, which does not form when NH4Y is immersed in liquid water and which leads to reduced acidity in HY even when dehydration conditions much milder than those typically practiced are employed. Na+-exchanged HY prepared via room-temperature aqueous dissolution demonstrates that Brønsted acid sites are lost in excess of the theoretical maximum that is possible from sodium titration. The structural impact of low-temperature aqueous-phase ion-exchange methods complicates the interpretation of subsequent data and likely explains the wide variation in reported acid site concentrations and catalytic activity of HY zeolites with high-Al content.

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.

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.

Dynamic Pt-OH-·H2O-Ag species mediate coupled electron and proton transfer for catalytic hydride reduction of 4-nitrophenol at the confined nanoscale interface

Generally, the catalytic transformation of 4-nitrophenol (4-NP) to 4-aminophenol (4-AP) at heterogeneous metal surfaces follows a Langmuir-Hinshelwood (L-H) mechanism when sodium borohydride (NaBH₄) is used as the sacrificial reductant. Herein, with Pt-Ag bimetallic nanoparticles confined in dendritic mesoporous silica nanospheres (DMSNs) as a model catalyst, we demonstrated that the conversion of 4-NP did not pass through the direct hydrogen transfer route with the hydride equivalents being supplied by borohydride via the bimolecular L-H mechanism, since Fourier transform infrared (FTIR) spectroscopy with the use of isotopically labeled reactants (NaBD₄ and D₂O) showed that the final product of 4-AP was composed of protons (or deuterons) that originated from the solvent water (or heavy water). Combined characterization by X-ray photoelectron spectroscopy (XPS), ¹H nuclear magnetic resonance (NMR) and the optical excitation and photoluminescence spectrum evidenced that the surface hydrous hydroxide complex bound to the metal surface (also called structural water molecules, SWs), due to the space overlap of p orbitals of two O atoms in SWs, could form an ensemble of dynamic interface transient states, which provided the alternative electron and proton transfer channels for selective transformation of 4-NP. The cationic Pt species in the Ag-Pt bimetallic catalyst mainly acts as a dynamic adsorption center to temporally anchor SWs and related reactants, and not as the active site for hydrogen activation.