Identification of the driving forces in methanol-to-olefin conversion by modeling the zeolite cage and contents

The rapidly increasing demand of oil-based chemicals calls for the development of new technologies based on other natural sources. Among these emerging alternatives, the methanol-to-olefin process (MTO) in acidic zeolites is one of the most promising. However, unraveling the reaction mechanism of such an extremely complex catalytic process like MTO conversion has been a challenging task from both experimental and theoretical viewpoint. For over 30 years the actual mechanism has been one of the most discussed topics in heterogeneous catalysis.[1] Instead of plainly following direct routes,[2-3] the MTO process has experimentally been found to proceed through a hydrocarbon pool mechanism, in which organic reaction centers act as cocatalysts inside the zeolite pores, adding a whole new level of complexity to this issue.[4-5] Therefore, a more detailed understanding of the elementary reaction steps can be obtained with the complementary assistance of theoretical modeling. In this work, a complete supramolecular complex of both the zeolite framework and the co-catalytic hydrocarbon pool species is modeled through state-of-the-art quantum chemical techniques [6-7]. This approach provides a more detailed understanding of the crucial interactions between the zeolite framework and its contents, which form the driving forces for successful methanol-to-olefin conversion. [1] Stocker, M., Microporous Mesoporous Mater. 29 (1999) 3. [2] Song, W.G., Marcus, D.M., Fu, H., Ehresmann, J.O., Haw, J.F., J. Am. Chem. Soc. 124 (2002) 3844. [3] Lesthaeghe, D., Van Speybroeck, V., Marin, G.B., Waroquier, M., Angew. Chem. Int. Ed. 45 (2006) 1714. [4] Dessau, R. M., J. Catal. 99 (1986) 111. [5] Dahl, I.M., Kolboe, S., Catal. Lett. 20 (1993) 329. [6] Lesthaeghe, D., De Sterck, B., Van Speybroeck, V., Marin, G.B., Waroquier, M., Angew. Chem. Int. Ed. 46 (2007) 1311. [7] McCann, D.M., Lesthaeghe, D., Kletnieks, P.W., Guenther, D.R., Hayman, M.J., Van Speybroeck, V., Waroquier, M., Haw, J.F. Angew. Chem. Int. Ed. 47 (2008) 5179

Theoretical identification of the interactions between the zeolite framework and the hydrocarbon pool co-catalyst in methanol-to-olefin conversion

The rapidly increasing demand of oil-based chemicals calls for the development of new technologies based on other natural sources. Among these emerging alternatives, the methanol-to-olefin process (MTO) in acidic zeolites is one of the most promising. However, unraveling the reaction mechanism of such an extremely complex catalytic process like MTO conversion has been a challenging task from both experimental and theoretical viewpoint. For over 30 years the actual mechanism has been one of the most discussed topics in heterogeneous catalysis.[1] Instead of plainly following direct routes,[2-3] the MTO process has experimentally been found to proceed through a hydrocarbon pool mechanism, in which organic reaction centers act as homogeneous co-catalysts inside the heterogeneous acid catalyst, adding a whole new level of complexity to this issue.[4-5] Therefore, a more detailed understanding of the elementary reaction steps can be obtained with the complementary assistance of theoretical modeling. In this work, a complete supramolecular complex of both the zeolite framework and the co-catalytic hydrocarbon pool species is modeled through state-of-the-art quantum chemical techniques [6-7]. This approach provides a more detailed understanding of the crucial interactions between the zeolite framework and its contents, which form the driving forces for successful methanol-to-olefin conversion

