Accelerating GW-Based Energy Level Alignment Calculations for Molecule-Metal Interfaces Using a Substrate Screening Approach.

The physics of electronic energy level alignment at interfaces formed between molecules and metals can in general be accurately captured by the ab initio GW approach. However, the computational cost of such GW calculations for typical interfaces is significant, given their large system size and chemical complexity. In the past, approximate self-energy corrections, such as those constructed from image-charge models together with gas-phase molecular level corrections, have been used to compute level alignment with good accuracy. However, these approaches often neglect dynamical effects of the polarizability and require the definition of an image plane. In this work, we propose a new approximation to enable more efficient GW-quality calculations of interfaces, where we greatly simplify the calculation of the noninteracting polarizability, a primary bottleneck for large heterogeneous systems. This is achieved by first computing the noninteracting polarizability of each individual component of the interface, e.g., the molecule and the metal, without the use of large supercells, and then using folding and spatial truncation techniques to efficiently combine these quantities. Overall this approach significantly reduces the computational cost for conventional GW calculations of level alignment without sacrificing the accuracy. Moreover, this approach captures both dynamical and nonlocal polarization effects without the need to invoke a classical image-charge expression or to define an image plane. We demonstrate our approach by considering a model system of benzene at relatively low coverage on the aluminum (111) surface. Although developed for such interfaces, the method can be readily extended to other heterogeneous interfaces.