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

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

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

Stability of single gold atoms on defective and doped diamond surfaces

Polycrystalline boron-doped diamond (BDD) is widely used as a working electrode material in electrochemistry, and its properties, such as its stability, make it an appealing support material for nanostructures in electrocatalytic applications. Recent experiments have shown that electrodeposition can lead to the creation of stable small nanoclusters and even single gold adatoms on the BDD surfaces. We investigate the adsorption energy and kinetic stability of single gold atoms adsorbed onto an atomistic model of BDD surfaces by using density functional theory. The surface model is constructed using hybrid quantum mechanics/molecular mechanics embedding techniques and is based on an oxygen-terminated diamond (110) surface. We use the hybrid quantum mechanics/molecular mechanics method to assess the ability of different density functional approximations to predict the adsorption structure, energy, and barrier for diffusion on pristine and defective surfaces. We find that surface defects (vacancies and surface dopants) strongly anchor adatoms on vacancy sites. We further investigated the thermal stability of gold adatoms, which reveals high barriers associated with lateral diffusion away from the vacancy site. The result provides an explanation for the high stability of experimentally imaged single gold adatoms on BDD and a starting point to investigate the early stages of nucleation during metal surface deposition

Dataset to support : FHI-aims benchmark for BEEF-vdW, vdW-DF2 and mBEEFvdW

We have performed a benchmark of the BEEF-vdW, mBEEF-vdW and vdW-DF2 implementations in the all-electron electronic structure code in FHIaims against the S22 benchmark set

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

Highlights from faraday discussion on designing nanoparticle systems for catalysis, London, UK, May 2018

The 2018 Faraday Discussion on “Designing Nanoparticle Systems for Catalysis” brought together leading scientists to discuss the current state-of-the-art in the fields of computational chemistry, characterization techniques, and nanomaterial synthesis, and to debate the challenges and opportunities going forward for rational catalyst design. The meeting was a vivid discussion of how the communities accummulate knowledge and on how innovativeness can be combined to have a stronger scientific impact. In the following, we provide an overview of the meeting structure, including plenaries, papers, discussion points and breakout sessions, and we hope to show, to the wider scientific community, that there is great value in continued international discussion and scientific collaboration in these fields

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

Dehydrogenation and dehydration of formic acid over orthorhombic molybdenum carbide

The dehydrogenation and dehydration of formic acid is investigated on the β-Mo2C (100) catalyst surface using time independent density functional theory. The energetics of the two mechanisms are calculated, and the thermochemistry and kinetics are discussed using the transition state theory. Subsequently, microkinetic modelling of the system is conducted, considering the batch reactor model. The potential energy landscape of the reaction shows a thermodynamically favourable cleavage of Hsingle bondCOOH to form CO; however, the kinetics show that the dehydrogenation mechanism is faster and CO2 is continuously formed. The effect of HCOOH adsorption on the surface is also analysed, in a temperature-programmed desorption, with the conversion proceeding at under 350 K and desorption of CO2 is observed with a selectivity of about 100 %, in line with the experimental reports

Hydrodeoxygenation of guaiacol over orthorhombic molybdenum carbide: a DFT and microkinetic study

The hydrodeoxygenation of guaiacol is modelled over a (100) β-Mo2C surface using density functional theory and microkinetic simulations. The thermochemistry of the process shows that the demethoxylation of the guaiacol, to form phenol, will be the initial steps, with a reaction energy of 29 kJ/mol (i.e. endothermic) and a highest activation barrier of 112 kJ/mol. Subsequently, the dehydroxylation of the phenol, which has a rate-determining activation barrier of 145 kJ/mol, will lead to the formation of benzene, with an overall reaction energy for conversion from guaiacol of -91 kJ/mol (i.e. exothermic). The sp2 and sp hybridized carbon atoms of the molecular functional groups are found to dissociate on the surface with minimum energy barriers, while the hydrogenation of the adsorbed molecules requires higher energy. The microkinetic modelling, which is performed considering typical reaction conditions of 500 to 700 K, and a partial pressure ratio H2:guaiacol of 1, shows quick formation and accumulation of phenol on the surface with increasing temperature, although high temperatures mitigate the guaiacol adsorption step. Based on simulated temperature programmed desorption (TPD), maximum conversion of guaiacol can be expected at 70% surface coverage of this species