Large-scale GW software development

Electronic excitations are important in understanding and designing many functional materials. In terms of {\it ab initio} methods, the GW and Bethe-Saltpeter Equation (GW-BSE) beyond DFT methods have proved successful in describing excited states in many materials. However, the heavy computational loads and large memory requirements have hindered their routine applicability by the materials physics community. We summarize some of our collaborative efforts to develop a new software framework designed for GW calculations on massively parallel supercomputers. Our GW code is interfaced with the plane-wave pseudopotential {\it ab initio} molecular dynamics software ``OpenAtom'' which is based on the Charm++ parallel library [1]. The computation of the electronic polarizability is one of the most expensive parts of any GW calculation. We describe our strategy that uses a real-space representation to avoid the large number of fast Fourier transforms (FFTs) common to most GW methods. We also describe an eigendecomposition of the plasmon modes from the resulting dielectric matrix that enhances efficiency. [1] Bohm {\it et al.}, IBM J. RES. \& DEV. vol. 52 no. 1/2, 2008