SrTiO3(001)(2x1) reconstructions: First-principles calculations of surface energy and atomic structure compared with scanning tunneling microscopy images

(1x1) and (2x1) reconstructions of the (001) SrTiO3 surface were studied using the first-principles full-potential linear muffin-tin orbital method. Surface energies were calculated as a function of TiO2 chemical potential, oxygen partial pressure and temperature. The (1x1) unreconstructed surfaces were found to be energetically stable for many of the conditions considered. Under conditions of very low oxygen partial pressure the (2x1) Ti2O3 reconstruction [Martin R. Castell, Surf. Sci. 505, 1 (2002)] is stable. The question as to why STM images of the (1x1) surfaces have not been obtained was addressed by calculating charge densities for each surface. These suggest that the (2x1) reconstructions would be easier to image than the (1x1) surfaces. The possibility that the presence of oxygen vacancies would destabilise the (1x1) surfaces was also investigated. If the (1x1) surfaces are unstable then there exists the further possibility that the (2x1) DL-TiO2 reconstruction [Natasha Erdman Nature (London) 419, 55 (2002)] is stable in a TiO2-rich environment and for p(O2)>10(-18) atm.

Structural and chemical embrittlement of grain boundaries by impurities: A general theory and first-principles calculations for copper

First-principles calculations of the Sigma 5(310)[001] symmetric tilt grain boundary in Cu with Bi, Na, and Ag substitutional impurities provide evidence that in the phenomenon of Bi embrittlement of Cu grain boundaries electronic effects do not play a major role; on the contrary, the embrittlement is mostly a structural or "size" effect. Na is predicted to be nearly as good an embrittler as Bi, whereas Ag does not embrittle the boundary in agreement with experiment. While we reject the prevailing view that "electronic" effects (i.e., charge transfer) are responsible for embrittlement, we do not exclude the role of chemistry. However, numerical results show a striking equivalence between the alkali metal Na and the semimetal Bi, small differences being accounted for by their contrasting "size" and "softness" (defined here). In order to separate structural and chemical effects unambiguously if not uniquely, we model the embrittlement process by taking the system of grain boundary and free surfaces through a sequence of precisely defined gedanken processes; each of these representing a putative mechanism. We thereby identify three mechanisms of embrittlement by substitutional impurities, two of which survive in the case of embrittlement or cohesion enhancement by interstitials. Two of the three are purely structural and the third contains both structural and chemical elements that by their very nature cannot be further unraveled. We are able to take the systems we study through each of these stages by explicit computer simulations and assess the contribution of each to the net reduction in intergranular cohesion. The conclusion we reach is that embrittlement by both Bi and Na is almost exclusively structural in origin; that is, the embrittlement is a size effect.

An understanding of chemoselective hydrogenation on crotonaldehyde over Pt(111) in the free energy landscape: The microkinetics study based on first-principles calculations

Microkinetic model is developed in the free energy landscape based on density functional theory (DFT) to quantitatively investigate the reaction mechanism of chemoselective partial hydrogenation of crotonaldehyde to crotyl alcohol over Pt(1 1 1) at the temperature of 353 K. Three different methods (mobile, immobile and collision theory models) were carried out to obtain free energy barrier of adsorption/desorption processes. The results from mobile and collision theory models are similar. The calculated TOFs from both models are close to the experiment value. However, for the immobile model, in which the free energy barrier of desorption approaches the energy barrier, the calculated TOF is 2 orders of magnitude lower than the other models. The difficulty of adsorption/ desorption may be overestimated in the immobile model. In addition, detailed analyses show that for the surface hydrogenation elementary steps, the entropy and internal energy effects are small under the reaction condition, while the zero-point-energy (ZPE) correction is significant, especially for the multi-step hydrogenation reaction. The total energy with the ZPE correction approaches to the full free energy calculation for the surface reaction under the reaction condition. (c) 2011 Elsevier B.V. All rights reserved.

Identifying the trend of reactivity for sp2 materials: an electron delocalization model from first principles calculations

The reactivity of sp² carbon materials is studied using the adsorption and dissociation of O₂ on graphene and graphene oxide as model systems. The reactions on the basal plane, zigzag and armchair edges of graphene and graphene oxide with different oxygen-containing groups are calculated using first principles calculations. Two Brønsted-Evans- Polanyi relationships are identified and an electron delocalization model is suggested to understand the general trend of reactivity for sp² carbon materials.

NO Oxidation on Platinum Group Metals Oxides:First Principles Calculations Combined with Microkinetic Analysis

