Structure and solvation of confined water and water-ethanol clusters within microporous Brønsted acids and their effects on ethanol dehydration catalysis

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.

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.

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.

New Aspects of Alcohol-Alcohol and Alcohol-Water Interactions: Crystallographic and Quantum Chemical Studies of Antiparallel O-H/O-H Interactions

New modes of interaction, antiparallel O-H/O-H interactions of alcohol-alcohol dimers and alcohol-water dimers, were studied by analyzing data in the Cambridge Structural Database (CSD) and by calculating potential energy surfaces at a very accurate quantum chemical CCSD(T)/CBS level. The data reveal the existence of antiparallel interactions in crystal structures and significant interaction energies. Data from the CSD for alcohol-alcohol dimers show 49.2% of contacts with classical hydrogen bonds and 10.1% of contacts with antiparallel interactions, while for alcohol-water dimers, 59.4% of contacts are classical hydrogen bonds and only 0.6% of contacts are antiparallel interactions. The calculations were performed on methanol, ethanol, and n-propanol dimers. Classical hydrogen-bonded alcohol-alcohol and alcohol-water dimers have interaction energies of up to -6.2 kcal/mol and up to -5.5 kcal/mol, respectively. Antiparallel interactions in alcohol-alcohol and alcohol-water dimers have interaction energies of up to -4.7 kcal/mol and up to -4.4 kcal/mol, respectively. Symmetry-adapted perturbation theory analysis for antiparallel interactions shows their electrostatic nature.

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.

Heterogeneous Ice Nucleation in Model Crystalline Porous Organic Polymers: Influence of Pore Size on Immersion Freezing

Heterogeneous ice nucleation activity is affected by aerosol particle composition, crystallinity, pore size, and surface area. However, these surface properties are not well understood, regarding how they act to promote ice nucleation and growth to form ice clouds. Therefore, synthesized materials for which surface properties can be tuned were examined in immersion freezing mode in this study. To establish the relationship between particle surface properties and efficiency of ice nucleation, materials, here, covalent organic frameworks (COFs), with different pore diameters and degrees of crystallinity (ordering), were characterized. Results showed that out of all the highly crystalline COFs, the sample with a pore diameter between 2 and 3 nm exhibited the most efficient ice nucleation activity. We posit that the highly crystalline structures with ordered pores have an optimal pore diameter where the ice nucleation activity is maximized and that the not highly crystalline structures with nonordered pores have more sites for ice nucleation. The results were compared and discussed in the context of other synthesized porous particle systems. Such studies give insight into how material features impact ice nucleation activity.

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.

Concepts Relevant for the Kinetic Analysis of Reversible Reaction Systems

The net rate of a reversible chemical reaction is the difference between unidirectional rates of traversal along forward and reverse reaction paths. In a multistep reaction sequence, the forward and reverse trajectories, in general, are not the microscopic reverse of one another; rather, each unidirectional route is comprised of distinct rate-controlling steps, intermediates, and transition states. Consequently, traditional descriptors of rate (e.g., reaction orders) do not reflect intrinsic kinetic information but instead conflate unidirectional contributions determined by (i) the microscopic occurrence of forward/reverse reactions (i.e., unidirectional kinetics) and (ii) the reversibility of reaction (i.e., nonequilibrium thermodynamics). This review aims to provide a comprehensive resource of analytical and conceptual tools which deconvolute the contributions of reaction kinetics and thermodynamics to disambiguate unidirectional reaction trajectories and precisely identify rate- and reversibility-controlling molecular species and steps in reversible reaction systems. The extrication of mechanistic and kinetic information from bidirectional reactions is accomplished through equation-based formalisms (e.g., De Donder relations) grounded in principles of thermodynamics and interpreted in the context of theories of chemical kinetics developed in the past 25 years. The aggregate of mathematical formalisms detailed herein is general to thermochemical and electrochemical reactions and encapsulates a diverse body of scientific literature encompassing chemical physics, thermodynamics, chemical kinetics, catalysis, and kinetic modeling.

