Impact of long-range electrostatic and dispersive interactions on theoretical predictions of adsorption and catalysis in zeolites

Level-Shifted Embedded Cluster Method for Modeling the Chemistry of Metal Oxides

The embedded cluster method has been used extensively in the study of the chemical and physical properties of metal oxides. This method has been a popular tool due to its relatively high accuracy and low computational cost. An even more promising option may entail integrating the embedded cluster method with the combined quantum mechanical and molecular mechanical (QM/MM) approach, thereby enabling further consideration of interactions within the entire system for superior results. We aim to accurately model the chemistry of metal oxides using this combined scheme. Here, using the prototypical MgO(100) surface as a test system, with Mg9O14 as the cluster in the quantum mechanical region, we show that the embedded cluster with untailored boundary effective core potentials (ECPs) can have frontier orbital energy levels that substantially deviate from the quantum mechanical reference results. This occurs even when Mg9O9, which retains the stoichiometry of MgO, is used as the cluster in the quantum mechanical region. As a result, the chemical properties of the embedded cluster models differ from those of the quantum mechanical reference model. To address this issue, we propose a new variant of the embedded cluster method called the level-shifted embedded cluster (LSEC) method, which allows the energy levels to be shifted to match the reference levels by tuning the boundary ECPs. Our validation calculations on the adsorption of various adsorbates with different properties on the MgO(100) surface show that the overall performance of QM/MM with the LSEC method is excellent for the adsorption energies, geometries, and charge properties. The excellent performance holds for both the nonstoichiometric and stoichiometric clusters (i.e., Mg9O14 and Mg9O9, respectively), demonstrating the robustness of the LSEC method. We expect that the LSEC method can be combined with QM/MM or used separately for future chemical studies of metal oxides and other ionically bonded systems.

Zeolite encapsulated organometallic complexes as model catalysts

Heterogeneities in the structure of active centers in metal-containing porous materials are unavoidable and complicate the description of chemical events occurring along reaction coordinates at the atomic level. Metal containing zeolites include sites of varied local coordination and secondary confining environments, requiring careful titration protocols to quantify the predominant active sites. Hybrid organometallic-zeolite catalysts are useful well-defined platform materials for spectroscopic, kinetic, and computational studies of heterogeneous catalysis that avoid the complications of conventional metal-containing porous materials. Such materials have been synthesized and studied previously, but catalytic applications were mostly limited to liquid-phase oxidation and electrochemical reactions. The hydrothermal stability, time-on-stream stability, and utility of these materials in gas-phase oxidation reactions are under-studied. The potential applications for single-site heterogeneous catalysts in fundamental research are abundant and motivate future synthetic, spectroscopic, kinetic, and computational studies.

Effect of the QM Size, Basis Set, and Polarization on QM/MM Interaction Energy Decomposition Analysis

Herein, an Energy Decomposition Analysis (EDA) scheme extended to the framework of QM/MM calculations in the context of electrostatic embeddings (QM/MM-EDA) including atomic charges and dipoles is applied to assess the effect of the QM region size on the convergence of the different interaction energy components, namely, electrostatic, Pauli, and polarization, for cationic, anionic, and neutral systems interacting with a strong polar environment (water). Significant improvements are found when the bulk solvent environment is described by a MM potential in the EDA scheme as compared to pure QM calculations that neglect bulk solvation. The predominant electrostatic interaction requires sizable QM regions. The results reported here show that it is necessary to include a surprisingly large number of water molecules in the QM region to obtain converged values for this energy term, contrary to most cluster models often employed in the literature. Both the improvement of the QM wave function by means of a larger basis set and the introduction of polarization into the MM region through a polarizable force field do not translate to a faster convergence with the QM region size, but they lead to better results for the different interaction energy components. The results obtained in this work provide insight into the effect of each energy component on the convergence of the solute-solvent interaction energy with the QM region size. This information can be used to improve the MM FFs and embedding schemes employed in QM/MM calculations of solvated systems.

Carbocation chemistry confined in zeolites: spectroscopic and theoretical characterizations

Acid-catalyzed reactions inside zeolites are one type of broadly applied industrial reactions, where carbocations are the most common intermediates of these reaction processes, including methanol to olefins, alkene/aromatic alkylation, and hydrocarbon cracking/isomerization. The fundamental research on these acid-catalyzed reactions is focused on the stability, evolution, and lifetime of carbocations under the zeolite confinement effect, which greatly affects the efficiency, selectivity and deactivation of zeolite catalysts. Therefore, a profound understanding of the carbocations confined in zeolites is not only beneficial to explain the reaction mechanism but also drive the design of new zeolite catalysts with ideal acidity and cages/channels. In this review, we provide both an in-depth understanding of the stabilization of carbocations by the pore confinement effect and summary of the advanced characterization methods to capture carbocations in zeolites, including UV-vis spectroscopy, solid-state NMR, fluorescence microscopy, IR spectroscopy and Raman spectroscopy. Also, we clarify the relationship between the activity and stability of carbocations in zeolite-catalyzed reactions, and further highlight the role of carbocations in various hydrocarbon conversion reactions inside zeolites with diverse frameworks and varying acidic properties.

