Fast iterative solution of the Bethe–Salpeter eigenvalue problem using low-rank and QTT tensor approximation

Peter Benner, Sergey Dolgov, Venera Khoromskaia, Boris N. Khoromskij

Research output: Contribution to journalArticle

11 Citations (Scopus)
127 Downloads (Pure)

Abstract

In this paper, we propose and study two approaches to approximate the solution of the Bethe–Salpeter equation (BSE) by using structured iterative eigenvalue solvers. Both approaches are based on the reduced basis method and low-rank factorizations of the generating matrices. We also propose to represent the static screen interaction part in the BSE matrix by a small active sub-block, with a size balancing the storage for rank-structured representations of other matrix blocks. We demonstrate by various numerical tests that the combination of the diagonal plus low-rank plus reduced-block approximation exhibits higher precision with low numerical cost, providing as well a distinct two-sided error estimate for the smallest eigenvalues of the Bethe–Salpeter operator. The complexity is reduced to O(Nb 2) in the size of the atomic orbitals basis set, Nb, instead of the practically intractable O(Nb 6) scaling for the direct diagonalization. In the second approach, we apply the quantized-TT (QTT) tensor representation to both, the long eigenvectors and the column vectors in the rank-structured BSE matrix blocks, and combine this with the ALS-type iteration in block QTT format. The QTT-rank of the matrix entities possesses almost the same magnitude as the number of occupied orbitals in the molecular systems, Nob, hence the overall asymptotic complexity for solving the BSE problem by the QTT approximation is estimated by O(log⁡(No)No 2). We confirm numerically a considerable decrease in computational time for the presented iterative approaches applied to various compact and chain-type molecules, while supporting sufficient accuracy.

Original languageEnglish
Pages (from-to)221-239
Number of pages19
JournalJournal of Computational Physics
Volume334
Early online date5 Jan 2017
DOIs
Publication statusPublished - 1 Apr 2017

Keywords

  • Bethe–Salpeter equation
  • Hartree–Fock calculus
  • Low-rank matrix
  • Model reduction
  • Quantized-TT format
  • Structured eigensolvers
  • Tensor decompositions

Fingerprint Dive into the research topics of 'Fast iterative solution of the Bethe–Salpeter eigenvalue problem using low-rank and QTT tensor approximation'. Together they form a unique fingerprint.

  • Projects

    Cite this