Preprint 14/2016

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

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

Contact the author: Please use for correspondence this email.
Submission date: 09. Feb. 2016
published in: Journal of computational physics, 334 (2017), p. 221-239 
DOI number (of the published article): 10.1016/j.jcp.2016.12.047
Bibtex
MSC-Numbers: 65F30, 65F50, 65N35, 65F10
Keywords and phrases: Bethe-Salpeter equation, Hartree-Fock calculus, model reduction, quantized-TT format, tensor decompositions, structured eigensolvers, low-rank matrix

Abstract:
In this paper, we study and implement the structural iterative eigensolvers for the large-scale eigenvalue problem in the Bethe-Salpeter equation (BSE) based on the reduced basis approach via low-rank factorizations in generating matrices, introduced in the previous paper. The approach reduces numerical costs down to 𝒪(Nb2) in the size of atomic orbitals basis set, Nb, instead of practically intractable 𝒪(Nb6) complexity scaling for the direct diagonalization of the BSE matrix. As an alternative to rank approximation of the static screen interaction part of the BSE matrix, we propose to restrict it to a small active sub-block, with a size balancing the storage for rank-structured representations of other matrix blocks. We demonstrate that the enhanced reduced-block approximation exhibits higher precision within the controlled numerical cost, providing as well a distinct two-sided error estimate for the BSE eigenvalues. It is shown that further reduction of the asymptotic computational cost is possible due to ALS-type iteration in block tensor train (TT) format applied to the quantized-TT (QTT) tensor representation of both long eigenvectors and rank-structured matrix blocks. The QTT-rank of these entities possesses almost the same magnitude as the number of occupied orbitals in the molecular systems, No, hence the overall asymptotic complexity for solving the BSE problem can be estimated by 𝒪(log(No)No2). We confirm numerically a considerable decrease in computational time for the presented iterative approach applied to various compact and chain-type molecules, while supporting sufficient accuracy.

The manuscript is located here arxiv.org/abs/1602.02646

07.10.2017, 01:42