QM/MM study of the stability of dimethyl ether in zeolites H-ZSM-5 and H-Y

The methanol-to-hydrocarbons (MTH) process transforms C1 carbon sources to higher hydrocarbons, but details of the mechanism that leads to the formation of the first carbon–carbon bond remain unclear. Here, we present a computational investigation of how a crucial intermediate, dimethyl ether (DME), interacts with different zeolite catalysts (H-ZSM-5, H-Y) to gain insight into the initial stages in the MTH process. We use QM/MM computational simulations to model the conversion of methanol to DME in H-ZSM-5, which is a well characterised and important reaction intermediate. We analyse and compare the stability of DME on several acid sites in H-ZSM-5 and H-Y, and show that the more acidic and open “intersection sites” in the H-ZSM-5 framework are able to bond strongest with DME, with complete deprotonation of the acid site occurring. The conversion of methanol to DME in H-ZSM-5 is calculated as requiring a higher activation energy than framework methoxylation, which indicates that a stepwise (indirect) mechanism, through a methoxy intermediate, is the most likely route to DME formation during the initiation of the MTH process

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

QM/MM study of the reactivity of zeolite bound methoxy and carbene groups

The conversion of methanol-to-hydrocarbons (MTH) is known to occur via an autocatalytic process in zeolites, where framework-bound methoxy species play a pivotal role, especially during catalyst induction. Recent NMR and FT-IR experimental studies suggest that methoxylated zeolites are able to produce hydrocarbons by a mechanism involving carbene migration and association. In order to understand these observations, we have performed QM/MM computational investigations on a range of reaction mechanisms for the reaction of zeolite bound methoxy and carbene groups, which are proposed to initiate hydrocarbon formation in the MTH process. Our simulations demonstrate that it is kinetically unfavourable for methyl species to form on the framework away from the zeolite acid site, and both kinetically and thermodynamically unfavourable for methyl groups to migrate through the framework and aggregate around an acid site. Formation of carbene moieties was considered as an alternative pathway to the formation of C–C bonds; however, the reaction energy for conversion of a methyl to a carbene is unfavourable. Metadynamics simulations help confirm further that methyl species at the framework acid sites would be more reactive towards formed C_{2+} species, rather than inter-framework migration, and that the role of carbenes in the formation of the first C–C bond will be via a concerted type of mechanism rather than stepwise

Magnetic coupling constants for MnO as calculated using hybrid density functional theory

The properties of MnO have been calculated using generalised gradient approximation (GGA-) and hybrid (h-) density functional theory (DFT), specifically variants of the popular PBE and PBESol exchange–correlation functionals. The GGA approaches are shown to be poor at reproducing experimental magnetic coupling constants and rhombohedral structural distortions, with the PBESol functional performing worse than PBE. In contrast, h-DFT results are in reasonable agreement with experiment. Calculation of the Néel temperatures using the mean-field approximation gives overestimates relative to experiment, but the discrepancies are as low as 15 K for the PBE0 approach and, generally, the h-DFT results are significant improvements over previous theoretical studies. For the Curie–Weiss temperature, larger disparities are observed between the theoretical results and previous experimental results

Methanol loading dependent methoxylation in zeolite H-ZSM-5

We evaluate the effect of the number of methanol molecules per acidic site of H-ZSM-5 on the methoxylation reaction at room temperature by applying operando diffuse reflectance infrared Fourier transformed spectroscopy (DRIFTS) and mass spectrometry (MS), which capture the methoxylation reaction by simultaneously probing surface adsorbed species and reaction products, respectively. To this end, the methanol loading in H-ZSM-5 (Si/Al z 25) pores is systematically varied between 32, 16, 8 and 4 molecules per unit cell, which corresponds to 8, 4, 2 and 1 molecules per Brønsted acidic site, respectively. The operando DRIFTS/MS data show that the room temperature methoxylation depends on the methanol loading: the higher the methanol loading, the faster the methoxylation. Accordingly, the reaction is more than an order of magnitude faster with 8 methanol molecules per Brønsted acidic site than that with 2 molecules, as evident from the evolution of the methyl rock band of the methoxy species and of water as a function of time. Significantly, no methoxylation is observed with #1 molecule per Brønsted acidic site. However, hydrogen bonded methanol occurs across all loadings studied, but the structure of hydrogen bonded methanol also depends on the loading. Methanol loading of #1 molecule per acidic site leads to the formation of hydrogen bonded methanol with no proton transfer (i.e. neutral geometry), while loading $2 molecules per acidic site results in a hydrogen bonded methanol with a net positive charge on the adduct (protonated geometry). The infrared vibrational frequencies of methoxy and hydrogen bonded methanol are corroborated by Density Functional Theory (DFT) calculations. Both the experiments and calculations reflect the methoxy bands at around 940, 1180, 2868–2876 and 2980–2973 cm�1 which correspond to n(C–O), r(CH3), ns(C–H) and nas(C–H), respectively

Computational investigation of CO adsorbed on Au-x, Ag-x and (AuAg)(x) nanoclusters (x=1-5, 147) and monometallic Au and Ag low-energy surfaces

