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

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

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

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

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

Evaluating the Role of Anharmonic Vibrations in Zeolite β Materials

The characterization of zeolitic materials is often facilitated by spectroscopic analysis of vibrations, which informs about the bonding character of the substrate and any adsorbents. Computational simulations aid the interpretation of the spectra but often ignore anharmonic effects that can affect the spectral characteristics significantly. Here, the impact of anharmonicity is demonstrated with a combination of dynamical and static simulations applied to the structures formed during the synthesis of Sn-BEA via solid-state incorporation (SSI): the initial siliceous BEA (Si-β), aluminosilicate BEA (H-β), dealuminated BEA (deAl-β), and Sn-BEA (Sn-β). Heteroatom and defect-containing BEA are shown to have strong anharmonic vibrational contributions, with atomic and elemental resolution highlighting particularly the prevalence for H atoms (H-β, deAl-β) as well as localization to heteroatoms at defect sites. We simulate the vibrational spectra of BEA accounting for anharmonic contributions and observe an improved agreement with experimental data compared to harmonic methods, particularly at wavenumbers below 1500 cm–1. The results demonstrate the importance of incorporating anharmonic effects in simulations of vibrational spectra, with consequences toward future characterization and application of zeolitic materials

Hydride pinning pathway in the hydrogenation of CO2 into formic acid on dimeric tin dioxide

Capture of CO2 and its conversion into organic feedstocks are increasingly needed as society moves towards a renewable energy economy. Here, a hydride‐assisted selective reduction pathway is proposed for the conversion of CO2 to formic acid (FA) over SnO2 monomers and dimers. Our density functional theory calculations infer a strong chemisorption of CO2 on SnO2 clusters forming a carbonate structure, whereas heterolytic cleavage of H2 provides a new pathway for the selective reduction of CO2 to formic acid at low overpotential. Among the two investigated pathways for reduction of CO2 to HCOOH, the hydride pinning pathway is found promising with a unique selectivity for HCOOH. The negatively‐charged hydride forms on the cluster during the dissociation of H2 and facilitates the formation of a formate intermediate, which determines the selectivity for FA over the alternative CO and H2 evolution reaction. It is confirmed that SnO2 clusters exhibit a different catalytic behaviour from their surface equivalents, thus offering promise for future work investigating the reduction of CO2 to FA via a hydride pinning pathway at low overpotential and CO2 capturing

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