Modelling metal centres, acid sites and reaction mechanisms in microporous catalysts.

We discuss the role of QM/MM (embedded cluster) computational techniques in catalytic science, in particular their application to microporous catalysis. We describe the methodologies employed and illustrate their utility by briefly summarising work on metal centres in zeolites. We then report a detailed investigation into the behaviour of methanol at acidic sites in zeolites H-ZSM-5 and H-Y in the context of the methanol-to-hydrocarbons/olefins process. Studying key initial steps of the reaction (the adsorption and subsequent methoxylation), we probe the effect of framework topology and Brønsted acid site location on the energetics of these initial processes. We find that although methoxylation is endothermic with respect to the adsorbed system (by 17-56 kJ mol(-1) depending on the location), there are intriguing correlations between the adsorption/reaction energies and the geometries of the adsorbed species, of particular significance being the coordination of methyl hydrogens. These observations emphasise the importance of adsorbate coordination with the framework in zeolite catalysed conversions, and how this may vary with framework topology and site location, particularly suited to investigation by QM/MM techniques.

[1]  C. Catlow,et al.  Computational studies of the reaction of methanol at aluminosilicate Brønsted acid sites , 1996 .

[2]  Ian J. Bush,et al.  The GAMESS-UK electronic structure package: algorithms, developments and applications , 2005 .

[3]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[4]  S. C. Rogers,et al.  QUASI: A general purpose implementation of the QM/MM approach and its application to problems in catalysis , 2003 .

[5]  J. Sauer,et al.  MOLECULAR MECHANICS POTENTIAL FOR SILICA AND ZEOLITE CATALYSTS BASED ON AB INITIO CALCULATIONS 2 , 1995 .

[6]  P. M. Rodger,et al.  DL_POLY: Application to molecular simulation , 2002 .

[7]  C. Catlow,et al.  On the structure and coordination of the oxygen-donating species in Ti↑MCM-41/TBHP oxidation catalysts: a density functional theory and EXAFS study , 2002 .

[8]  Travis V. Harris,et al.  The ONIOM Method and Its Applications. , 2015, Chemical reviews.

[9]  S. H. Vosko,et al.  Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis , 1980 .

[10]  J. Sauer,et al.  Ab initio molecular dynamics simulation of methanol interacting with acidic zeolites of different framework structure , 2000 .

[11]  Richard J. Hall,et al.  Adsorption energies of NH3 and NH4+ in zeolites. An embedded cluster model including electron correlation , 1995 .

[12]  P. Sherwood,et al.  Embedded Cluster Model for the ab Initio Study of Broensted Acidity in Zeolites , 1994 .

[13]  Comparison of embedded cluster models to study zeolite catalysis: Proton transfer reactions in acidic chabazite , 1997 .

[14]  C. Catlow,et al.  Density functional theory calculations of adsorption and reactivity of methanol at alumino-silicate Brønsted acid centres , 1997 .

[15]  E. Iglesia,et al.  The Strength of Brønsted Acid Sites in Microporous Aluminosilicates , 2015 .

[16]  Marek Sierka,et al.  Structure and reactivity of silica and zeolite catalysts by a combined quantum mechanics[ndash ]shell-model potential approach based on DFT , 1997 .

[17]  P. Sherwood,et al.  Hybrid QM/MM embedding approach for the treatment of localized surface states in ionic materials , 2004 .

[18]  C. Catlow,et al.  Room temperature methoxylation in zeolites: insight into a key step of the methanol-to-hydrocarbons process. , 2016, Chemical communications.

[19]  Wang,et al.  Accurate and simple analytic representation of the electron-gas correlation energy. , 1992, Physical review. B, Condensed matter.

[20]  Alexey A. Sokol,et al.  ChemShell—a modular software package for QM/MM simulations , 2014 .

[21]  M. Payne,et al.  Methanol Adsorption in ZeolitesA First-Principles Study , 1996 .

[22]  K. Lillerud,et al.  Conversion of methanol to hydrocarbons: how zeolite cavity and pore size controls product selectivity. , 2012, Angewandte Chemie.

[23]  P. Sherwood,et al.  Adsorption of water and methanol on zeolite Bro/nsted acid sites: An ab initio, embedded cluster study including electron correlation , 1996 .

[24]  A. Maiti,et al.  DFT study of methanol conversion to hydrocarbons in a zeolite catalyst , 2003 .

[25]  C. Catlow,et al.  The architecture of catalytically active centers in titanosilicate (TS-1) and related selective-oxidation catalysts , 2000 .

[26]  Arthur H. Edwards,et al.  Embedding techniques for irradiation-induced defects in crystalline SiO 2 , 2001 .

[27]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[28]  M. Frisch,et al.  Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields , 1994 .

[29]  P. Sherwood,et al.  Formation of heteroatom active sites in zeolites by hydrolysis and inversion. , 2006, Angewandte Chemie.

[30]  A. Schäfer,et al.  Absolute acidities and site specific properties of zeolite catalysts modelled by advanced computational chemistry technology , 1999 .

[31]  S. Blaszkowski,et al.  Theoretical Study of the Mechanism of Surface Methoxy and Dimethyl Ether Formation from Methanol Catalyzed by Zeolitic Protons , 1997 .

