Electronic-structure-based material descriptors: (in)dependence on self-interaction and Hartree-Fock exchange

Rational design of improved transition metal based materials mostly relies on their electronic structure descriptors, typically estimated by density functional theory and so unduly affected by self-interaction or static correlation errors. Here we show for all 30 transition metals that original or width-corrected d-band centers, and Hilbert transform highest peak descriptors are unaffected by self-interaction, while poor treatment of static correlation by hybrid functionals leads to an unbalanced description. Thus, descriptors have a general validity unbiased by a specific computational method.

A study on adatom transport through (√3 × √3)-R30°-CH3S self-assembled monolayers on Au(111) using first principles calculations

Self-assembled monolayers on Au(111) have outstanding chemical, electrical, and optical properties, and Au adatoms seem to play a key role in these properties. Still, the fundamental understanding of adatom transport inside the self-assembled structure is very thin. In this paper we use first-principles calculations to reveal new details about the migration mechanism of Au adatoms in the presence of a CH3S self-assembled structure on Au(111). We study the inclusion of Au adatoms inside a well-packed (√3 × √3)-R30°-CH3S self-assembled lattice and present atomistic models supporting adatom migration by means of a hopping mechanism between pairs of CH3S species. Our calculations reveal that the transport of Au adatoms is slowed down inside the molecular network where the kinetic barrier for adatom migration is larger than on the clean Au surface. We attribute the hindered mobility of Au adatoms to the fact that adatom transport involves the breaking and making of Au-S bonds. Our results form a basis for further understanding the role played by defect transport in the properties of molecular assemblies.

Long range coupling between defect centres in inorganic nanostructures: valence alternation pairs in nanoscale silica

Valence alternation pair (VAP) states are formed by a closed-shell combination of two space- and charge-separated topological defect centres. These pairs of defects, although historically invoked to explain the electronic properties of bulk inorganic glassy materials (e.g., amorphous silicon dioxide) via the concept of negative-U defects, have more recently been found in a number of theoretical studies of silica surfaces and nanoscale silica clusters. Using density functional theory we systematically probe the structure and internal stability of VAPs in a number of silica nanoclusters with respect to the separation of the two constituent defect centres. We find that VAP states in nanosilica are strongly stabilised by the attractive electrostatic interaction between their separated oppositely charged component defects such that VAPs can persist up to an internal separation of a least 1.5 nanometres. Beyond this distance VAPs become unstable with respect to an open-shell combination of topological defects, virtually indistinguishable from two isolated open-shell defect centres. Finally, we theoretically analyse the possibility of experimental observation of VAP states through their infra-red vibrational spectra.

Mechanisms of Defect Generation and Clustering in CH3S Self-Assembled Monolayers on Au(111)

Periodic density functional calculations probe that step edges play a key role as source of defects during self-assembly. It is shown that the self-assembly process strongly reduces the energy required to strip an atom from the gold surface, locally increasing the concentration of surface defects. The thermodynamic driving force for the atom stripping is considerably more favorable along step-edge lines within the self-assembly than on the higher-coordinated terrace sites. Furthermore, the clustering of surface defects is considered, and we probe that the formation of aggregates of vacancies in the form of vacancy pits significantly stabilizes the self-assembly on the terraces of gold, where the role of the step edges is expected to be less significant. The high stability of pit-like structures arises from a balance between the corrugation and the enhanced bonding of defect-rich substrates. Our results demonstrate the important role that step edges play during assembly and could be very valuable for discovering defect-free assembled structures.

Evidence for the formation of different energetically similar atomic structures in Ag(111)-(square root[7] x square root[7])-R19.1 degrees-CH3S

The atomic structure and thermodynamic stability of Ag(111)(sqrt[7]xsqrt[7])-R19.1 degrees -CH3S has been studied by means of density functional calculations and atomistic first principles thermodynamics. The unreconstructed model and two recently proposed reconstructions have been considered. It is found that, in spite of significant differences in the atomic structure, the different surface models have a very similar surface free energy. It is claimed that the different ordered phases can coexist and that the appearance of one or another depends on the external preparation conditions.

Bulk and surface oxygen vacancy formation and diffusion in single crystals, ultrathin films, and metal grown oxide structures

The neutral oxygen vacancy (OV) energy formation for bulk, subsurface sites at different depths from the surface and various surface sites has been estimated for single crystals, unsupported ultrathin films of MgO, CaO, and BaO, and MgO ultrathin films supported on Ag(001). From the calculated energy barriers for diffusion through the surface and from the surface to the bulk it is found that diffusion is a hindered event, especially for MgO. Nevertheless, diffusion from the terrace to step edges is largely favored while diffusion through terrace sites is less likely and surface to bulk has a very low probability. It is argued that this explains recent scanning tunneling microscopy images for MgO thin films supported on Ag(001) showing OV populating preferentially the step edge sites.

Role of Kinetics in the Selective Surface Oxidations of Transition Metal Carbides

The different oxidation behavior of TiC and VC(100) surfaces by molecular oxygen has been investigated by density functional theory with a slab model. From the thermodynamic stability of the final states that involve dissociated O(2), one cannot well explain the experimental observations. Two different oxidation pathways of TiC and VC(100) surfaces have been explored in this work, and the results indicate that two channels share the same precursor state. However, from the precursor, only the pathway leading to the formation of a C-O bond is energetically feasible for the TiC(100) surface, while on VC(100) the O atoms tend to occupy the metal surface sites due to a smaller energy barrier for this channel. Further band structure calculations reveal that the additional d electron of V atom favors the stability of the molecularly adsorbed species. The oxidation mechanism unveiled from the present calculations clearly evidences that the kinetic effects introduced by one additional d electron of the V atom play a crucial role in explaining the different surface chemistry between TiC and VC (100) surfaces.

