Phase separation and surface segregation in ceria-zirconia solid solutions

Using a combination of density functional theory calculations and statistical mechanics, we show that a wide range of intermediate compositions of ceria – zirconia solid solutions are thermodynamically metastable with respect to phase separation into Ce-rich and Zr-rich oxides. We estimate that the maximum equilibrium concentration of Zr in CeO2 at 1373 K is ~2%, and therefore equilibrated samples with higher Zr content are expected to exhibit heterogeneity at the atomic scale. We also demonstrate that in the vicinity of the (111) surface, cation redistribution at high temperatures will occur with significant Ce enrichment of the surface, which we attribute to the more covalent character of Zr-O bonds compared to Ce-O bonds. Although the kinetic barriers for cation diffusion normally prevent the decomposition/segregation of ceria-zirconia solid solutions in typical catalytic applications, the separation behaviour described here can be expected to occur in modern three-way catalytic converters, where very high temperatures are reached.

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.

Spin polarization, orbital occupation and band gap opening in vanadium dioxide: the effect of screened Hartree–Fock exchange

The metal–insulator transition of VO2 so far has evaded an accurate description by density functional theory. The screened hybrid functional of Heyd, Scuseria and Ernzerhof leads to reasonable solutions for both the low-temperature monoclinic and high-temperature rutile phases only if spin polarization is excluded from the calculations. We explore whether a satisfactory agreement with experiment can be achieved by tuning the fraction of Hartree Fock exchange (a) in the density functional. It is found that two branches of locally stable solutions exist for the rutile phase for 12:5% 6 a 6 20%. One is metallic and has the correct stability as compared to the monoclinic phase, the other is insulating with lower energy than the metallic branch. We discuss these observations based on the V 3d orbital occupations and conclude that a ¼ 10% is the best possible choice for spin-polarized VO2 calculations.

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.

First-principles study of the inversion thermodynamics and electronic structure of FeM2X4 (thio)spinels (M = Cr, Mn, Co, Ni; X = O, S)

FeM2X4 spinels, where M is a transition metal and X is oxygen or sulfur, are candidate materials for spin filters, one of the key devices in spintronics. We present here a computational study of the inversion thermodynamics and the electronic structure of these (thio)spinels for M = Cr, Mn, Co, Ni, using calculations based on the density functional theory with on-site Hubbard corrections (DFT+U). The analysis of the configurational free energies shows that different behaviour is expected for the equilibrium cation distributions in these structures: FeCr2X4 and FeMn2S4 are fully normal, FeNi2X4 and FeCo2S4 are intermediate, and FeCo2O4 and FeMn2O4 are fully inverted. We have analyzed the role played by the size of the ions and by the crystal field stabilization effects in determining the equilibrium inversion degree. We also discuss how the electronic and magnetic structure of these spinels is modified by the degree of inversion, assuming that this could be varied from the equilibrium value. We have obtained electronic densities of states for the completely normal and completely inverse cation distribution of each compound. FeCr2X4, FeMn2X4, FeCo2O4 and FeNi2O4 are half-metals in the ferrimagnetic state when Fe is in tetrahedral positions. When M is filling the tetrahedral positions, the Cr-containing compounds and FeMn2O4 are half-metallic systems, while the Co and Ni spinels are insulators. The Co and Ni sulfide counterparts are metallic for any inversion degree together with the inverse FeMn2S4. Our calculations suggest that the spin filtering properties of the FeM2X4 (thio)spinels could be modified via the control of the cation distribution through variations in the synthesis conditions.

The importance of anisotropic Coulomb interaction in LaMnO3

In low-temperature anti-ferromagnetic LaMnO3, strong and localized electronic interactions among Mn 3d electrons prevent a satisfactory description from standard local density and generalized gradient approximations in density functional theory calculations. Here we show that the strong on-site electronic interactions are described well only by using direct and exchange corrections to the intra-orbital Coulomb potential. Only DFT+U calculations with explicit exchange corrections produce a balanced picture of electronic, magnetic and structural observables in agreement with experiment. To understand the reason, a rewriting of the functional form of the +U corrections is presented that leads to a more physical and transparent understanding of the effect of these correction terms. The approach highlights the importance of Hund’s coupling (intra-orbital exchange) in providing anisotropy across the occupation and energy eigenvalues of the Mn d states. This intra-orbital exchange is the key to fully activating the Jahn-Teller distortion, reproducing the experimental band gap and stabilizing the correct magnetic ground state in LaMnO3. The best parameter values for LaMnO3 within the DFT(PBEsol)+U framework are determined to be U = 8 eV and J = 1.9 eV.

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.

Characterization and spectroscopic analysis of phenol-ethanol hydrogen bonded clusters

Different hydrogen bonded clusters involving phenol and ethanol are studied theoretically using MP2/aug-cc-pVDZ. Nine different 1: 1 clusters are obtained and analyzed according to their stability and spectroscopic properties. Different isomeric forms of ethanol are considered. Attention is also devoted to the spectral shift of the characteristic pi -> pi* transition of phenol. Using TDHF, CIS, CIS(D) and TDB3LYP in aug-cc-pVDZ basis set, all results agree that a red shift is obtained when phenol is the hydrogen donor and a blue shift is obtained in the opposite case. These results are used to rationalize the red shift observed for phenol in liquid ethanol. (C) 2010 Elsevier B.V. All rights reserved.