Conceptual and Practical Aspects of Metal-Organic Frameworks for Solid-Gas Reactions

The presence of site-isolated and well-defined metal sites has enabled the use of metal-organic frameworks (MOFs) as catalysts that can be rationally modulated. Because MOFs can be addressed and manipulated through molecular synthetic pathways, they are chemically similar to molecular catalysts. They are, nevertheless, solid-state materials and therefore can be thought of as privileged solid molecular catalysts that excel in applications involving gas-phase reactions. This contrasts with homogeneous catalysts, which are overwhelmingly used in the solution phase. Herein, we review theories dictating gas phase reactivity within porous solids and discuss key catalytic gas-solid reactions. We further treat theoretical aspects of diffusion within confined pores, the enrichment of adsorbates, the types of solvation spheres that a MOF might impart on adsorbates, definitions of acidity/basicity in the absence of solvent, the stabilization of reactive intermediates, and the generation and characterization of defect sites. The key catalytic reactions we discuss broadly include reductive reactions (olefin hydrogenation, semihydrogenation, and selective catalytic reduction), oxidative reactions (oxygenation of hydrocarbons, oxidative dehydrogenation, and carbon monoxide oxidation), and C-C bond forming reactions (olefin dimerization/polymerization, isomerization, and carbonylation reactions).

Boosting Catalytic Selectivity through a Precise Spatial Control of Catalysts at Pickering Droplet Interfaces

Exploration of new methodologies to tune catalytic selectivity is a long-sought goal in catalytic community. In this work, oil-water interfaces of Pickering emulsions are developed to effectively regulate catalytic selectivity of hydrogenation reactions, which was achieved via a precise control of the spatial distribution of metal nanoparticles at the droplet interfaces. It was found that Pd nanoparticles located in the inner interfacial layer of Pickering droplets exhibited a significantly enhanced selectivity for p-chloroaniline (up to 99.6%) in the hydrogenation of p-chloronitrobenzene in comparison to those in the outer interfacial layer (63.6%) in pure water (68.5%) or in pure organic solvents (46.8%). Experimental and theoretical investigations indicated that such a remarkable interfacial microregion-dependent catalytic selectivity was attributed to the microenvironments of the coexistence of water and organic solvent at the droplet interfaces, which could provide unique interfacial hydrogen-bonding interactions and solvation effects so as to alter the adsorption patterns of p-chloronitrobenzene and p-chloroaniline on the Pd nanoparticles, thereby avoiding the unwanted contact of C-Cl bonds with the metal surfaces. Our strategy of precise spatial control of catalysts at liquid-liquid interfaces and the unprecedented interfacial effect reported here not only provide new insights into the liquid-liquid interfacial reactions but also open an avenue to boost catalytic selectivity.

Maximum Impact of Ionic Strength on Acid-Catalyzed Reaction Rates Induced by a Zeolite Microporous Environment

The intracrystalline ionic environment in microporous zeolite can remarkably modify the excess chemical potential of adsorbed reactants and transition states, thereby influencing the catalytic turnover rates. However, a limit of the rate enhancement for aqueous-phase dehydration of alcohols appears to exist for zeolites with high ionic strength. The origin of such limitation has been hypothesized to be caused by the spatial constraints in the pores via, e.g., size exclusion effects. It is demonstrated here that the increase in turnover rate as well as the formation of a maximum and the rate drop are intrinsic consequences of the increasingly dense ionic environment in zeolite. The molecularly sized confines of zeolite create a unique ionic environment that monotonically favors the formation of alcohol-hydronium ion complexes in the micropores. The zeolite microporous environment determines the kinetics of catalytic steps and tailors the impact of ionic strength on catalytic rates.

Influence of Substrate Concentration on Kinetic Parameters of Ethanol Dehydration in MFI and CHA Zeolites and Relation of These Kinetic Parameters to Acid–Base Properties

