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

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

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

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

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

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

Embedded-Cluster Calculations in a Numeric Atomic Orbital Density-Functional Theory Framework

We integrate the all-electron electronic structure code FHI-aims into the general ChemShell package for solid-state embedding (QM/MM) calculations. A major undertaking in this integration is the implementation of pseudopotential functionality into FHI-aims to describe cations at the QM/MM boundary through effective core potentials and therewith prevent spurious overpolarization of the electronic density. Based on numeric atomic orbital basis sets, FHI-aims offers particularly efficient access to exact exchange and second order perturbation theory, rendering the established QM/MM setup an ideal tool for hybrid and double-hybrid level DFT calculations of solid systems. We illustrate this capability by calculating the reduction potential of Fe in the Fe-substituted ZSM-5 zeolitic framework and the reaction energy profile for (photo-)catalytic water oxidation at TiO2(110).Comment: 12 pages, 4 figure

Modelling the chemistry of Mn-doped MgO for bulk and (100) surfaces

We have investigated the energetic properties of Mn-doped MgO bulk and (100) surfaces using a QM/MM embedding computational method, calculating the formation energy for doped systems, as well as for surface defects, and the subsequent effect on chemical reactivity. Low-concentration Mn doping is endothermic for isovalent species in the bulk but exothermic for higher oxidation states under p-type conditions, and compensated by electrons going to the Fermi level rather than cation vacancies. The highest occupied dopant Mn 3d states are positioned in the MgO band gap, about 4.2 eV below the vacuum level. Surface Mn-doping is more favourable than subsurface doping, and marginally exothermic on a (100) surface at high O2 pressures. For both types of isovalent Mn-doped (100) surfaces, the formation energy for catalytically important oxygen defects is less than for pristine MgO, with F0 and F2+-centres favoured in n- and p-type conditions, respectively. In addition, F+-centres are stabilised by favourable exchange coupling between the Mn 3d states and the vacancy-localised electrons, as verified through calculation of the vertical ionisation potential. The adsorption of CO2 on to the pristine and defective (100) surface is used as a probe of chemical reactivity, with isovalent subsurface Mn dopants mildly affecting reactivity, whereas isovalent surface-positioned Mn strongly alters the chemical interactions between the substrate and adsorbate. The differing chemical reactivity, when compared to pristine MgO, justifies further detailed investigations for more varied oxidation states and dopant species

Deep vs shallow nature of oxygen vacancies and consequent n -type carrier concentrations in transparent conducting oxides

The source of n -type conductivity in undoped transparent conducting oxides has been a topic of debate for several decades. The point defect of most interest in this respect is the oxygen vacancy, but there are many conflicting reports on the shallow versus deep nature of its related electronic states. Here, using a hybrid quantum mechanical/molecular mechanical embedded cluster approach, we have computed formation and ionization energies of oxygen vacancies in three representative transparent conducting oxides: In 2 O 3 , SnO 2 , and ZnO. We find that, in all three systems, oxygen vacancies form well-localized, compact donors. We demonstrate, however, that such compactness does not preclude the possibility of these states being shallow in nature, by considering the energetic balance between the vacancy binding electrons that are in localized orbitals or in effective-mass-like diffuse orbitals. Our results show that, thermodynamically, oxygen vacancies in bulk In 2 O 3 introduce states above the conduction band minimum that contribute significantly to the observed conductivity properties of undoped samples. For ZnO and SnO 2 , the states are deep, and our calculated ionization energies agree well with thermochemical and optical experiments. Our computed equilibrium defect and carrier concentrations, however, demonstrate that these deep states may nevertheless lead to significant intrinsic n -type conductivity under reducing conditions at elevated temperatures. Our study indicates the importance of oxygen vacancies in relation to intrinsic carrier concentrations not only in In 2 O 3 , but also in SnO 2 and ZnO

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 technique

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