A systematic density functional theory study of the electronic structure of bulk and (001) surface of transition-metals carbides

A systematic study of the bulk and surface geometrical and electronic properties of a series of transition-metal carbides (TMC with TM = Ti, V, Zr, Nb, Mo, Hf, Ta, and W) by first-principles methods is presented. It is shown that in these materials the chemical bonding is strongly covalent, the cohesive energies being directly related to the bonding-antibonding gap although the shift of the center of the C(2s) band related peak in the density of states with respect to diamond indicates that some metal to carbon charge transfer does also take place. The (001) face of these metal carbides exhibits a noticeable surface rumpling which grows along the series. It is shown that neglecting surface relaxation results in very large errors on the surface energy and work function. The surface formation induces a significant shift of electronic energy levels with respect to the corresponding values in the bulk. The extent and nature of the shift can be understood from simple bonding-antibonding arguments and is enhanced by the structural rippling of this surface.

The interaction of CO2 with sodium-promoted W(011)

The activation of CO2 by interaction with Na atoms on tungsten was studied in a joint experimental/theoretical effort combining MIES, UPS (HeII) and first principles calculations. Experimentally, both the adsorption of Na on tungsten, followed by CO2 exposure to the Na-modified surface at 80 K, and the adsorption of CO2 on tungsten, followed by Na exposure to the CO2 covered substrate, were studied. Below about 120 K CO2 physisorbs on pure W(011), and the distance between the three main spectral features is as for gas phase CO2 (E(B) = 8.4, 12.1, 14.1 eV). When offered to a Na monolayer (ML) deposited onto W, CO2 is converted into a chemisorbed species. The spectral pattern is different from physisorbed CO2, and the three spectral features are shifted towards lower binding energies (E(B) = 6.3, 10.7, 13.9 eV). The chemisorption continues until all available Na species are converted into Na+ species. Additional CO2 offered to the system becomes physisorbed on top of the chemisorbed species. When a CO2 monolayer, physisorbed on tungsten at 80 K, is exposed to Na, the interaction leads initially to a decrease of the surface work function and to a rigid, global shift of all CO2 induced features towards larger binding energies by about 2 eV. Only beyond a minimum Na coverage of about 0.5 ML, chemisorbed species can be detected. We conclude that, initially, transfer of the Na(3s) electron to the tungsten substrate takes place. Above 0.5 ML Na coverage, back donation of charge to CO2 takes place whereby the physisorbed carbon dioxide species become converted into chemisorbed ones. The experimental results are interpreted with the help of first principle calculations carried out on suitable slab models. The structures and surface binding mode of the chemisorbed CO2 species are described. The calculated density of states for the most stable situations is in qualitative agreement with experimental data.

First principles analysis of the stability and diffusion of oxygen vacancies in metal oxides

Oxygen vacancies in metal oxides are known to determine their chemistry and physics. The properties of neutral oxygen vacancies in metal oxides of increasing complexity (MgO, CaO, alpha-Al2O3, and ZnO) have been studied using density functional theory. Vacancy formation energies, vacancy-vacancy interaction, and the barriers for vacancy migration are determined and rationalized in terms of the ionicity, the Madelung potential, and lattice relaxation. It is found that the Madelung potential controls the oxygen vacancy properties of highly ionic oxides whereas a more complex picture arises for covalent ZnO.

Putting error bars on theAb Initiotheoretical estimates of the magnetic coupling constants: The parent compounds of superconducting cuprates as a case study

The influence of the basis set size and computational method in the calculation of the magnetic coupling constant J is evaluated using a series of cuprate superconductor parent compounds as a case study. The variational DDCI method and an iterative modification, the IDDCI method, are tested, as well as the perturbative CASPT2 method, with two different reference wave functions. Results show that the DDCI magnetic coupling constant is in rather good agreement with the experiment, although it shows a moderate basis set dependency. The IDDCI results are less dependent on the size of the basis set, but slightly overestimate the magnetic coupling constant. CASPT2 results are nearly independent of the chosen basis set. With a minimal active space values are obtained that are about 20% smaller than the DDCI results. The experimental coupling constant can be reproduced when an extended reference wave function is used.

Accurate prediction of large antiferromagnetic interactions in high- T(c) HgBa2Ca(n-1)Cu(n)O(2n+2+delta) ( n = 2,3) superconductor parent compounds

The in-plane nearest-neighbor Heisenberg magnetic coupling constant, J, of La2CuO4, Nd2CuO4, Sr2CuO2Cl2, YBa2Cu3O6, and undoped HgBa(2)Ca(n-1)Cu(n)O(2n+2+delta) ( n = 1,2,3) is calculated from accurate ab initio configuration interaction calculations. For the first four compounds, the theoretical J values are in quantitative agreement with experiment. For the Hg-based compounds the predicted values are -135 meV ( n = 1) and approximately -160 meV ( n = 2,3), the latter being much larger than in previous cases and, for n = 3, increasing with pressure. Nevertheless, the physics governing J in all these layered cuprates appears to be the same. Moreover, calculations suggest a possible relationship between J and T(c).