| Popis: |
We present a massively parallel, GPU-accelerated implementation of the Bethe-Salpeter equation (BSE) for the calculation of the vertical excitation energies (VEEs) and optical absorption spectra of condensed and molecular systems, starting from single-particle eigenvalues and eigenvectors obtained with density functional theory. The algorithms adopted here enable one to circumvent the slowly converging sums over empty and occupied states and the inversion of large dielectric matrices, through a density matrix perturbation theory approach and a low-rank decomposition of the screened Coulomb interaction, respectively. In addition, we achieve computational savings by exploiting the nearsightedness of the density matrix of semiconductors and insulators, and we scale our calculations to thousands of GPUs with a hierarchical loop- and data-distribution strategy. We demonstrate the efficacy of our method by computing the VEEs of several spin defects in wide band-gap materials, and we show that supercells with up to 1000 atoms are necessary to obtain converged results. We discuss the validity of common approximations, e.g., the solution of the BSE with truncated sums over empty and occupied states. We then apply our GW-BSE implementation to a diamond lattice with 1727 atoms to study the symmetry breaking of triplet states caused by the interaction of a point defect with an extended line defect. |