Density functional theory calculations have been performed to investigate the use of CO as a probe molecule for the determination of the structure and composition of Au, Ag and AuAg nanoparticles. For very small nanoclusters (x = 1 − 5), the CO vibrational frequencies can be directly correlated to CO adsorption strength, whereas larger 147-atom nanoparticles show a strong energetic preference for CO adsorption at a vertex position but the highest wavenumbers are for the bridge positions. We also studied CO adsorption on Au and Ag (100) and (111) surfaces, for a 1 monolayer coverage, which proves to be energetically favourable on atop only and bridge positions for Au (100) and atop for Ag (100); vibrational frequencies of the CO molecules red-shift to lower wavenumbers as a result of increased metal coordination. We conclude that CO vibrational frequencies cannot be solely relied upon in order to obtain accurate compositional analysis, but we do propose that elemental rearrangement in the core@shell nanoclusters, from Ag@Au (or Au@Ag) to an alloy, would result in a shift in the CO vibrational frequencies that indicate changes in the surface composition

Bulk ionization potentials and band alignments from three-dimensional periodic calculations as demonstrated on rocksalt oxides

The position of the band edges of a material plays a key role in determining the properties for a range of applications, but fundamental band bending is an interface-dependent property that cannot be quantified without knowledge of bulk electron energy levels. We present a method for calculating the bulk position of the valence band maximum, and therefore the bulk ionization potential, from periodic plane wave calculations as shown for a range of rocksalt ionic oxides. We demonstrate that, for the popular “slab alignment” technique, explicit consideration of any surface induced electronic polarization is necessary to calculate accurate bulk ionization potentials. Our proposed method to quantify these surface effects, using polarizable-shell based interatomic potentials, is very computationally affordable, and our updated slab alignment method yields much improved agreement with the available experimental data

Segregation effects on the properties of (AuAg)₁₄₇.

AuAg nanoclusters are promising supported co-catalysts for photocatalytic hydrogen reduction. However, beyond the quantum regime (N > 100) little is known about how the electronic properties of these nanoparticles are affected by chemical ordering. We investigate the effects of chemical ordering on the properties of 147-atom cuboctahedral AuAg nanoclusters, using empirical potentials coupled with an atomic-swap basin-hopping search to optimise the elemental distribution, with the lowest energy arrangements then reminimised using Density Functional Theory (DFT). Force-field calculations show Au atoms preferentially occupy sub-surface positions in the bimetallic structures, which results in the formation of a pseudo-onion structure for Ag-rich compositions. At the DFT-level, however, an Ag core surrounded by an Au shell (Ag@Au) is energetically favoured, as electron density can be drawn more readily when Au atoms are positioned on the nanocluster surface, thus resulting in a partial negative charge. Core@shell configurations are analogous to structures that can be chemically synthesised, and further detailed electronic analysis is discussed in the context of nanocluster applications to co-catalysed photocatalysis

Morphological Features and Band Bending at Nonpolar Surfaces of ZnO

We employ hybrid density functional calculations to analyze the structure and stability of the (101̅0) and (112̅0) ZnO surfaces, confirming the relative stability of the two surfaces. We then examine morphological features, including steps, dimer vacancies, and grooves, at the main nonpolar ZnO surface using density functional methods. Calculations explain why steps are common on the (101̅0) surface even at room temperature, as seen in experiment. The surface structure established has been used to obtain the definitive ionization potential and electron affinity of ZnO in good agreement with experiment. The band bending across the surface is analyzed by the decomposition of the density of states for each atomic layer. The upward surface band bending at the (101̅0) surface affects mostly the valence band by 0.32 eV, which results in the surface band gap closing by 0.31 eV; at the (112̅0) surface, the valence band remains flat and the conduction band bends up by 0.18 eV opening the surface band gap by 0.12 eV

Tailoring Gold Nanoparticle Characteristics and the Impact on Aqueous-Phase Oxidation of Glycerol

Poly(vinyl alcohol) (PVA)-stabilized Au nanoparticles (NPs) were synthesized by colloidal methods in which temperature variations (−75 to 75 °C) and mixed H2O/EtOH solvent ratios (0, 50, and 100 vol/vol) were used. The resulting Au NPs were immobilized on TiO2 (P25), and their catalytic performance was investigated for the liquid phase oxidation of glycerol. For each unique solvent system, there was a systematic increase in the average Au particle diameter as the temperature of the colloidal preparation increased. Generation of the Au NPs in H2O at 1 °C resulted in a high observed activity compared with current Au/TiO2 catalysts (turnover frequency = 915 h–1). Interestingly, Au catalysts with similar average particle sizes but prepared under different conditions had contrasting catalytic performance. For the most active catalyst, aberration-corrected high angle annular dark field scanning transmission electron microscopy analysis identified the presence of isolated Au clusters (from 1 to 5 atoms) for the first time using a modified colloidal method, which was supported by experimental and computational CO adsorption studies. It is proposed that the variations in the populations of these species, in combination with other solvent/PVA effects, is responsible for the contrasting catalytic properties