Solvers for the Hermitian and the pseudo-Hermitian Bethe-Salpeter equation in the Yambo code: Implementation and Performance

Abstract

We analyze the performance of two strategies in solving the structured eigenvalue problem deriving from the Bethe-Salpeter equation (BSE) in condensed matter physics. The BSE matrix is constructed with the Yambo code, and the two strategies are implemented by interfacing Yambo with the ScaLAPACK and ELPA libraries for direct diagonalization, and with the SLEPc library for the iterative approach. We consider both the Hermitian (Tamm-Dancoff approximation) and pseudo-Hermitian forms, addressing dense matrices of three different sizes. A description of the implementation is also provided, with details for the pseudo-Hermitian case. Timing and memory utilization are analyzed on both CPU and GPU clusters. Our results demonstrate that it is now feasible to handle dense BSE matrices of the order of 10^5.

Submitted to SciPost Physics Codebases

References in corpus (7)

• GPAW: An open Python package for electronic-structure calculations
• Physical origin of giant excitonic and magneto-optical responses in two-dimensional ferromagnetic insulators
• Speeding up the solution of the Bethe-Salpeter equation by a double-grid method and Wannier interpolation
• Exciton lifetime and optical linewidth profile via exciton-phonon interactions: Theory and first-principles calculations for monolayer MoS${}_2$
• SIMPLE code: optical properties with optimal basis functions
• Excitonic response in TMD heterostructures from first-principles: impact of stacking, twisting, and interlayer distance
• Excitons under strain: light absorption and emission in strained hexagonal boron nitride