Lattice relaxations at substitutional impurities in semiconductors

First-principles calculations of hybrid inorganic-organic interfaces: from state-of-the-art to best practice

The computational characterization of inorganic-organic hybrid interfaces is arguably one of the technically most challenging applications of density functional theory. Due to the fundamentally different electronic properties of the inorganic and the organic components of a hybrid interface, the proper choice of the electronic structure method, of the algorithms to solve these methods, and of the parameters that enter these algorithms is highly non-trivial. In fact, computational choices that work well for one of the components often perform poorly for the other. As a consequence, default settings for one materials class are typically inadequate for the hybrid system, which makes calculations employing such settings inefficient and sometimes even prone to erroneous results. To address this issue, we discuss how to choose appropriate atomistic representations for the system under investigation, we highlight the role of the exchange-correlation functional and the van der Waals correction employed in the calculation and we provide tips and tricks how to efficiently converge the self-consistent field cycle and to obtain accurate geometries. We particularly focus on potentially unexpected pitfalls and the errors they incur. As a summary, we provide a list of best practice rules for interface simulations that should especially serve as a useful starting point for less experienced users and newcomers to the field.

Final-State Simulations of Core-Level Binding Energies at Metal-Organic Hybrid Interfaces: Artifacts Caused by Spurious Collective Electrostatic Effects

Core-level energies are frequently calculated to explain the X-ray photoelectron spectra of metal-organic hybrid interfaces. The current paper describes how such simulations can be flawed when modeling interfaces between physisorbed organic molecules and metals. The problem occurs when applying periodic boundary conditions to correctly describe extended interfaces and simultaneously considering core hole excitations in the framework of a final-state approach to account for screening effects. Since the core hole is generated in every unit cell, an artificial dipole layer is formed. In this work, we study methane on an Al(100) surface as a deliberately chosen model system for hybrid interfaces to evaluate the impact of this computational artifact. We show that changing the supercell size leads to artificial shifts in the calculated core-level energies that can be well beyond 1 eV for small cells. The same applies to atoms at comparably large distances from the substrate, encountered, for example, in extended, upright-standing adsorbate molecules. We also argue that the calculated work function change due to a core-level excitation can serve as an indication for the occurrence of such an artifact and discuss possible remedies for the problem.

Space-charge transfer in hybrid inorganic-organic systems

We discuss density functional theory calculations of hybrid inorganic-organic systems that explicitly include the global effects of doping (i.e., position of the Fermi level) and the formation of a space-charge layer. For the example of tetrafluoro-tetracyanoquinodimethane on the ZnO(0001[over ¯]) surface we show that the adsorption energy and electron transfer depend strongly on the ZnO doping. The associated work function changes are large, for which the formation of space-charge layers is the main driving force. The prominent doping effects are expected to be quite general for charge-transfer interfaces in hybrid inorganic-organic systems and important for device design.

Concentration of vacancies at metal-oxide surfaces: case study of MgO(100)

We investigate the effects of doping on the formation energy and concentration of oxygen vacancies at a metal-oxide surface, using MgO(100) as an example. Our approach employs density-functional theory, where the performance of the exchange-correlation functional is carefully analyzed, and the functional is chosen according to a condition on density-functional theory ionization energies. The approach is further validated by coupled-cluster calculations, including single, double, and perturbative triple substitutions, for embedded clusters. We demonstrate that the concentration of oxygen vacancies at a doped oxide surface is largely determined by the formation of a macroscopically extended space-charge region.