Assessing the stability of Pd-exchanged sites in zeolites with the aid of a high throughput quantum chemistry workflow

Cation exchanged-zeolites are functional materials with a wide range of applications from catalysis to sorbents. They present a challenge for computational studies using density functional theory due to the numerous possible active sites. From Al configuration, to placement of extra framework cation(s), to potentially different oxidation states of the cation, accounting for all these possibilities is not trivial. To make the number of calculations more tractable, most studies focus on a few active sites. We attempt to go beyond these limitations by implementing a workflow for a high throughput screening, designed to systematize the problem and exhaustively search for feasible active sites. We use Pd-exchanged CHA and BEA to illustrate the approach. After conducting thousands of explicit DFT calculations, we identify the sites most favorable for the Pd cation and discuss the results in detail. The high throughput screening identifies many energetically favorable sites that are non-trivial. Lastly, we employ these results to examine NO adsorption in Pd-exchanged CHA, which is a promising passive NOx adsorbent (PNA) during the cold start of automobiles. The results shed light on critical active sites for NOx capture that were not previously studied.

Invertase adsorption with polymers functionalized by aspartic acid

Today, the separation and purification processes are highly preferred over the affinity interactions in the scientific world. Among the materials used for this purpose, magnetic particles and cryogels are very popular. Both polymeric structures have their advantages and disadvantages. In this study, poly(2-Hydroxyethyl methacrylate-N-methacryloyl-L-aspartic acid), poly(HEMA-MAsp), magnetic microparticles, and cryogels were synthesized, and adsorption performances of both polymeric structures were investigated by using invertase from aqueous systems. Invertase (β-fructofuranoside fructohydrolase, EC 3.2.1.26) is a commercially important enzyme used in the food industry to obtain the product called invert sugar, which consists of a mixture of equivalent amounts of glucose and fructose. Therefore, it was preferred as a model enzyme in adsorption studies of polymeric structures. According to the results, 104.1 mg g−1 and 135.5 mg g−1 of adsorption capacity values were obtained for cryogel and magnetic microparticle forms, respectively. Increasing temperature slightly reduced the adsorption capacity of both polymeric structures. In the adsorption/desorption cycle studies performed five times with poly(HEMA-MAsp) polymers, both forms were found to have high reusable properties. It was determined that the activity of invertase immobilized on polymeric structures was preserved at a rate of 83.6% for the particle form and 89.2% for the cryogel form.

Monomolecular Cracking of Propane: Effect of Zeolite Confinement and Acidity

The effect of zeolite pore geometry and intrinsic acidity on the activation energy of propane monomolecular cracking was investigated for six topologically distinct zeolites with different pore sizes. Periodic density functional theory calculations were used to calculate the activation energy, while cluster models were used to calculate deprotonation energies. The computed intrinsic activation energies showed a smaller variation with topology than the adsorption energies. No correlation was found between the computed deprotonation and ammonia adsorption energies at the acid site and the intrinsic activation energy. Detailed analysis of the computed structures and properties suggests that acid sites with different pore topologies impose geometrical constraints on the ion-pair formed by the ammonium molecule, which differs significantly from those that affect the propane reaction.

Second-Sphere Lattice Effects in Copper and Iron Zeolite Catalysis

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

Computational Modeling Predicts the Stability of Both Pd+ and Pd2+ Ion-Exchanged into H-CHA

Passive NOx adsorbers (PNA) using Pd/zeolites have emerged as a promising solution for the reduction of cold-start emissions from vehicle exhaust. However, the nature of the active sites and the mechanisms underlying NOx adsorption in Pd/zeolites remain a subject of ongoing investigation. In this study, we employ quantum chemical simulations to investigate the structure of Pd species in cation-exchange sites at isolated Al and Al pairs in the 6-ring and 8-ring of the CHA framework, before the introduction of NOx. Our calculations show that the speciation of Pd in these exchange sites strongly depends on the precise Al arrangement within the framework, as well as the operating conditions. Ionically dispersed Pd is found to be the most favorable species over a wide range of oxidizing and reducing conditions. Small oligomers of PdO and metallic Pd do not appear to be competitive at either isolated Al or Al pairs. Notably, our calculations show that ion exchange sites other than next-next-nearest neighbor Al pairs in the 6-ring will be preferentially occupied by Pd+ instead of Pd2+. The stability of Pd+ in the zeolite environment is an interesting contrast with its rareness in molecular Pd compounds. Nonetheless, a detailed analysis of the electronic structure shows that predicted Pd oxidation states are consistent with chemical intuition for all complexes investigated in this study. We also discuss the potential ambiguity in Pd characterization provided by typical experimental techniques such as XANES, EXAFS and UV-VIS, and highlight the need for additional EPR spectroscopy studies to further elucidate the initial Pd speciation in zeolites for PNA applications.