The catalytic activity of zeolites is often related to their acid–base properties. In this work, the relationship between the value of apparent activation energy of ethanol dehydration, measured in a fixed bed reactor and by means of a temperature-programmed surface reaction (TPSR) depending on the amount of ethanol in the zeolite lattice and the value of activation energy of H/D exchange as a measure of acid–base properties of MFI and CHA zeolites, was studied. Tests in a fixed bed reactor were unable to provide reliable reaction kinetics data due to internal diffusion limitations and rapid catalyst deactivation. Only the TPSR method was able to provide activation energy values comparable to the activation energy values obtained from the H/D exchange rate measurements. In addition, for CHA zeolite, it has been shown that the values of ethanol dehydration activation energies depend on the amount of ethanol in the CHA framework, and this effect can be attributed to the substrate clustering effects supporting the deprotonation of zeolite Brønsted centers.

Ab initio molecular dynamics with enhanced sampling in heterogeneous catalysis

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

Dynamics of Heterogeneous Catalytic Processes at Operando Conditions

The rational design of high-performance catalysts is hindered by the lack of knowledge of the structures of active sites and the reaction pathways under reaction conditions, which can be ideally addressed by an in situ/operando characterization. Besides the experimental insights, a theoretical investigation that simulates reaction conditions-so-called operando modeling-is necessary for a plausible understanding of a working catalyst system at the atomic scale. However, there is still a huge gap between the current widely used computational model and the concept of operando modeling, which should be achieved through multiscale computational modeling. This Perspective describes various modeling approaches and machine learning techniques that step toward operando modeling, followed by selected experimental examples that present an operando understanding in the thermo- and electrocatalytic processes. At last, the remaining challenges in this area are outlined.

Quantifying Effects of Active Site Proximity on Rates of Methanol Dehydration to Dimethyl Ether over Chabazite Zeolites through Microkinetic Modeling

Control of the spatial proximity of Brønsted acid sites within the zeolite framework can result in materials with properties that are distinct from materials synthesized through conventional crystallization methods or available from commercial sources. Recent experimental evidence has shown that turnover rates of different acid-catalyzed reactions increase with the fraction of proximal sites in chabazite (CHA) zeolites. The catalytic conversion of oxygenates is an important research area, and the dehydration of methanol to dimethyl ether (DME) is a well-studied reaction as part of methanol-to-olefin chemistry catalyzed by solid acids. Published experimental data have shown that DME formation rates (per acid site) increase systematically with the fraction of proximal acid sites in the six-membered ring of CHA. Here, we probe the effect of acid site proximity in CHA on methanol dehydration rates using electronic structure calculations and microkinetic modeling to identify the primary causes of this chemistry and their relationship to the local structure of the catalyst at the nanoscale. We report a density functional theory-parametrized microkinetic model of methanol dehydration to DME, catalyzed by acidic CHA zeolite with direct comparison to experimental data. Effects of proximal acid sites on reaction rates were captured quantitatively for a range of operating conditions and catalyst compositions, with a focus on total paired acid site concentration and reactant clustering to form higher nuclearity complexes. Next-nearest neighbor paired acid sites were identified as promoting the formation of methanol trimer clusters rather than the inhibiting tetramer or pentamer clusters, resulting in large increases in the rate for DME production due to the lower energy barriers present in the concerted methanol trimer reaction pathway. The model framework developed in this study can be extended to other zeolite materials and reaction chemistries toward the goal of rational design and development of next-generation catalytic materials and chemical processes.

Impact of the Local Concentration of Hydronium Ions at Tungstate Surfaces for Acid-Catalyzed Alcohol Dehydration