By combining density functional theory calculation and microkinetic analysis, NO oxidation on the platinum group metal oxides (PtO(2), IrO(2), OsO(2)) is investigated, aiming at shedding light on the activities of metal oxides and exploring the activity variations of metal oxides compared to their corresponding metals. A microkinetic model, taking into account the possible low diffusion of surface species on metal oxide surfaces, is proposed for NO oxidation. The resultant turnover frequencies of NO oxidation show that under the typical experimental condition, T = 600 K, p(O2) = 0.1 atm, p(NO) = 3 x 10(-4) atm, p(NO2) = 1.7 x 10(-4) atm; (i) IrO(2)(110) exhibits higher activity than PtO(2)(110) and OsO(2)(110), and (ii) compared to the corresponding metallic Pt, Ir, and Os, the activity of PtO(2) to catalyze NO oxidation is lower, but interestingly IrO(2) and OsO(2) exhibit higher activities. The reasons for the activity differences between the metals and oxides are addressed. Moreover, other possible reaction pathways of NO oxidation on PtO(2)(110), involving O(2) molecule (NO + O(2) -> OONO) and lattice bridge-O(2c), are also found to give low activities. The origin of the Pt catalyst deactivation is also discussed.

Activity and coke formation of nickel and nickel carbide in dry reforming: A deactivation scheme from density functional theory

Dry reforming is a promising reaction to utilise the greenhouse gases CO2 and CH4. Nickel-based catalysts are the most popular catalysts for the reaction, and the coke formation on the catalysts is the main obstacle to the commercialisation of dry reforming. In this study, the whole reaction network of dry reformation on both flat and stepped nickel catalysts (Ni(111) and Ni(211)) as well as nickel carbide (flat: Ni3C(001); stepped: Ni3C(111)) is investigated using density functional theory calculations. The overall reaction energy profiles in the free energy landscape are obtained, and kinetic analyses are utilised to evaluate the activity of the four surfaces. By careful examination of our results, we find the following regarding the activity: (i) flat surfaces are more active than stepped surfaces for the dry reforming and (ii) metallic nickel catalysts are more active than those of nickel carbide, and therefore, the phase transformation from nickel to nickel carbide will reduce the activity. With respect to the coke formation, the following is found: (i) the coke formation probability can be measured by the rate ratio of CH oxidation pathway to C oxidation pathway (r(CH)/r(C)) and the barrier of CO dissociation, (ii) on Ni(111), the coke is unlikely to form, and (iii) the coke formations on the stepped surfaces of both nickel and nickel carbide can readily occur. A deactivation scheme, using which experimental results can be rationalised, is proposed. 

First-principles study of ferroelectricity and isotope effects in H-bonded KH2PO4 crystals

By means of extensive first-principles calculations we studied the ferroelectric phase transition and the associated isotope effect in KH2PO4 (KDP). Our calculations revealed that the spontaneous polarization of the ferroelectric phase is due to electronic charge redistributions and ionic displacements which are a consequence of proton ordering, and not vice versa. The experimentally observed double-peaked proton distribution in the paraelectric phase cannot be explained by a dynamics of only protons. This requires, instead, collective displacements within clusters that include also the heavier ions. These tunneling clusters can explain the recent evidence of tunneling obtained from Compton scattering measurements. The sole effect of mass change upon deuteration is not sufficient to explain the huge isotope effect. Instead, we find that structural modifications deeply connected with the chemistry of the H bonds produce a feedback effect on tunneling that strongly enhances the phenomenon. The resulting influence of the geometric changes on the isotope effect agrees with experimental data from neutron scattering. Calculations under pressure allowed us to analyze the issue of universality in the disappearance of ferroelectricity upon compression. Compressing DKDP so that the distance between the two peaks in the deuteron distribution is the same as for protons in KDP, corresponds to a modification of the underlying double-well potential, which becomes 23 meV shallower. This energy difference is what is required to modify the O-O distance in such a way as to have the same distribution for protons and deuterons. At the high pressures required experimentally, the above feedback mechanism is crucial to explain the magnitude of the geometrical effect.

Stability of Sr adatom model structures for SrTiO3(001) surface reconstructions

We report results of first-principles calculations on the thermodynamic stability of different Sr adatom structures that have been proposed to explain some of the observed reconstructions of the (001) surface of strontium titanate (Kubo and Nozoye 2003 Surf Sci. 542 177). From surface free energy calculations, a phase diagram is constructed indicating the range of conditions over which each structure is most stable. These results are compared with Kubo and Nozoye's experimental observations. It is concluded that low Sr adatom coverage structures can only be explained if the surface is far from equilibrium. Intermediate coverage structures are stable only if the surface is in or very nearly in equilibrium with the strontium oxide.

A stabilization mechanism of zirconia based on oxygen vacancies only

The microscopic mechanism leading to stabilization of cubic and tetragonal forms of zirconia (ZrO2) is analyzed by means of a self-consistent tight-binding model. Using this model, energies and structures of zirconia containing different vacancy concentrations are calculated, equivalent in concentration to the charge compensating vacancies associated with dissolved yttria (Y2O3) in the tetragonal and cubic phase fields (3.2 and 14.4% mol, respectively). The model is shown to predict the large relaxations around an oxygen vacancy, and the clustering of vacancies along the 111 directions, in good agreement with experiments and first principles calculations. The vacancies alone are shown to explain the stabilization of cubic zirconia, and the mechanism is analyzed. (C) 2002 Acta Materialia Inc. Published by Elsevier Science Ltd. All rights reserved.

