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.

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.

Strong second harmonic generation in SiC, ZnO, GaN two-dimensional hexagonal crystals from first-principles many-body calculations

The second harmonic generation (SHG) intensity spectrum of SiC, ZnO, GaN two-dimensional hexagonal crystals is calculated by using a real-time first-principles approach based on Green's function theory [Attaccalite et al., Phys. Rev. B: Condens. Matter Mater. Phys. 2013 88, 235113]. This approach allows one to go beyond the independent particle description used in standard first-principles nonlinear optics calculations by including quasiparticle corrections (by means of the GW approximation), crystal local field effects and excitonic effects. Our results show that the SHG spectra obtained using the latter approach differ significantly from their independent particle counterparts. In particular they show strong excitonic resonances at which the SHG intensity is about two times stronger than within the independent particle approximation. All the systems studied (whose stabilities have been predicted theoretically) are transparent and at the same time exhibit a remarkable SHG intensity in the range of frequencies at which Ti:sapphire and Nd:YAG lasers operate; thus they can be of interest for nanoscale nonlinear frequency conversion devices. Specifically the SHG intensity at 800 nm (1.55 eV) ranges from about 40-80 pm V(-1) in ZnO and GaN to 0.6 nm V(-1) in SiC. The latter value in particular is 1 order of magnitude larger than values in standard nonlinear crystals.

Catalytic Role of Metal Oxides in Gold-Based Catalysts: A First Principles Study of CO Oxidation on TiO2 Supported Au

CO oxidation on TiO2 supported Au has been studied using density functional theory calculations. Important catalytic roles of the oxide have been identified: (i) CO oxidation occurs at the interface between Au and the oxide with a very small barrier; and (ii) O-2 adsorption at the interface is the key step in the reaction. The physical origin of the oxide promotion effect has been further investigated: The oxide enhances electron transfer from the Au to the antibonding states of O-2, giving rise to (i) strong ionic bonding between the adsorbed O-2, Au, and the Ti cation; and (ii) a significant activation of O-2 towards CO oxidation.

First-principles study of oxygenates on Co surfaces in Fischer-Tropsch synthesis

Extensive density function theory calculations are performed to study the mechanism of the formation of aldehyde and alcohol on Co surfaces in Fischer-Tropsch synthesis, a challenging issue in heterogeneous catalysis. Three possible pathways for the production of formaldehyde and methanol on flat and stepped Co(0001) surfaces are investigated: (i) CO + 4H -> CHO + 3H -> CH2O + 2H -> CH3O + H -> CH3OH; (ii) CO + 4H -> COH + 3H -> CHOH + 2H -> CH2OH + H -> CH3OH; and (iii) the coupling reactions of CH2 + O -> CH2O and CH3 + OH -> CH3OH. It is found that these pathways are generally favored at step sites, and the preferred mechanism is pathway (i) via CHO. Furthermore, the three traditional chain growth mechanisms in Fischer-Tropsch synthesis are semi quantitatively compared and discussed. Our results suggest that the two mechanisms involving oxygenate intermediates (the CO-insertion and hydroxycarbene mechanisms) are less important than the carbene mechanism in the production of long chain hydrocarbons. However, the CO-insertion mechanism may be responsible for the production of long-chain oxygenates.

A first principles study of methanol decomposition on Pd(111): Mechanisms for O-H bond scission and C-O bond scission

There is some dispute as to whether methanol decomposition occurs by O-H bond scission or C-O bond scission. By carrying out density functional theory calculations, we investigate both scenario of the reaction pathways of methanol decomposition on a Pd(111) surface. It is shown that the O-H bond scission pathway is much more energetically favorable than the C-O bond scission pathway. The high reaction barrier in the latter case is found to be due to the poor bonding abilities of CH3 and OH with the surface at the reaction sites. (C) 2001 American Institute of Physics.

Stepwise addition reactions in ammonia synthesis: A first principles study

Catalytic ammonia synthesis is believed to proceed via dissociation of N-2 and H-2 with subsequent stepwise addition reactions from an adsorbed nitrogen atom to NH3. The first step, N-2 dissociation, has been thoroughly studied. However, little is known about the microscopic details of the stepwise addition reactions. To shed light on these stepwise addition reactions, density functional theory calculations with the generalized gradient approximation are employed to investigate NHx (x=1,3) formation on Ru(0001). Transition states and reaction barriers are determined in each elementary step. It is found that the reaction barriers for stepwise addition reactions are rather high, for example, the barrier for NH hydrogenation is calculated to be 1.28 eV, which is comparable with that of N-2 dissociation. In addition, one of the stepwise addition reactions on a stepped surface is also considered. The reaction barrier is found to be much higher than that of N-2 dissociation on the same stepped surface, which indicates the importance of stepwise addition reactions in ammonia synthesis. (C) 2001 American Institute of Physics.