Tungstate domains supported on ZrO₂, Al₂O₃, TiO₂, and activated carbon drastically influence the hydronium-ion-catalyzed aqueous-phase dehydration of alcohols. For all catalysts, the rate of cyclohexanol dehydration normalized to the concentration of Brønsted acid sites (turnover frequencies, TOFs) was lower for monotungstates than for polytungstates and larger crystallites of WO₃. TOFs were constant when reaching or exceeding the monolayer coverage of tungstate, irrespective of the specific nature of surface structures that continuously evolve with the surface W loading. However, the TOFs with polytungstates and large WO₃ crystallites depend strongly on the underlying support (e.g., WO_x/C catalysts are 10-50-fold more active than WO_x/Al₂O₃ catalysts). The electrical double layer (EDL) surrounding the negatively charged WO_x domains contains hydrated hydronium ions, whose local concentrations change with the support. This varying concentration of interfacial hydronium ions ("local ionic strength") impacts the excess chemical potential of the reacting alcohols and induces the marked differences in the TOFs. Primary H/D kinetic isotope effects (∼3), together with the substantially positive entropy of activation (111-195 J mol^-1 K^-1), indicate that C-H(D) bond cleavage is involved in the kinetically relevant step of an E1-type mechanistic sequence, regardless of the support identity. The remarkable support dependence of the catalytic activity observed here for the aqueous-phase dehydration of cycloalkanols likely applies to a broad set of hydronium-ion-catalyzed organic reactions sensitive to ionic strength.

Influence of solvent structure and hydrogen bonding on catalysis at solid-liquid interfaces

Solvent molecules interact with reactive species and alter the rates and selectivities of catalytic reactions by orders of magnitude. Specifically, solvent molecules can modify the free energies of liquid phase and surface species via solvation, participating directly as a reactant or co-catalyst, or competitively binding to active sites. These effects carry consequences for reactions relevant for the conversion of renewable or recyclable feedstocks, the development of distributed chemical manufacturing, and the utilization of renewable energy to drive chemical reactions. First, we describe the quantitative impact of these effects on steady-state catalytic turnover rates through a rate expression derived for a generic catalytic reaction (A → B), which illustrates the functional dependence of rates on each category of solvent interaction. Second, we connect these concepts to recent investigations of the effects of solvents on catalysis to show how interactions between solvent and reactant molecules at solid-liquid interfaces influence catalytic reactions. This discussion demonstrates that the design of effective liquid phase catalytic processes benefits from a clear understanding of these intermolecular interactions and their implications for rates and selectivities.

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

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

Insights into protonation for cyclohexanol/water mixtures at the zeolitic Brønsted acid site

Proton transfer from Brønsted acid sites (BASs) to alcohol molecules ignites the acid-catalyzed alcohol dehydration reactions. For aqueous phase dehydration reactions in zeolites, the coexisting water molecules around BASs in the zeolite pores significantly affect the alcohol dehydration activity. In the present work, proton transfer processes among the BASs of H-BEA zeolites, the adsorbed cyclohexanol and surrounding water clusters with different sizes up to 8 water molecules were investigated using ab initio molecular dynamics (AIMD) simulations combined with the multiple-walker well-tempered metadynamics algorithm. The plausible proton locations and proton transfer processes were characterized using two/three-dimensional free energy landscapes. The strong proton affinity makes the protonated cyclohexanol stable species until a water trimer is formed. The proton either is shared between protonated cyclohexanol and the water trimer or remains with the water trimer (H7O3+). With a further increase in water concentrations, the proton prefers to remain with the water clusters.

Kinetics of H2 Adsorption at the Metal-Support Interface of Au/TiO2 Catalysts Probed by Broad Background IR Absorbance

H₂ adsorption on Au catalysts is weak and reversible, making it difficult to quantitatively study. We demonstrate H₂ adsorption on Au/TiO₂ catalysts results in electron transfer to the support, inducing shifts in the FTIR background. This broad background absorbance (BBA) signal is used to quantify H₂ adsorption; adsorption equilibrium constants are comparable to volumetric adsorption measurements. H₂ adsorption kinetics measured with the BBA show a lower E_app value (23 kJ mol^-1 ) for H₂ adsorption than previously reported from proxy H/D exchange (33 kJ mol^-1 ). We also identify a previously unreported H-O-H bending vibration associated with proton adsorption on electronically distinct Ti-OH metal-support interface sites, providing new insight into the nature and dynamics of H₂ adsorption at the Au/TiO₂ interface.

