The electronic properties and the function of hybrid inorganic-organic systems (HIOS) are intimately linked to their interface geometry. Here we show that the inclusion of the many-body collective response of the substrate electrons inside the inorganic bulk enables us to reliably predict the HIOS geometries and energies. This is achieved by the combination of dispersion-corrected density-functional theory (the DFT+ van der Waals approach) [Phys. Rev. Lett. 102, 073005 (2009)], with the Lifshitz-Zaremba-Kohn theory for the nonlocal Coulomb screening within the bulk. Our method yields geometries in remarkable agreement (≈0.1 Å) with normal incidence x-ray standing wave measurements for the 3, 4, 9, 10-perylene-tetracarboxylic acid dianhydride (C(24)O(6)H(8), PTCDA) molecule on Cu(111), Ag(111), and Au(111) surfaces. Similarly accurate results are obtained for xenon and benzene adsorbed on metal surfaces.

The adsorption of benzene on metal surfaces is an important benchmark system for hybrid inorganic/organic interfaces. The reliable determination of the interface geometry and binding energy presents a significant challenge for both theory and experiment. Using the Perdew-Burke-Ernzerhof (PBE), PBE+vdW (van der Waals) and the recently developed PBE+vdW surf (densityfunctional theory with vdW interactions that include the collective electronic response of the substrate) methods, we calculated the structures and energetics for benzene on transition-metal surfaces: Cu, Ag, Au, Pd, Pt, Rh and Ir. Our calculations demonstrate that vdW interactions increase the binding energy by more than 0.70 eV for physisorbed systems (Cu, Ag and Au) and by an even larger amount for strongly bound systems (Pd, Pt, Rh and Ir). The collective response of the substrate electrons captured via the vdW surf method plays a significant role for most substrates, shortening the equilibrium distance by 0.25 Å for Cu and decreasing the binding energy by 0.27 eV for Rh. The reliability of our results is assessed by comparison with calculations using the random-phase approximation including renormalized single excitations, often hard to interpret, or even lacking. For example, due to the relative difficulty of controlling and measuring weakly bound systems, no experimental adsorption height has been reported so far for Bz physisorbed on noble metals. The only experimentally deduced adsorption height, to the best of our knowledge, was determined for the disordered Bz chemisorbed on the Pt(111) surface at a coverage close to or less than one [36].The binding energy, which reflects the strength of the interaction between an adsorbate and the substrate, is another key parameter for the description of HIOS. Experimental binding energies are mainly obtained by temperature-programmed desorption (TPD) [22][23][24][25]27] and microcalorimetry measurements [32,[37][38][39]. TPD is the most extensively used method for determining the kinetic and thermodynamic parameters of desorption processes and decomposition reactions. The desorbing molecular species are selected by their mass, while the amount of adsorbate is determined by integrating the peaks of the desorption spectrum. The Redhead formula is typically used to calculate the adsorption energy based on three parameters: the desorption temperature, the heating rate and a pre-exponential factor [40]. The wide range of empirical pre-exponential factors (10 13 -10 19 s −1 ) that are typically used for molecular desorption may cause a notable uncertainty in the determined binding energy [41][42][43]. TPD experiments have been carried out to study the interaction of Bz with the Cu [13], Ag [24] and Au surfaces [29]. However, special attention must be paid to the interpretation of TPD spectra for the Pd, Pt and Rh surfaces, because the adsorbed Bz molecules may decompose during heating, in particular at low coverage [44]. Here, we revisit the adsorption energies from the measured TPD spectra for Bz on Cu(111), A...

Adsorption geometry and stability of organic molecules on surfaces are key parameters that determine the observable properties and functions of hybrid inorganic/organic systems (HIOSs). Despite many recent advances in precise experimental characterization and improvements in rst-principles electronic structure methods, reliable databases of structures and energetics for large adsorbed molecules are largely amiss. In this review, we present such a database for a range of molecules adsorbed on metal single-crystal surfaces. The systems we analyze include noble-gas atoms, conjugated aromatic molecules, carbon nanostructures, and heteroaromatic compounds adsorbed on ve dierent metal surfaces. The overall objective is to establish a diverse benchmark dataset that enables an assessment of current and future electronic structure methods, and motivates further experimental studies that provide ever more reliable data. Specically, the benchmark structures and energetics from experiment are here compared with the recently developed van der Waals (vdW) inclusive density-functional theory (DFT) method, DFT+vdW^surf. In comparison to 23 adsorption heights and 17 adsorption energies from experiment we nd a mean average deviation of 0.06 Å and 0.16 eV, respectively. This conrms the DFT+vdW^surf method as an accurate and ecient approach to treat HIOSs. A detailed discussion identies remaining challenges to be addressed in future development of electronic structure methods, for which the here presented benchmark database may serve as an important reference

Modeling the adsorption of atoms and molecules on surfaces requires ecient electronic structure methods that are able to capture both covalent and non-covalent interactions in a reliable manner. In order to tackle this problem, we have developed a method within density-functional theory (DFT) to model screened van der Waals interactions (vdW) for atoms and molecules on surfaces (the so-called DFT+vdWsurf method). The relatively high accuracy of the DFT+vdWsurf method in the calculation of both adsorption distances and energies, as well as the high degree of its reliability across wide range of adsorbates, indicates the importance of the collective electronic ects within the extended substrate for the calculation of the vdW energy tail. We examine in detail the theoretical background of the method and assess its performance for adsorption phenomena including the physisorption of Xe on selected close-packed transition metal surfaces and 3,4,9,10-perylene-tetracarboxylic acid dianhydride (PTCDA) on Au(111). We also address the performance of DFT+vdWsurf in the case of non-close-packed surfaces by studying the adsorption of Xe on Cu(110) and the interfaces formed by the adsorption of a PTCDA monolayer on the Ag(111), Ag(100), and Ag(110) surfaces. We conclude by discussing outstanding challenges in the modeling of vdW interactions for studying atomic and molecular adsorbates on inorganic substrates

A correct description of electronic exchange and correlation effects for molecules in contact with extended (metal) surfaces is a challenging task for first-principles modeling. In this work we demonstrate the importance of collective van der Waals dispersion effects beyond the pairwise approximation for organic-inorganic systems on the example of atoms, molecules, and nanostructures adsorbed on metals. We use the recently developed many-body dispersion (MBD) approach in the context of density-functional theory [Phys. Rev. Lett. 108, 236402 (2012); J. Chem. Phys. 140, 18A508 (2014)] and assess its ability to correctly describe the binding of adsorbates on metal surfaces. We briefly review the MBD method and highlight its similarities to quantum-chemical approaches to electron correlation in a quasiparticle picture. In particular, we study the binding properties of xenon, 3,4,9,10-perylene-tetracarboxylic acid (PTCDA), and a graphene sheet adsorbed on the Ag(111) surface. Accounting for MBD effects we are able to describe changes in the anisotropic polarizability tensor, improve the description of adsorbate vibrations, and correctly capture the adsorbate-surface interaction screening. Comparison to other methods and experiment reveals that inclusion of MBD effects improves adsorption energies and geometries, by reducing the overbinding typically found in pairwise additive dispersion-correction approaches.