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

Miki Bonacci,Junfeng Qiao,Nicola Spallanzani,Antimo Marrazzo,Giovanni Pizzi,Elisa Molinari,Daniele Varsano,Andrea Ferretti,Deborah Prezzi
DOI: https://doi.org/10.1038/s41524-023-01027-2
2023-01-18
Abstract: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, a novel 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.
Materials Science
What problem does this paper attempt to address?
### What problems does this paper attempt to solve? This paper aims to solve the challenges in the automation of many - body perturbation theory (MBPT) calculations and high - throughput screening, especially in GW approximation and Bethe - Salpeter equation (BSE) calculations. Specifically, the authors raise the following problems and attempt to solve them: 1. **Automation and high - throughput calculations**: - **Problem**: Current material discovery depends on the combination of experimental and computational methods. Among computational methods, high - throughput (HT) computational screening is an important tool for accelerating the discovery of new materials. However, most existing automated workflows are mainly based on density functional theory (DFT). For many - body perturbation theories (such as GW - BSE) that require a more accurate description of excited - state properties, automation and high - throughput calculations still face huge challenges. - **Solution**: The authors propose a new set of algorithms and implementation schemes to extend the automated method to many - body perturbation theory calculations, especially GW and BSE simulations. 2. **Efficient convergence algorithms**: - **Problem**: In GW and BSE calculations, the selection of parameters and convergence is a complex and time - consuming process. Due to the interdependence between these parameters, the traditional method of converging parameters one by one is not only cumbersome but also computationally expensive. - **Solution**: The authors develop a new efficient convergence algorithm that can reduce the consumption of computational resources while ensuring accuracy. This algorithm achieves rapid convergence by automatically selecting appropriate parameter combinations and has been integrated into the AiiDA framework. 3. **Automated band - structure interpolation**: - **Problem**: The calculation of quasiparticle band - structure usually requires a large amount of computational resources, especially when dealing with complex materials. To reduce the computational burden while maintaining high accuracy, an effective band - structure interpolation method is required. - **Solution**: The authors introduce an automatic GW band - structure interpolation scheme based on maximally - localized Wannier functions. This method can maintain high accuracy while significantly reducing the computational cost. 4. **Cross - software interoperability**: - **Problem**: Many - body perturbation theory calculations usually involve multiple steps, and each step may use different software tools, which increases the difficulty of automation and interoperability. - **Solution**: The authors achieve seamless interoperability between different software codes through the AiiDA platform, including from the initial DFT calculation to the final GW - BSE calculation and its post - processing. ### Summary The main goal of this paper is to develop efficient automated workflows and algorithms to achieve high - throughput screening in many - body perturbation theory (especially GW and BSE) calculations. Through these improvements, researchers can explore the electronic structures and optical properties of new materials more quickly and accurately, thus accelerating innovation and development in the field of materials science.