GPU-Accelerated Solution of the Bethe-Salpeter Equation for Large and Heterogeneous Systems

Abstract

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 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. Further computational savings are achieved by exploiting the nearsightedness of the density matrix of semiconductors and insulators to reduce the number of screened Coulomb integrals. We scale our calculations to thousands of GPUs with a hierarchical loop and data distribution strategy. The efficacy of our method is demonstrated by computing the VEEs of several spin defects in wide-band-gap materials, showing that supercells with up to 1000 atoms are necessary to obtain converged results. We discuss the validity of the common approximation that solves 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.

References in corpus (46)

• CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
• Optimized norm-conserving Vanderbilt pseudopotentials
• The nitrogen-vacancy colour centre in diamond
• Quantum ESPRESSO toward the exascale
• Optimization Algorithm for the Generation of ONCV Pseudopotentials
• Sensitivity Optimization for NV-Diamond Magnetometry
• BerkeleyGW: A Massively Parallel Computer Package for the Calculation of the Quasiparticle and Optical Properties of Materials and Nanostructures
• Ab initio calculation of excitonic effects in the optical spectra of semiconductors
• NWChem: Past, Present, and Future
• SIESTA: recent developments and applications
• Many-body perturbation theory calculations using the yambo code
• Properties of nitrogen-vacancy centers in diamond: group theoretic approach
• First-principles GW calculations for fullerenes, porphyrins, phtalocyanine, and other molecules of interest for organic photovoltaic applications
• Self-consistent hybrid functional for condensed systems
• Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems
• Polytype control of spin qubits in silicon carbide
• Turbo charging time-dependent density-functional theory with Lanczos chains
• Observation of an environmentally insensitive solid state spin defect in diamond
• The negatively charged nitrogen-vacancy centre in diamond: the electronic solution
• A new and efficient approach to time-dependent density-functional perturbation theory for optical spectroscopy
• Quasiparticle GW calculations for solids, molecules and 2D materials
• Optical excitations in organic molecules, clusters and defects studied by first-principles Green's function methods
• Phonon-Induced Population Dynamics and Intersystem Crossing in Nitrogen-Vacancy Centers
• New infrared emission of the NV centre in diamond: Zeeman and uniaxial stress studies
• Isolated spin qubits in SiC with a high-fidelity infrared spin-to-photon interface
• Quantum simulations of materials on near-term quantum computers
• Nonempirical Range-separated Hybrid Functionals for Solids and Molecules
• State-selective intersystem crossing in nitrogen-vacancy centers
• \emph{Ab initio} calculation of spin-orbit coupling for NV center in diamond exhibiting dynamic Jahn-Teller effect
• Theory of the optical spinpolarization loop of the nitrogen-vacancy center in diamond
• Green's function formulation of quantum defect embedding theory
• Roadmap on Electronic Structure Codes in the Exascale Era
• Ab Initio Bethe-Salpeter Equation Approach to Neutral Excitations in Molecules with Numeric Atom-Centered Orbitals
• A finite field approach to solving the Bethe Salpeter equation
• Electronic structure of the neutral silicon-vacancy center in diamond
• Large scale GW-BSE calculations with N3 scaling: excitonic effects in dye sensitised solar cells
• Excited state properties of point defects in semiconductors and insulators investigated with time-dependent density functional theory
• GPGPU Acceleration of All-Electron Electronic Structure Theory Using Localized Numeric Atom-Centered Basis Functions
• First-principles Studies of Strongly Correlated States in Defect Spin Qubits in Diamond
• Koopmans Meets Bethe-Salpeter: Excitonic Optical Spectra without GW
• Theoretical polarization of zero phonon lines in point defects
• A Finite-field Approach for $GW$ Calculations Beyond the Random Phase Approximation
• Stability and electronic structure of NV centers at dislocation cores in diamond
• Computational protocol to evaluate electron-phonon interactions within density matrix perturbation theory
• Quantum Vibronic Effects on the Excitation Energies of the Nitrogen-Vacancy Center in Diamond
• Solving the Bethe-Salpeter equation on massively parallel architectures

Cited by in corpus (4)

• First-Principles Framework for the Prediction of Intersystem Crossing Rates in Spin Defects: The Role of Electron Correlation
• A Guide to Molecular Properties from the Bethe-Salpeter Equation
• Towards dislocation-driven quantum interconnects
• Many-body perturbation theory with hybrid density functional theory starting points accelerated by adaptively compressed exchange