Charge self-consistent many-body corrections using optimized projected localized orbitals

In order for methods combining ab initio density-functional theory and many-body techniques to become routinely used, a flexible, fast, and easy-to-use implementation is crucial. We present an implementation of a general charge self-consistent scheme based on projected localized orbitals in the projector augmented wave framework in the Vienna Ab Initio Simulation Package. We give a detailed description on how the projectors are optimally chosen and how the total energy is calculated. We benchmark our implementation in combination with dynamical mean-field theory: first we study the charge-transfer insulator NiO using a Hartree-Fock approach to solve the many-body Hamiltonian. We address the advantages of the optimized against non-optimized projectors and furthermore find that charge self-consistency decreases the dependence of the spectral function-especially the gap-on the double counting. Second, using continuous-time quantum Monte Carlo we study a monolayer of SrVO₃, where strong orbital polarization occurs due to the reduced dimensionality. Using total-energy calculation for structure determination, we find that electronic correlations have a non-negligible influence on the position of the apical oxygens, and therefore on the thickness of the single SrVO₃ layer.

The cohesive energy of superheavy element copernicium determined from accurate relativistic coupled-cluster theory

The cohesive energy of bulk copernicium is accurately determined using the incremental method within a relativistic coupled-cluster approach. For the lowest energy structure of hexagonal close-packed (hcp) symmetry, we obtain a cohesive energy of -36.3 kJ mol^-1 (inclusion of uncertainties leads to a lower bound of -39.6 kJ mol^-1), in excellent agreement with the experimentally estimated sublimation enthalpy of -38 kJ mol^-1 [R. Eichler et al., Angew. Chem. Int. Ed., 2008, 47, 3262]. At the coupled-cluster singles, doubles and perturbative triples level of theory, we find that the hcp structure is energetically quasi-degenerate with both face-centred and body-centred cubic structures. These results provide a basis for testing various density-functionals, of which the PBEsol functional yields a cohesive energy of -34.1 kJ mol^-1 in good agreement with our coupled-cluster value.

Maximally localized Wannier functions in LaMnO3 within PBE + U, hybrid functionals and partially self-consistent GW: an efficient route to construct ab initio tight-binding parameters for eg perovskites

Using the newly developed VASP2WANNIER90 interface we have constructed maximally localized Wannier functions (MLWFs) for the e(g) states of the prototypical Jahn-Teller magnetic perovskite LaMnO(3) at different levels of approximation for the exchange-correlation kernel. These include conventional density functional theory (DFT) with and without the additional on-site Hubbard U term, hybrid DFT and partially self-consistent GW. By suitably mapping the MLWFs onto an effective e(g) tight-binding (TB) Hamiltonian we have computed a complete set of TB parameters which should serve as guidance for more elaborate treatments of correlation effects in effective Hamiltonian-based approaches. The method-dependent changes of the calculated TB parameters and their interplay with the electron-electron (el-el) interaction term are discussed and interpreted. We discuss two alternative model parameterizations: one in which the effects of the el-el interaction are implicitly incorporated in the otherwise 'noninteracting' TB parameters and a second where we include an explicit mean-field el-el interaction term in the TB Hamiltonian. Both models yield a set of tabulated TB parameters which provide the band dispersion in excellent agreement with the underlying ab initio and MLWF bands.

First-principles study of Cu2ZnSnS4 and the related band offsets for photovoltaic applications

First-principles calculations of the band offsets between Cu(2)ZnSnS(4) (CZTS) and XS (X = Cd, Zn) are performed. While the interface dipole contribution for the band offsets is calculated using the Perdew-Burke-Ernzerhof functional, the Heyd-Scuseria-Ernzerhof hybrid functional is employed to introduce the quasiparticle corrections to the band offsets. The calculated conduction band offset between CZTS and CdS is 0.2 eV, validating CdS for the buffer layer of the CZTS solar cell. The small conduction band offset stems from the band gap narrowing of CdS under the interface strain caused by the lattice misfit with CZTS. A large valence band offset over 0.9 eV between CZTS and ZnS indicates that precipitated ZnS is regarded as an inactive insulator phase in CZTS absorbers.

Cohesive properties and asymptotics of the dispersion interaction in graphite by the random phase approximation

