Adsorption of organic molecules at the TiO2(110) surface: the effect of van der Waals interactions

Understanding the interaction of organic molecules with TiO2 surfaces is important for a wide range of technological applications. While density functional theory (DFT) calculations can provide valuable insight about these interactions, traditional DFT approaches with local exchange-correlation functionals suffer from a poor description of non-bonding van der Waals (vdW) interactions. We examine here the contribution of vdW forces to the interaction of small organic molecules (methane, methanol, formic acid and glycine) with the TiO2 (110) surface, based on DFT calculations with the optB88-vdW functional. The adsorption geometries and energies at different configurations were also obtained in the standard generalized gradient approximation (GGA-PBE) for comparison. We find that the optB88-vdW consistently gives shorter surface adsorbate-to-surface distances and slightly stronger interactions than PBE for the weak (physisorbed) modes of adsorption. In the case of strongly adsorbed (chemisorbed) molecules both functionals give similar results for the adsorption geometries, and also similar values of the relative energies between different chemisorption modes for each molecule. In particular both functionals predict that dissociative adsorption is more favourable than molecular adsorption for methanol, formic acid and glycine, in general agreement with experiment. The dissociation energies obtained from both functionals are also very similar, indicating that vdW interactions do not affect the thermodynamics of surface deprotonation. However, the optB88-vdW always predicts stronger adsorption than PBE. The comparison of the methanol adsorption energies with values obtained from a Redhead analysis of temperature programmed desorption data suggests that optB88-vdW significantly overestimates the adsorption strength, although we warn about the uncertainties involved in such comparisons.

A DFT study of the structures, stabilities and redox behaviour of the major surfaces of magnetite Fe3O4

The renewed interest in magnetite (Fe3O4) as a major phase in different types of catalysts has led us to study the oxidation–reduction behaviour of its most prominent surfaces. We have employed computer modelling techniques based on the density functional theory to calculate the geometries and surface free energies of a number of surfaces at different compositions, including the stoichiometric plane, and those with a deficiency or excess of oxygen atoms. The most stable surfaces are the (001) and (111), leading to a cubic Fe3O4 crystal morphology with truncated corners under equilibrium conditions. The scanning tunnelling microscopy images of the different terminations of the (001) and (111) stoichiometric surfaces were calculated and compared with previous reports. Under reducing conditions, the creation of oxygen vacancies in the surface leads to the formation of reduced Fe species in the surface in the vicinity of the vacant oxygen. The (001) surface is slightly more prone to reduction than the (111), due to the higher stabilisation upon relaxation of the atoms around the oxygen vacancy, but molecular oxygen adsorbs preferentially at the (111) surface. In both oxidized surfaces, the oxygen atoms are located on bridge positions between two surface iron atoms, from which they attract electron density. The oxidised state is thermodynamically favourable with respect to the stoichiometric surfaces under ambient conditions, although not under the conditions when bulk Fe3O4 is thermodynamically stable with respect to Fe2O3. This finding is important in the interpretation of the catalytic properties of Fe3O4 due to the presence of oxidised species under experimental conditions.

The interaction of iron pyrite with oxygen, nitrogen and nitrogen oxides: a first-principles study

Sulphide materials, in particular MoS2, have recently received great attention from the surface science community due to their extraordinary catalytic properties. Interestingly, the chemical activity of iron pyrite (FeS2) (the most common sulphide mineral on Earth), and in particular its potential for catalytic applications, has not been investigated so thoroughly. In this study, we use density functional theory (DFT) to investigate the surface interactions of fundamental atmospheric components such as oxygen and nitrogen, and we have explored the adsorption and dissociation of nitrogen monoxide (NO) and nitrogen dioxide (NO2) on the FeS2(100) surface. Our results show that both those environmentally important NOx species chemisorb on the surface Fe sites, while the S sites are basically unreactive for all the molecular species considered in this study and even prevent NO2 adsorption onto one of the non-equivalent Fe–Fe bridge sites of the (1 1)–FeS2(100) surface. From the calculated high barrier for NO and NO2 direct dissociation on this surface, we can deduce that both nitrogen oxides species are adsorbed molecularly on pyrite surfaces.

Electronic structure of Pd multi-layers on Re(0001): the role of charge transfer