Sorption of Pb(II) from Aqueous Solutions by Acid-Modified Clinoptilolite-Rich Tuffs with Different Si/Al Ratios

The removal of Pb(II) from aqueous solutions by acid-modified clinoptilolite-rich tuff was investigated in this work. Clinoptilolite-rich tuff samples were treated using H2SO4 at different concentrations. Prior to and following acid treatment, the samples were characterized using X-ray diffraction (XRD), scanning electron microscopy (SEM), energy dispersive X-ray spectroscopy (EDS), and Fourier-transform infrared spectroscopy (FTIR). The pH of the point of zero charge (pHPZC) was also determined as part of this characterization. Batch studies were studied to investigate Pb(II) removal as a function of contact time, initial Pb(II) concentration, adsorbent dosage, and solution pH. The results of the XRD and SEM techniques showed that clinoptilolite is the main mineral of the non- and acid-treated natural zeolite samples. However, EDS analysis indicated that the Si/Al ratio increases as the exchangeable ions decrease with increasing acid concentrations. The optimum conditions for Pb(II) removal for samples with 4.37 ≤ Si/Al ≤ 7.9 were found to be as follows: Contact time of 60–360 min, pH: 6–8, and adsorbent dose of 6 mg g−1; whereas for acid-modified clinoptilolite-rich tuffs with 9.01 ≤ Si/Al ≤ 9.52, these conditions were as follows: Contact time of 1440 min, pH: 8–10, and adsorbent dose of 10 mg g−1. The experimental data were analyzed by kinetic and isotherms models. The results showed that the sorption of Pb(II) on samples with Si/Al ratios of 4.37, 5.31, and 7.91 were in agreement with the pseudo-second order and Langmuir isotherm with qm = 48.54, 37.04, and 14.99 mg g−1, respectively, while the kinetic data and isotherm for samples with 9.01 ≤ Si/Al ≤ 9.52 were found to fit the pseudo-first order and Freundlich model.

Mechanism and Kinetics of Propane Dehydrogenation and Cracking over Ga/H-MFI Prepared via Vapor-Phase Exchange of H-MFI with GaCl3

In this study, the mechanism and kinetics of C₃H₈ dehydrogenation and cracking are examined over Ga/H-MFI catalysts prepared via vapor-phase exchange of H-MFI with GaCl₃. The present study demonstrates that [GaH]²⁺ cations are the active centers for C₃H₈ dehydrogenation and cracking, independent of the Ga/Al ratio. For identical reaction conditions, [GaH]²⁺ cations in Ga/H-MFI exhibit a turnover frequency for C₃H₈ dehydrogenation that is 2 orders of magnitude higher and for C₃H₈ cracking, that is 1 order of magnitude higher than the corresponding turnover frequencies over H-MFI. C₃H₈ dehydrogenation and cracking exhibit first-order kinetics with respect to C₃H₈ over H-MFI, but both reactions exhibit first-order kinetics over Ga/H-MFI only at very low C₃H₈ partial pressures and zero-order kinetics at higher C₃H₈ partial pressures. H₂ inhibits both reactions over Ga/H-MFI. It is also found that the ratio of the rate of dehydrogenation to the rate of cracking over Ga/H-MFI is independent of C₃H₈ and H₂ partial pressures but weakly dependent on temperature. Measured activation enthalpies together with theoretical analysis are consistent with a mechanism in which both the dehydrogenation and cracking of C₃H₈ proceed over Ga/H-MFI via reversible, heterolytic dissociation of C₃H₈ at [GaH]²⁺ sites to form [C₃H₇-GaH]⁺-H⁺ cation pairs. The rate-determining step for dehydrogenation is the β-hydride elimination of C₃H₆ and H₂ from the C₃H₇ fragment. The rate-determining step for cracking is C-C bond attack of the same propyl fragment by the proximal Brønsted acid O-H group. H₂ inhibits both dehydrogenation and cracking over Ga/H-MFI via reaction with [GaH]²⁺ cations to form [GaH₂]⁺-H⁺ cation pairs.