Loading…

Efficient and accurate algorithms for solving the Bethe–Salpeter eigenvalue problem for crystalline systems

Optical properties of materials related to light absorption and scattering are explained by the excitation of electrons. The Bethe–Salpeter equation is the state-of-the-art approach to describe these processes from first principles (ab initio), i.e. without the need for empirical data in the model....

Full description

Saved in:
Bibliographic Details
Published in:Journal of computational and applied mathematics 2022-01, Vol.400, p.113650, Article 113650
Main Authors: Benner, Peter, Penke, Carolin
Format: Article
Language:English
Subjects:
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Optical properties of materials related to light absorption and scattering are explained by the excitation of electrons. The Bethe–Salpeter equation is the state-of-the-art approach to describe these processes from first principles (ab initio), i.e. without the need for empirical data in the model. To harness the predictive power of the equation, it is mapped to an eigenvalue problem via an appropriate discretization scheme. The eigenpairs of the resulting large, dense, structured matrix can be used to compute dielectric properties of the considered crystalline or molecular system. The matrix always shows a 2 × 2 block structure. Depending on exact circumstances and discretization schemes, one ends up with a matrix structure such as H1=AB−B−A∈ℂ2n×2n,A=AH,B=BH,orH2=AB−BH−AT∈ℂ2n×2norR2n×2n,A=AH,B=BT.H1 can be acquired for crystalline systems (see Sander et al. (2015)), H2 is a more general form found e.g. in Shao et al. (2016) and Penke et al. (2020), which can for example be used to study molecules. Additionally, certain definiteness properties may hold. In this work, we compile theoretical results characterizing the structure of H1 and H2 in the language of non-standard scalar products. These results enable us to develop a generalized perspective on the currently used direct solution approach for matrices of form H1. This new viewpoint is used to develop two alternative methods for solving the eigenvalue problem. Both have advantages over the method currently in use and are well suited for high performance environments and only rely on basic numerical linear algebra building blocks. The results are extended to hold even without the mentioned definiteness property, showing the usefulness of our new perspective.
ISSN:0377-0427
1879-1778
DOI:10.1016/j.cam.2021.113650