The exclusive (e,e′'p) reaction at high missing momenta

The reduced (e,e$'$p) cross section is calculated for kinematics that probe high missing momenta. The final-state interaction is handled within a non-relativistic many-body framework. One- and two-body nuclear currents are included. Electron distortion effects are treated in an exact distorted wave calculation. It is shown that at high missing momenta the calculated (e,e$'$p) cross sections exhibit a pronounced sensitivity to ground-state correlations of the RPA type and two-body currents. The role of these mechanisms is found to be relatively small at low missing momenta.Comment: 15 pages in REVtex with embedded psfigure

Monte Carlo simulations to understand 'breathing' phenomenon of metal organic frameworks

Metal Organic Frameworks (MOFs) are a new class of porous materials synthesized from metal clusters connected by organic linkers. One of the promising applications of MOFs is carbon capture from fuel gasses, where CO2 is adsorbed in the pores of the material. In this presentation, we explore framework flexibility as a possible mechanism for selective and reversible CO2 adsorption by means of Monte Carlo simulations. Most MOFs are fairly rigid structures, in the sense that they undergo small changes in volume when external stress is applied. Typical volume changes are of the order of a few percent only. Nevertheless, some MOF materials have an unexpectedly high flexibility and impressively shrink or swell under pressure, temperature or adsorption changes. A well-known example is MIL-53, a structure that shows volume changes of over 40%. In an adsorption experiment, the gas pressure is gradually increased while the amount of adsorbed material in the pores is measured. For MIL-53, the measured adsorption isotherm shows interesting features: when MIL-53 is brought into contact with a gas at increasing pressure, the framework's pores constrict, while at even higher pressures, the pores return to their original geometry. The process, referred to as "breathing", is reversible and shows hysteresis. Based on Monte Carlo runs, we have constructed a mean-field model to gain insight in the thermodynamics of the breathing. The model shows that the behavior is the result of the different factors at play in a (Nmof,μ,P,T) ensemble (constant amount of MOF material, constant gas chemical potential, constant gas pressure, constant temperature), i.e. the entropy, the pressure and the resistance given by the adsorbed particles. We further investigate how the MOFs' flexibility could be exploited to design an efficient pressure swing setup

The Monomer Electron Density Force Field (MEDFF) : a physically inspired model for noncovalent interactions

We propose a methodology to derive pairwise-additive noncovalent force fields from monomer electron densities without any empirical input. Energy expressions are based on the symmetry-adapted perturbation theory (SAPT) decomposition of interaction energies. This ensures a physically motivated force field featuring an electrostatic, exchange repulsion, dispersion, and induction contribution, which contains two types of parameters. First, each contribution depends on several fixed atomic parameters, resulting from a partitioning of the monomer electron density. Second, each of the last three contributions (exchange-repulsion, dispersion, and induction) contains exactly one linear fitting parameter. These three so-called interaction parameters in the model are initially estimated separately using SAPT reference calculations for the S66x8 database of noncovalent dimers. In a second step, the three interaction parameters are further refined simultaneously to reproduce CCSD(T)/CBS interaction energies for the same database. The limited number of parameters that are fitted to dimer interaction energies (only three) avoids ill-conditioned fits that plague conventional parameter optimizations. For the exchange repulsion and dispersion component, good results are obtained for all dimers in the S66x8 database using one single value for the associated interaction parameters. The values of those parameters can be considered universal and can also be used for dimers not present in the original database used for fitting. For the induction component such an approach is only viable for the dispersion dominated dimers in the S66x8 database. For other dimers (such as hydrogen-bonded complexes), we show that our methodology remains applicable. However, the interaction parameter needs to be determined on a case-specific basis. As an external validation:, the force field predicts interaction energies in good agreement with CCSD(T)/CBS values for dispersion dominated dimers extracted from an HIV-II protease crystal structure with a bound ligand (indinavir). Furthermore, experimental second virial coefficients of small alkanes and alkenes are well reproduced

Deactivation of the catalyst during the MTO process from a molecular modeling perspective

Currently, the industrially important conversion process of methanol to olefins (MTO) forms a key process for the production of higher valued products that can easily be transported, such as ethylene and propylene. Methanol can be made from natural gas or coal via synthesis gas. Unraveling the underlying reaction mechanism of the complex MTO process has already shown to be very challenging. Recent ab initio calculations, in combination with experimental data, are in strong support of the “hydrocarbon pool model” as opposed to a direct (C-C coupling) route [1, 2]. The hydrocarbon pool has been described as a catalytic scaffold inside the zeolite building, consisting of polymethylbenzenes and their cationic derivatives. The continued growth of these initially active carbonaceous species within acidic zeolites, such as H-ZSM-5 and H-SAPO-34, is an undesired side effect resulting from secondary reactions for which at present no computational data exist whatsoever. The presence of these large species – coke precursors - inside or at the external cups of the periodic structure leads to blockage of the pores or channels and ultimately to the deactivation of the catalyst. An improved in-depth understanding of the underlying reaction mechanisms of coke formation is therefore desperately needed. A main problem is the generally poor characterization of coke, despite the great number of techniques (gas chromatography, mass spectroscopy) that can be used for locating and identifying the deposits [3, 4]. Because of this, it is not clear whether benzenoid species consisting of 3 rings can already be regarded as coke as opposed to large aromatic species present in the hydrocarbon pool that still allow an active route. Within this contribution possible reaction routes leading to the formation of naphthalene- and/or phenanthrene-like species are studied from theoretical viewpoint within various industrially relevant zeolite topologies. For each of these elementary steps reaction rates are evaluated based on energies and frequencies originating from reliable ab initio data. The latter were obtained by taking into account a large portion of the zeolites, as to be representative for the actual topology

Modelling of Lewis-Acid Catalyzed Ring Opening of Oxanorbornenes in the Synthesis of Azaheterocyclic Phosphonates

Since the discovery of the biological activity of aminophosphonates, research started on the synthesis of more constraint azaheterocyclic phosphonates. We developed a route via an intramolecular Diels-Alder reaction towards α-aminophosphonates 1. The obtained oxanorbornene skeleton is a valuable synthetic intermediate that has been used in various natural product syntheses. An important synthetic transformation involves the cleavage of the oxygen bridge, used to construct substituted arenes and cyclohexenes. We wanted to investigate the ring opening of adducts 1 using different Lewis acids experimentally and to get more insight in the reaction pathways towards the different products via molecular modelling. In this presentation the results obtained with TiCl4 and FeCl3 catalyst are shown. One of the difficulties in studying transition metal catalysts is the determination of their proper spin state. The tetrahedral TiCl4 monomer has spin zero. The FeCl3 catalyst has a high-spin ground state. It prefers a half-filled d-shell and has multiplicity 6. The complexation of the Lewis acids with different binding sites was investigated at a B3LYP level of theory with a LanL2DZ pseudopotential for the transition metals. Bidentate coordination towards the most electronegative phosphonate oxygen and the oxygen bridge is favoured for both catalysts. The reaction pathways were evaluated at a TPSSh and a B3LYP level of theory. The role of dispersion interactions was evaluated using the Van der Waals correction term from B3LYP-D. The energy barrier for breaking the C-O bond with FeCl3 is larger than with TiCl4. This corresponds with the experimental observation that the titanium catalyzed reaction completes at 0°C and the reaction with the iron catalyst requires reflux conditions in CH2Cl2. The main difference between the TiCl4 and FeCl3 as Lewis acids in the opening of the oxanorbornene oxygen bridge is their way of stabilizing the oxide anion. When the C-O bond is broken, the bond between the alkoxide anion and the transition metal tightens. With TiCl4, the alkoxide replaces a chloride that adds to the allyl cation in a concerted way. With FeCl3, however, the carbocation is stabilized by the alkoxide anion itself and no chloride transfer occurs; instead, the bond with the phosphonate is broken. A plausible further reaction path towards the experimentally observed ketone involves a 1,2-hydride shift