Structure Preserving Parallel Algorithms for Solving the Bethe-Salpeter Eigenvalue Problem

Abstract

The Bethe-Salpeter eigenvalue problem is a dense structured eigenvalue problem arising from discretized Bethe-Salpeter equation in the context of computing exciton energies and states. A computational challenge is that at least half of the eigenvalues and the associated eigenvectors are desired in practice. We establish the equivalence between Bethe-Salpeter eigenvalue problems and real Hamiltonian eigenvalue problems. Based on theoretical analysis, structure preserving algorithms for a class of Bethe-Salpeter eigenvalue problems are proposed. We also show that for this class of problems all eigenvalues obtained from the Tamm-Dancoff approximation are overestimated. In order to solve large scale problems of practical interest, we discuss parallel implementations of our algorithms targeting distributed memory systems. Several numerical examples are presented to demonstrate the efficiency and accuracy of our algorithms.

Cited by in corpus (15)

• Many-body theory of positron binding in polyatomic molecules
• ELSI -- An Open Infrastructure for Electronic Structure Solvers
• Fully relativistic $GW$/Bethe-Salpeter calculations in BerkeleyGW: implementation, symmetries, benchmarking, and performance
• A reduced basis approach for calculation of the Bethe-Salpeter excitation energies using low-rank tensor factorizations
• Fast iterative solution of the Bethe-Salpeter eigenvalue problem using low-rank and QTT tensor approximation
• A structure preserving Lanczos algorithm for computing the optical absorption spectrum
• A Model Order Reduction Algorithm for Estimating the Absorption Spectrum
• High Performance Solution of Skew-symmetric Eigenvalue Problems with Applications in Solving the Bethe-Salpeter Eigenvalue Problem
• Energy-Specific Bethe-Salpeter Equation Implementation for Efficient Optical Spectrum Calculations
• Computing the density of states for optical spectra by low-rank and QTT tensor approximation
• Structured eigenvalue problems in electronic structure methods from a unified perspective
• Leveraging configuration interaction singles for qualitative descriptions of ground and excited states: state-averaging, linear-response, and spin-projection
• Study of Time Evolution for Approximation of Two-Body Spinless-Salpeter Equation
• Doubling algorithm for the discretized Bethe-Salpeter eigenvalue problem
• Solvers for the Hermitian and the pseudo-Hermitian Bethe-Salpeter equation in the Yambo code: Implementation and Performance