Simple, yet predictive bonding models are essential achievements of chemistry. In the solid state, in particular, they often appear in the form of visual bonding indicators. Because the latter require the crystal orbitals to be constructed from local basis sets, the application of the most popular density-functional theory codes (namely, those based on plane waves and pseudopotentials) appears as being ill-fitted to retrieve the chemical bonding information. In this paper, we describe a way to re-extract Hamilton-weighted populations from plane-wave electronic-structure calculations to develop a tool analogous to the familiar crystal orbital Hamilton population (COHP) method. We derive the new technique, dubbed "projected COHP" (pCOHP), and demonstrate its viability using examples of covalent, ionic, and metallic crystals (diamond, GaAs, CsCl, and Na). For the first time, this chemical bonding information is directly extracted from the results of plane-wave calculations.
Quantum-chemical computations of solids benefit enormously from numerically efficient plane-wave (PW) basis sets, and together with the projector augmented-wave (PAW) method, the latter have risen to one of the predominant standards in computational solid-state sciences. Despite their advantages, plane waves lack local information, which makes the interpretation of local densities-of-states (DOS) difficult and precludes the direct use of atom-resolved chemical bonding indicators such as the crystal orbital overlap population (COOP) and the crystal orbital Hamilton population (COHP) techniques. Recently, a number of methods have been proposed to overcome this fundamental issue, built around the concept of basis-set projection onto a local auxiliary basis. In this work, we propose a novel computational technique toward this goal by transferring the PW/PAW wavefunctions to a properly chosen local basis using analytically derived expressions. In particular, we describe a general approach to project both PW and PAW eigenstates onto given custom orbitals, which we then exemplify at the hand of contracted multiple-ζ Slater-type orbitals. The validity of the method presented here is illustrated by applications to chemical textbook examples-diamond, gallium arsenide, the transition-metal titanium-as well as nanoscale allotropes of carbon: a nanotube and the C60 fullerene. Remarkably, the analytical approach not only recovers the total and projected electronic DOS with a high degree of confidence, but it also yields a realistic chemical-bonding picture in the framework of the projected COHP method.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.