Electronic stopping power of H and He in Al and LiF from first principles

Non-linearities in the electronic stopping power of light projectiles in bulk Al and LiF are addressed from first principles using time-evolving time-dependent density functional theory. In the case of Al, the agreement of the calculations with experiments for H and He projectiles is fair, but a recently observed transition for He from one value of the electronic friction coefficient to a higher value at v ~ 0.3 a.u. is not reproduced by the calculations. For LiF, better accuracy is obtained as compared with previously published simulations, albeit the threshold remains overestimated. © 2013 Elsevier B.V.

A first principles study of CH3 dehydrogenation on Ni(111)

Density functional theory with gradient corrections and spin polarization has been used to study the dehydrogenation of CH3 on Ni(111), a crucial step in many important catalytic reactions. The reaction, CH3(ads)--> CH2(ads)+H-(ads), is about 0.5 eV endothermic with an activation energy of more than 1 eV. The overall reaction pathway is rather intriguing. The C moiety translates from a hcp to a fcc site during the course of the reaction. The transition state of the reaction has been identified. The CH3 species is highly distorted, and both C and the active H are centered nearly on top of a row of Ni atoms with a long C-H bond length of 1.80 Angstrom. The local density of states coupled with examination of the real space distribution of individual quantum states has been used to analyze the reaction pathway. (C) 2000 American Institute of Physics. [S0021-9606(00)30218-5].

First-Principles Study of the Electronic and Magnetic Properties of Defects in Carbon Nanostructures

Understanding the magnetic properties of graphenic nanostructures is instrumental in future spintronics applications. These magnetic properties are known to depend crucially on the presence of defects. Here we review our recent theoretical studies using density functional calculations on two types of defects in carbon nanostructures: Substitutional doping with transition metals, and sp$^3$-type defects created by covalent functionalization with organic and inorganic molecules. We focus on such defects because they can be used to create and control magnetism in graphene-based materials. Our main results are summarized as follows: i)Substitutional metal impurities are fully understood using a model based on the hybridization between the $d$ states of the metal atom and the defect levels associated with an unreconstructed D$_{3h}$ carbon vacancy. We identify three different regimes, associated with the occupation of distinct hybridization levels, which determine the magnetic properties obtained with this type of doping; ii) A spin moment of 1.0 $\mu_B$ is always induced by chemical functionalization when a molecule chemisorbs on a graphene layer via a single C-C (or other weakly polar) covalent bond. The magnetic coupling between adsorbates shows a key dependence on the sublattice adsorption site. This effect is similar to that of H adsorption, however, with universal character; iii) The spin moment of substitutional metal impurities can be controlled using strain. In particular, we show that although Ni substitutionals are non-magnetic in flat and unstrained graphene, the magnetism of these defects can be activated by applying either uniaxial strain or curvature to the graphene layer. All these results provide key information about formation and control of defect-induced magnetism in graphene and related materials.

Catalytic Water Formation on Platinum: A First-Principles Study

The study of catalytic behavior begins with one seemingly simple process, namely the hydrogenation of O to H2O on platinum. Despite the apparent simplicity its mechanism has been much debated. We have used density functional theory with,gradient corrections to examine microscopic reaction pathways for several elementary steps implicated in this fundamental catalytic process. We find that H2O formation from chemisorbed O and H atoms is a highly activated process. The largest barrier along this route, with a value of similar to1 eV, is the addition of the first H to O to produce OH. Once formed, however, OH groups are easily hydrogenated to H2O with a barrier of similar to0.2 eV. Disproportionation reactions with 1:1 and 2:1 stoichiometries of H2O and O have been examined as alternative routes for OH formation. Both stoichiometries of reaction produce OH groups with barriers that are much lower than that associated with the O + H reaction. H2O, therefore, acts as an autocatalyst in the overall H O formation process. Disproportionation with a 2:1 stoichiometry is thermodynamically and kinetically favored over disproportionation with a l:I stoichiometry. This highlights an additional (promotional) role of the second H2O molecule in this process. In support of our previous suggestion that the key intermediate in the low-temperature H2O formation reaction is a mixed OH and H2O overlayer we find that then is a very large barrier for the dissociation of the second H2O molecule in the 2:1 disproportionation process. We suggest that the proposed intermediate is then hydrogenated to H2O through a very facile proton transfer mechanism.