Self-consistent GW method: O(N) algorithm for polarizability and self energy

An efficient implementation of the self-consistent GW method in the FlapwMBPT code https://www.bnl.gov/cmpmsd/flapwmbpt/ is presented. It features the evaluation of polarizability and self-energy which scales linearly with respect to the system size. Altogether the computational time scaling was measured to be between linear and quadratic in the applications to silicon supercells with up to 72 atoms. Application to such materials as paracostibite CoSbS, supercells of La2CuO4 (up to 56 atoms) and SmB6, illustrate the potential of the approach in computational material science. Published by Elsevier B.V.