Projects per year
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 language | English |
---|---|
Pages (from-to) | 221-239 |
Number of pages | 19 |
Journal | Journal of Computational Physics |
Volume | 334 |
Early online date | 5 Jan 2017 |
DOIs | |
Publication status | Published - 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
- 2 Finished
-
Tensor product numerical methods for high-dimensional problems in probability and quantum calculations
Dolgov, S. (PI)
1/01/16 → 31/12/18
Project: Research council
-
Sergey Dolgov Fellowship - Tensor Product Numerical Methods for High-Dimensional Problems in Probablility and Quantum Calculations
Scheichl, R. (PI)
Engineering and Physical Sciences Research Council
1/01/16 → 31/12/18
Project: Research council