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.

Free energy and molecular dynamics calculations for the cubic-tetragonal phase transition in zirconia

The high-temperature cubic-tetragonal phase transition of pure stoichiometric zirconia is studied by molecular dynamics (MD) simulations and within the framework of the Landau theory of phase transformations. The interatomic forces are calculated using an empirical, self-consistent, orthogonal tight-binding model, which includes atomic polarizabilities up to the quadrupolar level. A first set of standard MD calculations shows that, on increasing temperature, one particular vibrational frequency softens. The temperature evolution of the free-energy surfaces around the phase transition is then studied with a second set of calculations. These combine the thermodynamic integration technique with constrained MD simulations. The results seem to support the thesis of a second-order phase transition but with unusual, very anharmonic behavior above the transition temperature.

Density functional calculations for polymers and clusters - progress and limitations

The applicability of density functional (DF) methods has progressed greatly since the first workshop of this series ten years ago. Applications that show both the successes and the limitations can be found in the fields of: (a) the structures of the isomers of atomic clusters. and (b) the structure of organic molecules and polymers, and their reactions with additional molecules. We shall review some of the results and the lessons to be learned from them. (C) 2001 Elsevier Science B.V. All rights reserved.

Optimal basis set for electronic structure calculations in periodic systems

An efficient method for calculating the electronic structure of systems that need a very fine sampling of the Brillouin zone is presented. The method is based on the variational optimization of a single (i.e., common to all points in the Brillouin zone) basis set for the expansion of the electronic orbitals. Considerations from k.p-approximation theory help to understand the efficiency of the method. The accuracy and the convergence properties of the method as a function of the optimal basis set size are analyzed for a test calculation on a 16-atom Na supercell.

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.

Interactions of 1-Methylimidazole with UO2(CH3CO2)2 and UO2(NO3)2: Structural, Spectroscopic and Theoretical Evidence for Imidazole Binding to the Uranyl Ion

The first definitive high-resolution single-crystal X-ray structure for the coordination of the 1-methylimidazole (Meimid) ligand to UO2(Ac)2 (Ac = CH3CO2) is reported. The crystal structure evidence is confirmed by IR, Raman, and UV-vis spectroscopic data. Direct participation of the nitrogen atom of the Meimid ligand in binding to the uranium center is confirmed. Structural analysis at the DFT (B3LYP) level of theory showed a conformational difference of the Meimid ligand in the free gas-phase complex versus the solid state due to small energetic differences and crystal packing effects. Energetic analysis at the MP2 level in the gas phase supported stronger Meimid binding over H2O binding to both UO2(Ac)2 and UO2(NO3)2. In addition, self-consistent reaction field COSMO calculations were used to assess the aqueous phase energetics of combination and displacement reactions involving H2O and Meimid ligands to UO2R2 (R = Ac, NO3). For both UO2(NO3)2 and UO2(Ac)2, the displacement of H2O by Meimid was predicted to be energetically favorable, consistent with experimental results that suggest Meimid may bind uranyl at physiological pH. Also, log(Knitrate/KAc) calculations supported experimental evidence that the binding stoichiometry of the Meimid ligand is dependent upon the nature of the reactant uranyl complex. These results clearly demonstrate that imidazole binds to uranyl and suggest that binding of histidine residues to uranyl could occur under normal biological conditions.

A quantitative determination of reaction mechanisms from density functional theory calculations: Fischer-Tropsch synthesis on flat and stepped cobalt surfaces

We systematically investigated the mechanism of the C-1 + C-1 coupling reactions using density functional theory. The activation energies of C-1 + C-1 coupling and carbon hydrogenation reactions on both flat and stepped surfaces were calculated and analyzed. Moreover, the coverages of adsorbed C-1 species were estimated, and the reaction rates of all possible C-1 + C-1 coupling pathways were quantitatively evaluated. The results suggest that the reactions of CH2 + CH2 and CH3 + C at steps are most likely to be the key C-1 + C-1 coupling steps in FT synthesis on Co catalysts. The reactions of C-2 + C-1 and C-3 + C-1 coupling also were studied; the results demonstrate that in addition to the pathways of RCH + CH2 and RCH2 + C, the coupling of RC + C and RC + CH also may contribute to the chain growth after C-1. (C) 2008 Elsevier Inc. All rights reserved.