The structural properties of graphite, such as the interlayer equilibrium distance, the elastic constant, and the net layer binding energy, are obtained using the adiabatic-connection fluctuation-dissipation theorem in the random phase approximation. Excellent agreement is found with the available experimental data; however, our computed binding energy of 48 meV per atom is somewhat smaller than the one obtained by quantum Monte Carlo methods. The asymptotic behavior of the interlayer dispersion interaction, previously derived from analytic approximations, is explicitly demonstrated to follow a d-3 behavior at very large distances.

Accurate surface and adsorption energies from many-body perturbation theory

Kohn-Sham density functional theory is the workhorse computational method in materials and surface science. Unfortunately, most semilocal density functionals predict surfaces to be more stable than they are experimentally. Naively, we would expect that consequently adsorption energies on surfaces are too small as well, but the contrary is often found: chemisorption energies are usually overestimated. Modifying the functional improves either the adsorption energy or the surface energy but always worsens the other aspect. This suggests that semilocal density functionals possess a fundamental flaw that is difficult to cure, and alternative methods are urgently needed. Here we show that a computationally fairly efficient many-electron approach, the random phase approximation to the correlation energy, resolves this dilemma and yields at the same time excellent lattice constants, surface energies and adsorption energies for carbon monoxide and benzene on transition-metal surfaces.

Polaronic hole trapping in doped BaBiO3

The present ab initio study shows that in BaBiO3, Bi3+ sites can trap two holes from the valence band to form Bi5+ cations. The trapping is accompanied by large local lattice distortions; therefore the composite particle consisting of the electronic hole and the local lattice phonon field forms a polaron. Our study clearly shows that even sp elements can trap carriers at lattice sites, if local lattice relaxations are sufficiently large to screen the localized hole. The derived model describes all relevant experimental results, and settles the issue of why hole-doped BaBiO3 remains semiconducting upon moderate hole doping.

Step-orientation-dependent oxidation: from 1D to 2D oxides

Using scanning tunneling microscopy and density functional theory, we have studied the initial oxidation of Rh(111) surfaces with two types of straight steps, having {100} and {111} microfacets. The one-dimensional (1D) oxide initially formed at the steps acts as a barrier impeding formation of the 2D oxide on the (111) terrace behind it. We demonstrate that the details of the structure of the 1D oxide govern the rate of 2D oxidation and discuss implications for oxidation of nanoparticles.

Accurate band gaps and dielectric properties from one-electron theories (abstract only)

For semiconductor modeling, a major shortcoming of density functional theory is that the predicted band gaps are usually significantly too small. It is generally argued that this shortcoming is related to the fact that density functional theory is a ground state theory, and as a result, one is not allowed to associate the one-electron energies with the energies of quasi-particles. Although this fundamental objection is certainly correct, the modeling of the positioning of donor and acceptor levels in semiconductors faces serious limitations with present density functionals. Several solutions to this problem have been suggested. A particular attractive and fairly simple one is the inclusion of a small fraction of the non-local exchange in the Hamiltonian (hybrid functionals). This approach leads to sensible band gaps for most semiconductors, but fails for ionic solids. A more reliable approach is via many-electron Green's function techniques, which have made tremendous advances in recent years. Here GW calculations in various flavors are presented for small gap and large gap systems, comprising typical semiconductors (Si, SiC, GaAs, GaN, ZnO, ZnS, CdS and AlP), small gap semiconductors (PbS, PbSe, PbTe), insulators (C, BN, MgO, LiF) and noble gas solids (Ar, Ne). The general finding is that single-shot G(0)W(0) calculations based on wavefunctions obtained from conventional density functional theory yield too small band gaps, whereas G(0)W(0) calculations following hybrid functional calculations tend to overestimate the band gaps by roughly the same amount. This is at first sight astonishing, since the hybrid functionals yield very good band gaps themselves. The contradiction is resolved by showing that the inclusion of the attractive electron-hole interactions (excitonic effects) is required to obtain good static and dynamic dielectric functions using hybrid functionals. The corrections are usually incorporated in GW calculations using 'vertex corrections', and in fact inclusion of these vertex corrections rectifies the predicted band gaps. Finally, in order to remove the dependence on the initial wavefunctions, self-consistent GW calculations are presented, again including an approximate treatment of vertex corrections. The results are in excellent agreement with experiment, with a few per cent deviation for all materials considered. We conclude that predictive band gap engineering is now possible with the theoretical description approaching experimental accuracy.

