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

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

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

A computational study of direct CO₂ hydrogenation to methanol on Pd surfaces

The reaction mechanism of direct CO2 hydrogenation to methanol is investigated in detail on Pd (111), (100) and (110) surfaces using density functional theory (DFT), supporting investigations into emergent Pd-based catalysts. Hydrogen adsorption and surface mobility are firstly considered, with high-coordination surface sites having the largest adsorption energy and being connected by diffusion channels with low energy barriers. Surface chemisorption of CO2, forming a partially charged CO2δ−, is weakly endothermic on a Pd (111) whilst slightly exothermic on Pd (100) and (110), with adsorption enthalpies of 0.09, −0.09 and −0.19 eV, respectively; the low stability of CO2δ− on the Pd (111) surface is attributed to negative charge accumulating on the surface Pd atoms that interact directly with the CO2δ− adsorbate. Detailed consideration for sequential hydrogenation of the CO2 shows that HCOOH hydrogenation to H2COOH would be the rate determining step in the conversion to methanol, for all surfaces, with activation barriers of 1.41, 1.51, and 0.84 eV on Pd (111), (100) and (110) facets, respectively. The Pd (110) surface exhibits overall lower activation energies than the most studied Pd (111) and (100) surfaces, and therefore should be considered in more detail in future Pd catalytic studies

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

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 ≈ 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 ν(C–O), ρ(CH3), νs(C–H) and νas(C–H), respectively

Controlling structural transitions in AuAg nanoparticles through precise compositional design

We present a study of the transitional pathways between highsymmetry structural motifs for AgAu nanoparticles, with a specific focus on controlling the energetic barriers through chemical design. We show that the barriers can be altered by careful control of the elemental composition and chemical arrangement, with core@shell and vertex-decorated arrangements being specifically influential on the barrier heights. We also highlight the complexity of the potential and free energy landscapes for systems where there are low-symmetry geometric motifs that are energetically competitive to the high-symmetry arrangements. In particular, we highlight that some core@shell arrangements preferentially transition through multistep restructuring of lowsymmetry truncated octahedra and rosette-icosahedra, instead of via the more straightforward square-diamond transformations, due to lower energy barriers and competitive energetic minima. Our results have promising implications for the continuing efforts in bespoke nanoparticle design for catalytic and plasmonic applications

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

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

Investigation of the Pd(1−x)Znx alloy phase diagram using ab initio modelling approaches

The identification of the stable phases in alloy materials is challenging because composition affects the structural stability of different intermediate phases. Computational simulation, via multiscale modelling approaches, can significantly accelerate the exploration of phase space and help to identify stable phases. Here, we apply such new approaches to understand the complex phase diagram of binary alloys of PdZn, with the relative stability of structural polymorphs considered through application of density functional theory coupled with cluster expansion (CE). The experimental phase diagram has several competing crystal structures, and we focus on three different closed-packed phases that are commonly observed for PdZn, namely the face-centred cubic (FCC), body-centred tetragonal (BCT) and hexagonal close packed (HCP), to identify their respective stability ranges. Our multiscale approach confirms a narrow range of stability for the BCT mixed alloy, within the Zn concentration range from 43.75% to 50%, which aligns with experimental observations. We subsequently use CE to show that the phases are competitive across all concentrations, but with the FCC alloy phase favoured for Zn concentrations below 43.75%, and that the HCP structure favoured for Zn-rich concentrations. Our methodology and results provide a platform for future investigations of PdZn and other close-packed alloy systems with multiscale modelling techniques