The electronic conductance of a benzene molecule connected to gold electrodes via thiol, thiolate, and amino anchoring groups is calculated using nonequilibrium Green functions in combination with the fully selfconsistent GW approximation. The calculated conductance of benzenedithiol and benzenediamine is five times lower than predicted by standard density functional theory (DFT) in very good agreement with experiments. In contrast, the widely studied benzenedithiolate structure is found to have a significantly higher conductance due to the unsaturated sulfur bonds. These findings suggest that more complex gold/thiolate structures where the thiolate anchors are chemically passivated by Au adatoms are responsible for the measured conductance. Analysis of the energy level alignment obtained with DFT, Hartree-Fock and GW reveals the importance of self-interaction corrections (exchange) on the molecule and dynamical screening at the metal-molecule interface. The main effect of the GW self-energy is to renormalize the level positions, however, its influence on the shape of molecular resonances also affects the conductance. Non-selfconsistent G0W0 calculations, starting from either DFT or Hartree-Fock, yield conductance values within 50% of the selfconsistent GW results.

The electronic conductance of a benzene molecule connected to gold electrodes via thiol, thiolate, or amino anchoring groups is calculated using nonequilibrium Green functions in combination with the fully self-consistent GW approximation for exchange and correlation. The calculated conductance of benzenedithiol and benzenediamine is one-fifth that predicted by standard density functional theory (DFT), in very good agreement with experiments. In contrast, the widely studied benzenedithiolate structure is found to have a significantly higher conductance due to the unsaturated sulfur bonds. These findings suggest that more complex gold-thiolate structures where the thiolate anchors are chemically passivated by Au adatoms are responsible for the measured conductance. Analysis of the energy level alignment obtained with DFT, Hartree-Fock, and GW reveals the importance of self-interaction corrections (exchange) on the molecule and dynamical screening at the metal-molecule interface. The main effect of the GW self-energy is to renormalize the level positions; however, its influence on the shape of molecular resonances also affects the conductance. Non-self-consistent G(0)W(0) calculations, starting from either DFT or Hartree-Fock, yield conductance values within 50% of the self-consistent GW results.

Electronic structure calculations have become an indispensable tool in many areas of materials science and quantum chemistry. Even though the Kohn-Sham formulation of the density-functional theory (DFT) simplifies the many-body problem significantly, one is still confronted with several numerical challenges. In this article we present the projector augmented-wave (PAW) method as implemented in the GPAW program package (https://wiki.fysik.dtu.dk/gpaw) using a uniform real-space grid representation of the electronic wavefunctions. Compared to more traditional plane wave or localized basis set approaches, real-space grids offer several advantages, most notably good computational scalability and systematic convergence properties. However, as a unique feature GPAW also facilitates a localized atomic-orbital basis set in addition to the grid. The efficient atomic basis set is complementary to the more accurate grid, and the possibility to seamlessly switch between the two representations provides great flexibility. While DFT allows one to study ground state properties, time-dependent density-functional theory (TDDFT) provides access to the excited states. We have implemented the two common formulations of TDDFT, namely the linear-response and the time propagation schemes. Electron transport calculations under finite-bias conditions can be performed with GPAW using non-equilibrium Green functions and the localized basis set. In addition to the basic features of the real-space PAW method, we also describe the implementation of selected exchange-correlation functionals, parallelization schemes, ΔSCF-method, x-ray absorption spectra, and maximally localized Wannier orbitals.

We calculate single-particle excitation energies for a series of 33 molecules using fully selfconsistent GW, one-shot G$_0$W$_0$, Hartree-Fock (HF), and hybrid density functional theory (DFT). All calculations are performed within the projector augmented wave (PAW) method using a basis set of Wannier functions augmented by numerical atomic orbitals. The GW self-energy is calculated on the real frequency axis including its full frequency dependence and off-diagonal matrix elements. The mean absolute error of the ionization potential (IP) with respect to experiment is found to be 4.4, 2.6, 0.8, 0.4, and 0.5 eV for DFT-PBE, DFT-PBE0, HF, G$_0$W$_0$[HF], and selfconsistent GW, respectively. This shows that although electronic screening is weak in molecular systems its inclusion at the GW level reduces the error in the IP by up to 50% relative to unscreened HF. In general GW overscreens the HF energies leading to underestimation of the IPs. The best IPs are obtained from one-shot G$_0$W$_0$ calculations based on HF since this reduces the overscreening. Finally, we find that the inclusion of core-valence exchange is important and can affect the excitation energies by as much as 1 eV.

On the basis of first-principles ${G}_{0}{W}_{0}$ calculations we systematically study how the electronic levels of a benzene molecule are renormalized by substrate polarization when physisorbed on different metallic and semiconducting surfaces. The polarization-induced reduction in the energy gap between occupied and unoccupied molecular levels is found to scale with the substrate density of states at the Fermi level (for metals) and substrate band gap (for semiconductors). These conclusions are further supported by self-consistent GW calculations on simple lattice models. By expressing the electron self-energy in terms of the substrate's joint density of states we relate the level shift to the surface electronic structure, thus providing a microscopic explanation of the trends in the GW and ${G}_{0}{W}_{0}$ calculations. While image charge effects are not captured by semilocal and hybrid exchange-correlation functionals, we find that error cancellations lead to remarkably good agreement between the ${G}_{0}{W}_{0}$ and Kohn-Sham energies for the occupied orbitals of the adsorbed molecule.

We present a computational screening study of ternary metal borohydrides for reversible hydrogen storage based on density functional theory. We investigate the stability and decomposition of alloys containing 1 alkali metal atom, Li, Na, or K (M(1)); and 1 alkali, alkaline earth or 3d/4d transition metal atom (M(2)) plus two to five (BH(4))(-) groups, i.e., M(1)M(2)(BH(4))(2-5), using a number of model structures with trigonal, tetrahedral, octahedral, and free coordination of the metal borohydride complexes. Of the over 700 investigated structures, about 20 were predicted to form potentially stable alloys with promising decomposition energies. The M(1)(Al/Mn/Fe)(BH(4))(4), (Li/Na)Zn(BH(4))(3), and (Na/K)(Ni/Co)(BH(4))(3) alloys are found to be the most promising, followed by selected M(1)(Nb/Rh)(BH(4))(4) alloys.