Understanding the origin of the properties of metal-supported metal thin films is important for the rational design of bimetallic catalysts and other applications, but it is generally difficult to separate effects related to strain from those arising from interface interactions. Here we use density functional (DFT) theory to examine the structure and electronic behavior of few-layer palladium films on the rhenium (0001) surface, where there is negligible interfacial strain and therefore other effects can be isolated. Our DFT calculations predict stacking sequences and interlayer separations in excellent agreement with quantitative low-energy electron diffraction experiments. By theoretically simulating the Pd core-level X-ray photoemission spectra (XPS) of the films, we are able to interpret and assign the basic features of both low-resolution and high-resolution XPS measurements. The core levels at the interface shift to more negative energies, rigidly following the shifts in the same direction of the valence d-band center. We demonstrate that the valence band shift at the interface is caused by charge transfer from Re to Pd, which occurs mainly to valence states of hybridized s-p character rather than to the Pd d-band. Since the d-band filling is roughly constant, there is a correlation between the d-band center shift and its bandwidth. The resulting effect of this charge transfer on the valence d-band is thus analogous to the application of a lateral compressive strain on the adlayers. Our analysis suggests that charge transfer should be considered when describing the origin of core and valence band shifts in other metal / metal adlayer systems.

The atomic structure of low-index surfaces of the intermetallic compound InPd

The intermetallic compound InPd (CsCl type of crystal structure with a broad compositional range) is considered as a candidate catalyst for the steam reforming of methanol. Single crystals of this phase have been grown to study the structure of its three low-index surfaces under ultra-high vacuum conditions, using low energy electron diffraction (LEED), X-ray photoemission spectroscopy (XPS), and scanning tunneling microscopy (STM). During surface preparation, preferential sputtering leads to a depletion of In within the top few layers for all three surfaces. The near-surface regions remain slightly Pd-rich until annealing to ∼580 K. A transition occurs between 580 and 660 K where In segregates towards the surface and the near-surface regions become slightly In-rich above ∼660 K. This transition is accompanied by a sharpening of LEED patterns and formation of flat step-terrace morphology, as observed by STM. Several superstructures have been identified for the different surfaces associated with this process. Annealing to higher temperatures (≥750 K) leads to faceting via thermal etching as shown for the (110) surface, with a bulk In composition close to the In-rich limit of the existence domain of the cubic phase. The Pd-rich InPd(111) is found to be consistent with a Pd-terminated bulk truncation model as shown by dynamical LEED analysis while, after annealing at higher temperature, the In-rich InPd(111) is consistent with an In-terminated bulk truncation, in agreement with density functional theory (DFT) calculations of the relative surface energies. More complex surface structures are observed for the (100) surface. Additionally, individual grains of a polycrystalline sample are characterized by micro-spot XPS and LEED as well as low-energy electron microscopy. Results from both individual grains and “global” measurements are interpreted based on comparison to our single crystals findings, DFT calculations and previous literature.

Asymmetric oxidation of vinyl- and ethynyl terthiophene ligands in triruthenium complexes

A series of ruthenium(II) complexes [{RuCl(CO)(PMe3)3(–CHvCH–)}nX], 1a–1c (1a: n = 3, X = 3,3’’- dimethyl-2,2’:3’,2’’-terthiophene; 1b: n = 2, X = 2,2’-bithiophene; 1c: n = 2, X = 2,3-bis(3-methylthiophen- 2-yl)benzothiophene) and [{Cp*(dppe)2Ru(–CuC–)}3X], 1d (X = 3,3’’-dimethyl-2,2’:3’,2’’- terthiophene), were prepared and characterized by 1H, 13C and 31P NMR. Their redox, spectroscopic and bonding properties were studied with a range of spectro-electrochemical methods in combination with density functional theory calculations. The first two anodic steps observed for 1a and 1d are largely localized on the lateral frameworks of the molecular triangle, the direct conjugation between them being precluded due to the photostable open form of the dithienyl ethene moiety. The third anodic step is then mainly localized on the centerpiece of the triangular structure, affecting both bithiophene laterals. The experimental IR and UV-vis-NIR spectroelectrochemical data and, largely, also DFT calculations account for this explanation, being further supported by direct comparison with the anodic behavior of reference diruthenium complexes 1b and 1c.

Quantitative structure of an acetate dye molecule analogue at the TiO2- acetic acid interface

The positions of atoms in and around acetate molecules at the rutile TiO2(110) interface with 0.1 M acetic acid have been determined with a precision of ±0.05 Å. Acetate is used as a surrogate for the carboxylate groups typically employed to anchor monocarboxylate dye molecules to TiO2 in dye-sensitised solar cells (DSSC). Structural analysis reveals small domains of ordered (2 x 1) acetate molecules, with substrate atoms closer to their bulk terminated positions compared to the clean UHV surface. Acetate is found in a bidentate bridge position, binding through both oxygen atoms to two five-fold titanium atoms such that the molecular plane is along the [001] azimuth. Density functional theory calculations provide adsorption geometries in excellent agreement with experiment. The availability of these structural data will improve the accuracy of charge transport models for DSSC.