Many-body energy decomposition of hydrogen-bonded glycine clusters in gas-phase

A detailed analysis of the many-body contribution to the interaction energies of the gas-phase hydrogen-bonded glycine clusters, (Gly)(N), N = 1-4 is presented. The energetics of the hydrogen-bonded dimer, trimer and tetramer complexes have been analyzed using density-functional theory. The magnitude of the two-through four-body energy terms have been calculated and compared. The relaxation energy and the two-body energy terms are the principal contributors to the total binding energy. Four-body contribution is negligible. However, the three-body contribution is found to be sizable and the formation of the cyclic glycine trimer presents geometric strains that make it less favorable. (C) 2010 Elsevier B.V. All rights reserved.

NMR Chemical Shielding and Spin-Spin Coupling Constants of Liquid NH(3): A Systematic Investigation using the Sequential QM/MM Method

The NMR spin coupling parameters, (1)J(N,H) and (2)J(H,H), and the chemical shielding, sigma((15)N), of liquid ammonia are studied from a combined and sequential QM/MM methodology. Monte Carlo simulations are performed to generate statistically uncorrelated configurations that are submitted to density functional theory calculations. Two different Lennard-Jones potentials are used in the liquid simulations. Electronic polarization is included in these two potentials via an iterative procedure with and without geometry relaxation, and the influence on the calculated properties are analyzed. B3LYP/aug-cc-pVTZ-J calculations were used to compute the V(N,H) constants in the interval of -67.8 to -63.9 Hz, depending on the theoretical model used. These can be compared with the experimental results of -61.6 Hz. For the (2)J(H,H) coupling the theoretical results vary between -10.6 to -13.01 Hz. The indirect experimental result derived from partially deuterated liquid is -11.1 Hz. Inclusion of explicit hydrogen bonded molecules gives a small but important contribution. The vapor-to-liquid shifts are also considered. This shift is calculated to be negligible for (1)J(N,H) in agreement with experiment. This is rationalized as a cancellation of the geometry relaxation and pure solvent effects. For the chemical shielding, U(15 N) Calculations at the B3LYP/aug-pcS-3 show that the vapor-to-liquid chemical shift requires the explicit use of solvent molecules. Considering only one ammonia molecule in an electrostatic embedding gives a wrong sign for the chemical shift that is corrected only with the use of explicit additional molecules. The best result calculated for the vapor to liquid chemical shift Delta sigma((15)N) is -25.2 ppm, in good agreement with the experimental value of -22.6 ppm.

Polarization and Spectral Shift of Benzophenone in Supercritical Water

Monte Carlo simulation and quantum mechanics calculations based on the INDO/CIS and TD-DFT methods were utilized to study the solvatochromic shift of benzophenone when changing the environment from normal water to supercritical (P = 340.2 atm and T = 673 K) condition. Solute polarization increases the dipole moment of benzophenone, compared to gas phase, by 88 and 35% in normal and supercritical conditions, giving the in-solvent dipole value of 5.8 and 4.2 D, respectively. The average number of solute-solvent hydrogen bonds was analyzed, and a large decrease of 2.3 in normal water to only 0.8 in the supercritical environment was found. By using these polarized models of benzophenone in the two different conditions of water, we performed MC simulations to generate statistically uncorrelated configurations of the solute surrounded by the solvent molecules and subsequent quantum mechanics calculations on these configurations. When changing from normal to supercritical water environment, INDO/CIS calculations explicitly considering all valence electrons of the 235 solvent water molecules resulted in a solvatochromic shift of 1425 cm(-1) for the most intense transition of benzophenone, that is, slightly underestimated in comparison with the experimentally inferred result of 1700 cm(-1). TD-B3LYP/6-311+G(2d,p) calculations on the same configurations but with benzophenone electrostatically embedded in the 320 water molecules resulted in a solvatochromic shift of 1715 cm(-1) for this transition, in very good agreement with the experimental result. When using the unpolarized model of the benzophenone, this calculated solvatochromic shift was only 640 cm(-1). Additional calculations were also made by using BHandHLYP/6-311+G(2d,p) to analyze the effect of the asymptotic decay of the exchange functional. This study indicates that, contrary to the general expectation, there is a sizable solute polarization even in the low-density regime of supercritical condition and that the inclusion of this polarization is important for a reliable description of the spectral shifts considered here.

Conformational behavior of different possible ways of oligoglycine formation in a solvent-free environment

The possible ways for glycine oligopeptide formation in gas phase, both in the extended P-strand like conformation and folded 2(7)-ribbon like conformations are analyzed using quantum chemical calculations. We focus on the sequential formation of peptide bond through upgradation of the immediate lower order molecule and observe the consequences in other related processes like oligoglycine formation through simultaneous peptide linkage of n glycine monomers and interchange of molecular conformation through peptide linkage. A comparison is made between the structures and binding energies obtained for both conformers. All binding energies are increased by the zero-point energy contribution. The role of electron correlation effects is briefly analyzed. The folded 2(7)-ribbon-like conformations in vacuo are found to be more stable in comparison to the extended structure. (c) 2007 Elsevier B.V. All rights reserved.