Hybrid functionals applied to extended systems

We present an overview of the description of structural, thermochemical, and electronic properties of extended systems using several well known hybrid Hartree-Fock/density-functional-theory functionals (PBE0, HSE03, and B3LYP). In addition we address a few aspects of the evaluation of the Hartree-Fock exchange interactions in reciprocal space, relevant to all methods that employ a plane wave basis set and periodic boundary conditions.

Accurate quasiparticle spectra from self-consistent GW calculations with vertex corrections

Self-consistent GW calculations, maintaining only the quasiparticle part of the Green's function G, are reported for a wide class of materials, including small gap semiconductors and large gap insulators. We show that the inclusion of the attractive electron-hole interaction via an effective nonlocal exchange correlation kernel is required to obtain accurate band gaps in the framework of self-consistent GW calculations. If these are accounted for via vertex corrections in W, the band gaps are found to be within a few percent of the experimental values.

Nanotemplate with holes: ultrathin alumina on Ni3Al(111)

We have determined the structure of the ultrathin (sqrt[67] x sqrt[67])R12.2 degrees aluminum oxide on Ni3Al(111) by a combination of scanning tunneling microscopy and density functional theory. In addition to other local defects, the main structural feature of the unit cell is a 0.4-nm-diameter hole reaching down to the metal substrate. Understanding the structure and metal growth on this oxide allows us to use it as a template for growing highly regular arrays of nanoparticles.

Relaxed core projector-augmented-wave method

We extend the full-potential projector-augmented-wave method beyond the frozen core approximation, i.e., include the self-consistent optimization of the core charge density, in such a manner that the valence wave functions remain orthogonal to the core. The method consists of an on-the-fly repseudization of the all-electron problem, solving for the self-consistent core charge density within a spherical approximation. The key ideas in our procedure are to keep the projector functions fixed throughout the electronic minimization and to derive the new pseudopartial waves from these original projector functions, at each step of the electronic minimization procedure. Results of relaxed core calculations for atomic interconfigurational energies, structural energy differences between bulk phases of Fe, atomization energies of a subset of Pople's G2-1 set, and the Rh 3d surface core level shifts for the (log3 x log3)-Rh(111) surface at 1/3 CO coverage are presented.

Initial oxidation of the Rh(110) surface: Ordered adsorption and surface oxide structures

The initial oxidation of the Rh(110) surface was studied by scanning tunneling microscopy, core level spectroscopy, and density functional theory. The experiments were carried out exposing the Rh(110) surface to molecular or atomic oxygen at temperatures in the 500-700 K range. In molecular oxygen ambient, the oxidation terminates at oxygen coverage close to a monolayer with the formation of alternating islands of the (10x2) one-dimensional surface oxide and (2x1)p2mg adsorption phases. The use of atomic oxygen facilitates further oxidation until a structure with a c(2x4) periodicity develops. The experimental and theoretical results reveal that the c(2x4) structure is a "surface oxide" very similar to the hexagonal O-Rh-O trilayer structures formed on the Rh(111) and Rh(100) substrates. Some of the experimentally found adsorption phases appear unstable in the phase diagram predicted by thermodynamics, which might reflect kinetic hindrance. The structural details, core level spectra, and stability of the surface oxides formed on the three basal planes are compared with those of the bulk RhO2 and Rh2O3.

Oxygen-deficient line defects in an ultrathin aluminum oxide film

A model for the straight antiphase domain boundary of the ultrathin aluminum oxide film on the NiAl(110) substrate is derived from scanning tunneling microscopy measurements and density-functional theory calculations. Although the local bonding environment of the perfect film is maintained, the structure is oxygen deficient and possesses a favorable adsorption site. The domain boundary exhibits a downwards band bending and three characteristic unoccupied electronic states, in excellent agreement with scanning tunneling spectroscopy measurements.

Screened hybrid density functionals applied to solids

Hybrid Fock exchange/density functional theory functionals have shown to be very successful in describing a wide range of molecular properties. For periodic systems, however, the long-range nature of the Fock exchange interaction and the resultant large computational requirements present a major drawback. This is especially true for metallic systems, which require a dense Brillouin zone sampling. Recently, a new hybrid functional [HSE03, J. Heyd, G. E. Scuseria, and M. Ernzerhof, J. Chem. Phys. 118, 8207 (2003)] that addresses this problem within the context of methods that evaluate the Fock exchange in real space was introduced. We discuss the advantages the HSE03 functional brings to methods that rely on a reciprocal space description of the Fock exchange interaction, e.g., all methods that use plane wave basis sets. Furthermore, we present a detailed comparison of the performance of the HSE03 and PBE0 functionals for a set of archetypical solid state systems by calculating lattice parameters, bulk moduli, heats of formation, and band gaps. The results indicate that the hybrid functionals indeed often improve the description of these properties, but in several cases the results are not yet on par with standard gradient corrected functionals. This concerns in particular metallic systems for which the bandwidth and exchange splitting are seriously overestimated.

Structure of Ag(111)-p(4 x 4)-O: no silver oxide

The structure of the oxygen-induced p(4 x 4) reconstruction of Ag(111) is determined by a combination of scanning tunneling microscopy, surface x-ray diffraction, core level spectroscopy, and density functional theory. We demonstrate that all previous models of this surface structure are incorrect and propose a new model which is able to explain all our experimental findings but has no resemblance to bulk silver oxide. We also shed some light on the limitations of current density functional theories and the potential role of van der Waals interactions in the stabilization of oxygen-induced surface reconstructions of noble metals.

Ab initio density-functional study of the bridging addition of acrylonitrile on the Si(100) surface

Acrylonitrile may react with the Si(100) surface in a number of ways. Fifteen different configurations have been identified. This study which complements an earlier study devoted to cycloadditions and bonds involving a single atom deals essentially with configurations involving cumulative double bonds. Factors affecting the relative stability of various mesomeric forms are shown to be the proximity to tetrahedral geometry and the compatibility of the dipole form with buckled structures. Tripods are shown to be a way to accommodate a priori unfavorable structures.

High-coverage oxygen structures on Rh111: adsorbate repulsion and site preference is not enough

A new O induced structure on Rh(111) displaying a (2 sqaureroot[3] x 2sqaureroot[3])R30 degrees periodicity with an oxygen coverage of 2/3 has been studied by high resolution core level spectroscopy, scanning tunneling microscopy, and density functional theory. Although O favors fcc hollow sites in all other known phases, it occupies both fcc and hcp sites in this structure, which cannot be explained by pairwise adsorbate repulsion only. Both the (2sqaureroot[3] x 2sqaureroot[3])R30 degrees and (2 x 2)-3O structures also exemplify that density-of-states contrast can lead to oxygen adatoms appearing as protrusions in scanning tunneling microscopy images.

Fullerene quantum gyroscope

We report the observation of quantized rotational states of a diatomic C2 unit in solid endohedral fullerene C(2)Sc(2)@C(84). The rotational transitions induce a periodic line pattern in the low energy Raman spectrum. The rotational constant B and the C-C distance were found to be 1.73 cm(-1) and 0.127 nm, respectively. Density functional calculations revealed an intrinsic rotational barrier of the order of only a few meV for the C2 unit. The Schrödinger equation involving the potential barrier was solved and the Raman tensor matrix elements were calculated, yielding good quantitative agreement with the experiment. To our best knowledge this is the first intrinsic rotational spectrum of a diatomic plane molecular rotor.

Planar vanadium oxide clusters: two-dimensional evaporation and diffusion on Rh(111)

The formation of novel vanadium oxide cluster molecules by oxidative two-dimensional evaporation from vanadium oxide nanostructures is reported on a Rh(111) metal surface. The structure and stability of the planar V6O12 clusters and the physical origin of their 2D evaporation process have been elucidated by high-resolution scanning tunneling microscopy (STM) and ab initio density functional theory calculations. The surface diffusion of the clusters has been followed in elevated-temperature STM experiments, and the diffusion parameters have been extracted, indicating diffusion by hopping of the entire surface stabilized cluster units.

Self-limited growth of a thin oxide layer on Rh(111)

The oxidation of the Rh(111) surface at oxygen pressures from 10(-10) mbar to 0.5 bar and temperatures between 300 and 900 K has been studied on the atomic scale using a multimethod approach of experimental and theoretical techniques. Oxidation starts at the steps, resulting in a trilayer O-Rh-O surface oxide which, although not thermodynamically stable, prevents further oxidation at intermediate pressures. A thick corundum like Rh2O3 bulk oxide is formed only at significantly higher pressures and temperatures.

Novel stabilization mechanism on polar surfaces: ZnO(0001)-Zn

The (1x1) terminated (0001)-Zn surface of wurtzite ZnO was investigated with scanning tunneling microscopy. The surface is characterized by the presence of nanosized islands with a size-dependent shape and triangular holes with single-height, O-terminated step edges. It is proposed that the resulting overall decrease of the surface Zn concentration stabilizes this polar surface. Ab initio calculations of test geometries predict triangularly shaped reconstructions over a wide range of oxygen and hydrogen chemical potentials. The formation of these reconstructions appears to be electrostatically driven.

Reversible dynamic behavior in catalyst systems: oscillations of structure and morphology

Exposing vanadium oxide nanoparticles on a Pd(111) surface to reducing conditions is shown to cause a spreading of the oxide over the metal until a reduced oxide phase covers the entire surface. Reoxidation reverses this process and oxide island structures and bare metal patches are reestablished. The physical origin of this wetting-dewetting process is revealed at the atomic level by in situ variable-temperature scanning tunneling microscopy and in terms of a surface oxide phase stability diagram, as calculated by density functional theory as a function of the chemical potential of oxygen and the vanadium concentration.

Two-dimensional oxide on Pd(111)

The oxidation of Pd(111) leads to an incommensurate surface oxide, which was studied by the use of scanning tunneling microscopy, surface x-ray diffraction, high resolution core level spectroscopy, and density functional calculations. A combination of these methods reveals a two-dimensional structure having no resemblance to bulk oxides of Pd. Our study also demonstrates how the atomic arrangement of a nontrivial incommensurate surface can be solved by molecular dynamics in a case where experimental techniques alone give no solution.

Phonon softening in metallic nanotubes by a Peierls-like mechanism

The radial dependency of the vibrational frequencies of single-wall carbon nanotubes in the G band (1500-1600 cm(-1)) is studied by density functional theory. In metallic nanotubes, a mode with A1 symmetry is found to be significantly softer than the corresponding mode in insulating tubes or graphite. The mechanism that leads to the mode softening is explored. It is reminiscent of the driving force inducing Peierls distortions. At ambient temperature, the energy gained by opening the gap is, however, not sufficient for a static lattice distortion. Instead the corresponding vibrational frequency is lowered.

Coexistence of atomic and molecular chemisorption states: H(2)/Pd(210)

A novel H(2) molecular adsorption state on metal surfaces has been detected by temperature-programmed desorption and electron energy loss spectroscopy experiments of the H(2)/Pd(210) system. The molecular nature of this state has been verified by isotope exchange experiments. This molecular state leads to a decrease of the surface work function while atomic hydrogen on Pd(210) causes an increase. Ab initio total-energy calculations have confirmed all experimental findings. Through these calculations the microscopic nature of this novel molecular adsorption state could be identified; it turns out that this state is stabilized by the presence of atomic hydrogen on the Pd(210) surface.

Novel interface-mediated metastable oxide phases: vanadium oxides on Pd(111)

In the growth process of ultrathin films of vanadium oxides on Pd(111), a sequence of novel oxide phases with layer-dependent structures and oscillating oxidation states has been detected experimentally and understood theoretically. These phases are interface mediated and metastable with respect to further oxide growth. Transformation into the stable oxide configuration occurs beyond a critical thickness, where energetics combined with kinetic limitations determine the oxide multilayer structure.

Adsorption sites and ligand effect for CO on an alloy surface: a direct view

CO adsorption on a PtCo(111) surface was studied by scanning tunneling microscopy. Comparison of images with chemical contrast of Pt and Co and images showing the CO molecules indicates that CO resides exclusively on top of Pt sites and never on Co. CO bonding is highly sensitive to the chemical environment. The probability to find CO on a Pt atom increases drastically with the number of its Co nearest neighbors. Ab initio calculations show that this ligand effect is due to different positions of the center of the Pt d band.