Relevance of heterometallic binding energy for metal underpotential deposition

We present first-principles calculations for a number of metals adsorbed on several different metallic substrates. Some of these systems are very relevant in electrochemistry, especially in the field of underpotential deposition phenomena. The present studies reveal the existence of a relationship between the excess binding energy and the surface energy difference between substrate and adsorbate. Comparisons with experimental underpotential shifts show that excess binding energies are systematically underestimated. By analyzing experimental information on different systems, we conclude that this discrepancy between our vacuum calculations and experiments carried out in an electrolytic solution is likely to be due to anion adsorption and/or solvent effects.

Broken symmetry and pseudogaps in ropes of carbon nanotubes

Since the discovery of carbon nanotubes, it has been speculated that these materials should behave like nanoscale wires with unusual electronic properties and exceptional strength. Recently, 'ropes' of close-packed single-wall nanotubes have been synthesized in high yield. The tubes in these ropes are mainly of the (10,10) type3, which is predicted to be metallic. Experiments on individual nanotubes and ropes indicate that these systems indeed have transport properties that qualify them to be viewed as nanoscale quantum wires at low temperature. It has been expected that the close-packing of individual nanotubes into ropes does not change their electronic properties significantly. Here, however, we present first-principles calculations which show that a broken symmetry of the (10,10) tube caused by interactions between tubes in a rope induces a pseudogap of about 0.1 eV at the Fermi level. This pseudogap strongly modifies many of the fundamental electronic properties: we predict a semimetal-like temperature dependence of the electrical conductivity and a finite gap in the infrared absorption spectrum. The existence of both electron and hole charge carriers will lead to qualitatively different thermopower and Hall-effect behaviours from those expected for a normal metal.

Are current-induced forces conservative?

The expression for the force on an ion in the presence of current can be derived from first principles without any assumption about its conservative character. However, energy functionals have been constructed that indicate that this force can be written as the derivative of a potential. On the other hand, there exist specific arguments that strongly suggest the contrary. We propose physical mechanisms that invalidate such arguments and demonstrate their existence with first-principles calculations. While our results do not constitute a formal resolution to the fundamental question of whether current-induced forces are conservative, they represent a substantial step forward in this direction.

Relative energetics and structural properties of zirconia using a self-consistent tight-binding model

We describe an empirical, self-consistent, orthogonal tight-binding model for zirconia, which allows for the polarizability of the anions at dipole and quadrupole levels and for crystal field splitting of the cation d orbitals, This is achieved by mixing the orbitals of different symmetry on a site with coupling coefficients driven by the Coulomb potentials up to octapole level. The additional forces on atoms due to the self-consistency and polarizabilities are exactly obtained by straightforward electrostatics, by analogy with the Hellmann-Feynman theorem as applied in first-principles calculations. The model correctly orders the zero temperature energies of all zirconia polymorphs. The Zr-O matrix elements of the Hamiltonian, which measure covalency, make a greater contribution than the polarizability to the energy differences between phases. Results for elastic constants of the cubic and tetragonal phases and phonon frequencies of the cubic phase are also presented and compared with some experimental data and first-principles calculations. We suggest that the model will be useful for studying finite temperature effects by means of molecular dynamics.

Electronic structures and bonding properties of chlorine-treated nitrogenated carbon nanotubes: X-ray absorption and scanning photoelectron microscopy studies

The electronic and bonding properties of nitrogenated carbon nanotubes (N-CNTs) exposed to chlorine plasma were investigated using C and N K-edge x-ray absorption near-edge structure (XANES) and scanning photoelectron microscopy (SPEM). The C and N K-edge XANES spectra of chlorine-treated N-CNTs consistently reveal the formation of pyridinelike N-CNTs by the observation of 1s ->pi(*)(e(2u)) antibonding and 1s ->pi(*)(b(2g)) bonding states. The valence-band photoemission spectra obtained from SPEM images indicate that chlorination of the nanotubes enhances the C-N bonding. First-principles calculations of the partial densities of states in conjunction with C K-edge XANES data identify the presence of C-Cl bonding in chlorine treated N-CNTs. (C) 2007 American Institute of Physics.

Origin of Low CO2 Selectivity on Platinum in the Direct Ethanol Fuel Cell

Calculated answer: First-principles calculations have been applied to calculate the energy barrier for the key step in CO formation on a Pt surface (see picture; Pt blue, Pt atoms on step edge yellow) to understand the low CO2 selectivity in the direct ethanol fuel cell. The presence of surface oxidant species such as O (brown bar) and OH (red bar) led to an increase of the energy barrier and thus an inhibition of the key step. © 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim