Towards high-throughput many-body perturbation theory: efficient algorithms and automated workflows

The automation of ab initio simulations is essential in view of performing high-throughput (HT) computational screenings oriented to the discovery of novel materials with desired physical properties. In this work, we propose algorithms and implementations that are relevant to extend this approach beyond density functional theory (DFT), in order to automate many-body perturbation theory (MBPT) calculations. Notably, an algorithm pursuing the goal of an efficient and robust convergence procedure for GW and BSE simulations is provided, together with its implementation in a fully automated framework. This is accompanied by an automatic GW band interpolation scheme based on maximally localized Wannier functions, aiming at a reduction of the computational burden of quasiparticle band structures while preserving high accuracy. The proposed developments are validated on a set of representative semiconductor and metallic systems.

Automated summary

In this work, algorithms and implementations were proposed to automate many-body perturbation theory (MBPT) calculations beyond density functional theory (DFT). An algorithm for efficient and robust convergence of GW and BSE simulations, along with its implementation in a fully automated framework, was provided. Additionally, an automatic GW band interpolation scheme based on maximally localized Wannier functions was introduced to reduce the computational burden of quasiparticle band structures while maintaining high accuracy. These developments were validated on representative semiconductor and metallic systems.