[32]  A. Corma,et al.  Formation of Surface Methoxy Groups on H-Zeolites from Methanol. A Quantum Chemical Study , 1995 .

[33]  C. Lamberti,et al.  The Role of Isolated Sites in Heterogeneous Catalysis: Characterization and Modeling , 2001 .

[34]  Lamberti,et al.  Neutron powder diffraction study of orthorhombic and monoclinic defective silicalite , 2000, Acta crystallographica. Section B, Structural science.

[35]  C. Catlow,et al.  Structure and dissociation mechanisms of methanol in ZSM-5 zeolite , 1989 .

[36]  M. Vincent,et al.  Computer simulation of zeolite structure and reactivity using embedded cluster methods , 1997 .

[37]  Joachim Sauer,et al.  Molecular mechanics potential for silica and zeolite catalysts based on ab initio calculations. 1. Dense and microporous silica , 1994 .

[38]  D. Nachtigallová,et al.  Calculations of the site specific stretching frequencies of CO adsorbed on Li+/ZSM-5 , 2004 .

[39]  David J. Tozer,et al.  Hybrid exchange-correlation functional determined from thermochemical data and ab initio potentials , 2001 .

[40]  A. Cheetham,et al.  Powder Neutron Diffraction and 29Si MAS NMR Studies of Siliceous Zeolite-Y , 1993 .

[41]  V. Van Speybroeck,et al.  Unraveling the reaction mechanisms governing methanol-to-olefins catalysis by theory and experiment. , 2013, Chemphyschem : a European journal of chemical physics and physical chemistry.

[42]  M. Sierka,et al.  Heterogeneity of Brønsted Acidic Sites in Faujasite Type Zeolites due to Aluminum Content and Framework Structure , 1998 .

[43]  J. Carruthers,et al.  An ab initio study of methanol adsorption in zeolites , 1993 .

[44]  A. Edwards,et al.  Calibration of embedded-cluster method for defect studies in amorphous silica , 2004 .

[45]  C. Catlow,et al.  Active sites for heterogeneous catalysis by functionalisation of internal and external surfaces , 2004 .

[46]  Notker Rösch,et al.  Elastic Polarizable Environment Cluster Embedding Approach for Covalent Oxides and Zeolites Based on a Density Functional Method , 2003 .

[47]  S. French,et al.  Formation of active sites in TS-1 by hydrolysis and inversion , 2007 .

[48]  Michel Waroquier,et al.  The Rise and Fall of Direct Mechanisms in Methanol-to-Olefin Catalysis: An Overview of Theoretical Contributions , 2007 .

[49]  D. Goldfarb.  A family of variable-metric methods derived by variational means , 1970 .

[50]  C. G. Broyden.  The convergence of single-rank quasi-Newton methods , 1970 .

[51]  J. Hutter,et al.  Ab initio molecular dynamics simulation of methanol adsorbed in chabazite , 1997 .

[52]  J. Sauer,et al.  Predicting Absolute and Site Specific Acidities for Zeolite Catalysts by a Combined Quantum Mechanics/Interatomic Potential Function Approach† , 1997 .

[53]  Paul Sherwood,et al.  Embedded-cluster calculations in a numeric atomic orbital density-functional theory framework. , 2014, The Journal of chemical physics.

[54]  Peter R. Taylor,et al.  The choice of Gaussian basis sets for molecular electronic structure calculations , 1981 .

[55]  Michel Waroquier,et al.  Advances in theory and their application within the field of zeolite chemistry. , 2015, Chemical Society reviews.

[56]  D. Shanno.  Conditioning of Quasi-Newton Methods for Function Minimization , 1970 .

[57]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[58]  S. Blaszkowski,et al.  The mechanism of dimethyl ether formation from methanol catalyzed by zeolitic protons , 1996 .

[59]  C. Catlow,et al.  Hybrid QM/MM Investigations into the Structure and Properties of Oxygen-Donating Species in TS-1 , 2008 .

[60]  S. Blaszkowski,et al.  Density functional theory calculations of the activation of methanol by a Bronsted zeolitic proton , 1995 .

[61]  M. Payne,et al.  Methanol in microporous materials from first principles , 1999 .

[62]  Weiguo Song,et al.  The Role of Methoxy Groups in Methanol to Olefin Conversion , 2008 .

[63]  R. Fletcher,et al.  A New Approach to Variable Metric Algorithms , 1970, Comput. J..

[64]  Alexey A. Sokol,et al.  QM/MM modelling of the TS-1 catalyst using HPCx , 2006 .

[65]  J. Sauer,et al.  Quantum chemical modeling of zeolite-catalyzed methylation reactions: toward chemical accuracy for barriers. , 2009, Journal of the American Chemical Society.

[66]  J. Sauer,et al.  Bridging and terminal hydroxyls. A structural chemical and quantum chemical discussion , 1984 .

[67]  P. Ratnasamy,et al.  Reactive oxygen species in epoxidation reactions over titanosilicate molecular sieves , 2003 .

[68]  W. Jost.  Diffusion and Electrolytic Conduction in Crystals (Ionic Semiconductors) , 1933 .

[69]  Clarence Dayton Chang,et al.  The conversion of methanol and other O-compounds to hydrocarbons over zeolite catalysts: II. Pressure effects , 1977 .

[70]  C. Maiden.  The New Zealand Gas-to-Gasoline Project , 1988 .