Modeling Proton Transfer in Zeolites: Convergence Behavior of Embedded and Constrained Cluster Calculations

Seamless Multilayer─A Novel Total Energy Partition Scheme for Embedded and Hybrid Calculations

In this paper, we provide general formulation of a multilayer approach, covering both additive and subtractive quantum mechanics/molecular mechanics (QM/MM) as special cases. After that, we suggest a novel definition of QM/MM total energy based on the consideration of a system divided into three layers. In a simplified form, it is E = E Q M ( 1 + 2 ) - E Q M ( 2 ) + E M M ( 2 + 3 ) , where layers 1, 2, and 3 represent inner QM, outer QM, and classical MM regions, respectively. The novel formulation is also not limited by only QM/MM combination of methods─in fact, any computational methods can be combined in a hybrid calculation. In this paper, we call the new approach seamless multilayer. Test calculations performed for silica and boric oxide show that the new approach requires no QM/MM interface parameterization as well as no or very simple correction terms for boundary atoms. This can greatly facilitate QM/MM studies of covalent inorganic solids. However, test calculations of α-Al2O3 show that for ionic compounds, the new method requires some additional development.

The application of QM/MM simulations in heterogeneous catalysis

The QM/MM simulation method is provenly efficient for the simulation of biological systems, where an interplay of extensive environment and delicate local interactions drives a process of interest through a funnel on a complex energy landscape. Recent advances in quantum chemistry and force-field methods present opportunities for the adoption of QM/MM to simulate heterogeneous catalytic processes, and their related systems, where similar intricacies exist on the energy landscape. Herein, the fundamental theoretical considerations for performing QM/MM simulations, and the practical considerations for setting up QM/MM simulations of catalytic systems, are introduced; then, areas of heterogeneous catalysis are explored where QM/MM methods have been most fruitfully applied. The discussion includes simulations performed for adsorption processes in solvent at metallic interfaces, reaction mechanisms within zeolitic systems, nanoparticles, and defect chemistry within ionic solids. We conclude with a perspective on the current state of the field and areas where future opportunities for development and application exist.

Comparison of Catalytic Properties of Vanadium Centers Introduced into BEA Zeolite and Present on (010) V2O5 Surface–DFT Studies

Vanadium-based catalysts, in which vanadium is present either as bulk V2O5 or as isolated species, are active in numerous oxidation reactions. In the present study, vanadium speciation and the possibility of its introduction in various forms (V=O, V–OH, V(=O)(–OH)) into the structurally different crystallographic positions in BEA zeolite was considered by means of Density Functional Theory (DFT). Out of nine nonequivalent positions, T2 and T3 positions are the most preferred. The former may accommodate V=O or V–OH, the latter V–OH or V(=O)(–OH). The structural and electronic properties of all possible centers present in the BEA zeolite are then compared with the characteristics of the same species on the most abundant (010) V2O5 surface. It is demonstrated that they exhibit higher nucleophilic character when introduced into the zeolite, and thus, may be more relevant for catalysis.

Cation pair formation in copper and palladium exchanged MFI zeolite frameworks – a theoretical study

[Pd–O–Pd]²⁺ as a mediator in proton transfer from Brønsted acid sites to reactants; [Cu–OH–Cu]²⁺ stabilized by hydrogen bonds to framework oxygen.

A theoretical study of nitric oxide adsorption and dissociation on copper-exchanged zeolites SSZ-13 and SAPO-34: the impact of framework acid-base properties

The adsorption of nitric oxide as dinitrosyls and the deNOx proton-mediated reaction mechanism are assessed using electronic structure methods and transition state theory. Dinitrosyls bind to copper cations either via a N-atom or via an O-atom, with N-binding being more stable. In their ground states, dinitrosyls reach a planar configuration with the metal cation. The two nitric oxide molecules are kept together in O-bonded dinitrosyls by the N-N bond and the adsorption complex obtains a cyclic planar structure, while N-bonded dinitrosyls have out-of-plane conformations with low energy barriers. An asymmetric structure ZCu(ON)(NO) with one N-bonded nitrosyl and the other O-bonded is of the lowest stability. The cyclic hyponitrite ZCu(ON)2 adsorption complex undergoes O-N bond breaking upon protonation of one oxygen atom and this lowers the energy barrier of the first reaction step of nitric oxide dissociation to yield N2O and a hydroxylated copper site ZCu(OH) by 45 kJ mol(-1) for Cu-SAPO-34 and by 46 kJ mol(-1) for Cu-SSZ-13. The more stable N-bonded dinitrosyl ZCu(NO)2 provides less favorable reaction which passes through the asymmetric ZCu(ON)(NO) intermediate structure. Brønsted acid sites facilitate the reversal of one nitrosyl group. The role of proton transfer from a Brønsted acid site to dinitrosyls is not limited to the initial step of facilitating the N-O bond cleavage, but it also contributes to the stabilization of intermediate oxygen species formed at the copper site as hydroxide ZCu(OH) and hydroperoxide, ZCuOOH. Without protonation, the unstable ZCuO intermediate causes structural deformation with strongly lengthened T-O bonds in the framework. The rate determining step is N2O decomposition to N2 and O2, whether starting with a ZCu(NO)2 or a ZCu(ON)2 adsorption complex, and Cu-SSZ-13 has a clear advantage with an energy barrier of 195 kJ mol(-1)vs. 265 kJ mol(-1) for Cu-SAPO-34. In the final step the Brønsted acid site is restored by proton transfer from the hydroperoxide ZCuOOH to the framework and molecular oxygen is released. The overall energy barrier for the proton-assisted reaction of ZCu(ON)2 decomposition for Cu-SSZ-13 is by 48 kJ mol(-1) lower than the barrier of the proton-free pathway.

