A reduced basis approach for calculation of the Bethe–Salpeter excitation energies by using low-rank tensor factorisations*. (17th April 2016)

• Record Type: Journal Article
• Authors: Benner, Peter; Khoromskaia, Venera; Khoromskij, Boris N.
• Is Part Of: Molecular physics. Volume 114:Number 7/8(2016)
• Page Start: 1148
• Page End: 1161
• Publication Date: 2016-04-17
• Subjects: Bethe–Salpeter equation -- Hartree–Fock equation -- two-electron integrals -- low-rank tensor decompositions -- model reduction -- reduced basis -- truncated Cholesky factorisation; 65F30 -- 65F50 -- 65N35 -- 65F10; Molecules -- Periodicals; Chemistry, Physical and theoretical -- Periodicals; Molécules -- Périodiques; Chimie physique et théorique -- Périodiques; 539.6.05
• DOI: 10.1080/00268976.2016.1149241
• Languages: English
• ISSNs: 0026-8976

Abstract

ABSTRACT: The Bethe–Salpeter equation (BSE) is a reliable model for estimating the absorption spectra in molecules and solids on the basis of accurate calculation of the excited states from first principles. Direct diagonalisation of the BSE matrix is practically intractable due to O ( N 6 ) complexity scaling in the size of the atomic orbital basis set, N . In this paper, we introduce and analyse a reduced basis approach to computation of the Bethe–Salpeter excitation energies which can lead to a relaxation of the numerical costs down to O ( N 3 ). The BSE operator is specified in terms of the two-electron integrals in the Hartree–Fock molecular orbital basis and the respective energies, calculated by the tensor-based solver described in previous works. The reduced basis method includes two steps. First, the diagonal plus low-rank approximation to fully populated blocks in the BSE matrix is calculated, enabling an easier partial eigenvalue solver for a large auxiliary system relying only on matrix–vector multiplications with rank-structured matrices. Second, a small subset of eigenvectors from the auxiliary eigenvalue problem is selected to build a projection of the exact BSE system onto this reduced basis set. Numerical tests demonstrate the ϵ-rank bounds for the blocks of the BSE matrix on examples of some compact molecules. The accuracy of the reduced basis approach vs. the effective matrix rank is illustrated. Abstract :