Kinetic effects of molecular clustering and solvation by extended networks in zeolite acid catalysis

Reactions catalyzed within porous inorganic and organic materials and at electrochemical interfaces commonly occur at high coverage and in condensed media, causing turnover rates to depend strongly on interfacial structure and composition, collectively referred to as "solvent effects". Transition state theory treatments define how solvation phenomena enter kinetic rate expressions, and identify two distinct types of solvent effects that originate from molecular clustering and from the solvation of such clusters by extended solvent networks. We review examples from the recent literature that investigate reactions within microporous zeolite catalysts to illustrate these concepts, and provide a critical appraisal of open questions in the field where future research can aid in developing new chemistry and catalyst design principles.

Ordered Hydrogen-Bonded Alcohol Networks Confined in Lewis Acid Zeolites Accelerate Transfer Hydrogenation Turnover Rates

The disruption of ordered water molecules confined within hydrophobic reaction pockets alters the energetics of adsorption and catalysis, but a mechanistic understanding of how nonaqueous solvents influence catalysis in microporous voids remains unclear. Here, we use kinetic analyses coupled with IR spectroscopy to study how alkanol hydrogen-bonding networks confined within hydrophobic and hydrophilic zeolite catalysts modify reaction free energy landscapes. Hydrophobic Beta zeolites containing framework Sn atoms catalyze the transfer hydrogenation reaction of cyclohexanone in a 2-butanol solvent 10× faster than their hydrophilic analogues. This rate enhancement stems from the ability of hydrophobic Sn-Beta to inhibit the formation of extended liquid-like 2-butanol oligomers and promote dimeric H-bonded 2-butanol networks. These different intraporous 2-butanol solvent structures manifest as differences in the activation and adsorption enthalpies and entropies that comprise the free energy landscape of transfer hydrogenation catalysis. The ordered H-bonding solvent network present in hydrophobic Sn-Beta stabilizes the transfer hydrogenation transition state to a greater extent than the liquid-like 2-butanol solvent present in hydrophilic Sn-Beta, giving rise to higher turnover rates on hydrophobic Sn-Beta. Additionally, reactant adsorption within hydrophobic Sn-Beta is driven by the breakup of intraporous solvent-solvent interactions, resulting in positive enthalpies of adsorption that are partially compensated by an increase in the solvent reorganization entropy. Collectively, these results emphasize the ability of the zeolite pore to regulate the structure of confined nonaqueous H-bonding solvent networks, which offers an additional dimension to modulate adsorption and reactivity.

Rigid Arrangements of Ionic Charge in Zeolite Frameworks Conferred by Specific Aluminum Distributions Preferentially Stabilize Alkanol Dehydration Transition States

Zeolite reactivity depends on the solvating environments of their micropores and the proximity of their Brønsted acid sites. Turnover rates (per H⁺ ) for methanol and ethanol dehydration increase with the fraction of H⁺ sites sharing six-membered rings of chabazite (CHA) zeolites. Density functional theory (DFT) shows that activation barriers vary widely with the number and arrangement of Al (1-5 per 36 T-site unit cell), but cannot be described solely by Al-Al distance or density. Certain Al distributions yield rigid arrangements of anionic charge that stabilize cationic intermediates and transition states via H-bonding to decrease barriers. This is a key feature of acid catalysis in zeolite solvents, which lack the isotropy of liquid solvents. The sensitivity of polar transition states to specific arrangements of charge in their solvating environments and the ability to position such charges in zeolite lattices with increasing precision herald rich catalytic diversity among zeolites of varying Al arrangement.