Adsorption in zeolites using mechanically embedded ONIOM clusters

We have explored mechanically embedded three-layer QM/QM/MM ONIOM models for computational studies of binding in Al-substituted zeolites.

First principle chemical kinetics in zeolites: the methanol-to-olefin process as a case study

To optimally design next generation catalysts a thorough understanding of the chemical phenomena at the molecular scale is a prerequisite. Apart from qualitative knowledge on the reaction mechanism, it is also essential to be able to predict accurate rate constants. Molecular modeling has become a ubiquitous tool within the field of heterogeneous catalysis. Herein, we review current computational procedures to determine chemical kinetics from first principles, thus by using no experimental input and by modeling the catalyst and reacting species at the molecular level. Therefore, we use the methanol-to-olefin (MTO) process as a case study to illustrate the various theoretical concepts. This process is a showcase example where rational design of the catalyst was for a long time performed on the basis of trial and error, due to insufficient knowledge of the mechanism. For theoreticians the MTO process is particularly challenging as the catalyst has an inherent supramolecular nature, for which not only the Brønsted acidic site is important but also organic species, trapped in the zeolite pores, must be essentially present during active catalyst operation. All these aspects give rise to specific challenges for theoretical modeling. It is shown that present computational techniques have matured to a level where accurate enthalpy barriers and rate constants can be predicted for reactions occurring at a single active site. The comparison with experimental data such as apparent kinetic data for well-defined elementary reactions has become feasible as current computational techniques also allow predicting adsorption enthalpies with reasonable accuracy. Real catalysts are truly heterogeneous in a space- and time-like manner. Future theory developments should focus on extending our view towards phenomena occurring at longer length and time scales and integrating information from various scales towards a unified understanding of the catalyst. Within this respect molecular dynamics methods complemented with additional techniques to simulate rare events are now gradually making their entrance within zeolite catalysis. Recent applications have already given a flavor of the benefit of such techniques to simulate chemical reactions in complex molecular environments.

EFFECT OF HIGH-LEVEL LAYER SELECTION IN ONIOM ON Na+ ADSORPTION IN ZSM-5 ZEOLITE

Here we present a systematical investigation on the effect of high-level layer selection in our own N-layered integrated molecular orbital + molecular mechanics (ONIOM) method. We used the adsorption of Na+ in the M7 site of ZSM-5 as an instructive example. We embedded n tetrahedral sites (nT= 7, 19, 22, 33) for high-level layer into a low-level layer represented by 75T (denoted as nT@75T). Our calculations showed that 7T@75T behaved poorly in both structure optimization and interaction energy calculation. Hence only including the first coordination circle of the ion (i.e. 7T), as commonly used in literature, is not sufficiently good. However, for n ≥ 19, both structure optimization and interaction energy calculation by ONIOM are in good agreement with the full 75T model calculation. Such a reliable ONIOM calculation needs only 1/10 to 1/5 computational time for each self-consistent field (SCF) iteration as compared to the full 75T model calculation, highlighting its pivotal role in zeolite simulation.

QM/MM Modeling of Enantioselective Pybox–Ruthenium- and Box–Copper-Catalyzed Cyclopropanation Reactions: Scope, Performance, and Applications to Ligand Design

An extensive comparison of full-QM (B3LYP) and QM/MM (B3LYP:UFF) levels of theory has been made for two enantioselective catalytic systems, namely, Pybox-Ru and Box-Cu complexes, in the cyclopropanation of alkenes (ethylene and styrene) with methyl diazoacetate. The geometries of the key reaction intermediates and transition structures calculated at the QM/MM level are generally in satisfactory agreement with full-QM calculated geometries. More importantly, the relative energies calculated at the QM/MM level are in good agreement with those calculated at the full-QM level in all cases. Furthermore, the QM/MM energies are often in better agreement with the stereoselectivity experimentally observed, and this suggests that QM/MM calculations can be superior to full-QM calculations when subtle differences in inter- and intramolecular interactions are important in determining the selectivity, as is the case in enantioselective catalysis. The predictive value of the model presented is validated by the explanation of the unusual enantioselectivity behavior exhibited by a new bis-oxazoline ligand, the stereogenic centers of which are quaternary carbon atoms.