Block circulant and Toeplitz structures in the linearized Hartree-Fock equation on finite lattices: tensor approach

Abstract

This paper introduces and analyses the new grid-based tensor approach to approximate solution of the elliptic eigenvalue problem for the 3D lattice-structured systems. We consider the linearized Hartree-Fock equation over a spatial $L_1\times L_2\times L_3$ lattice for both periodic and non-periodic problem setting, discretized in the localized Gaussian-type orbitals basis. In the periodic case, the Galerkin system matrix obeys a three-level block-circulant structure that allows the FFT-based diagonalization, while for the finite extended systems in a box (Dirichlet boundary conditions) we arrive at the perturbed block-Toeplitz representation providing fast matrix-vector multiplication and low storage size. The proposed grid-based tensor techniques manifest the twofold benefits: (a) the entries of the Fock matrix are computed by 1D operations using low-rank tensors represented on a 3D grid, (b) in the periodic case the low-rank tensor structure in the diagonal blocks of the Fock matrix in the Fourier space reduces the conventional 3D FFT to the product of 1D FFTs. Lattice type systems in a box with Dirichlet boundary conditions are treated numerically by our previous tensor solver for single molecules, which makes possible calculations on rather large $L_1\times L_2\times L_3$ lattices due to reduced numerical cost for 3D problems. The numerical simulations for both box-type and periodic $L\times 1\times 1$ lattice chain in a 3D rectangular "tube" with $L$ up to several hundred confirm the theoretical complexity bounds for the block-structured eigenvalue solvers in the limit of large $L$.

30 pages, 12 figures. arXiv admin note: substantial text overlap with arXiv:1408.3839

References in corpus (5)

• A new approach to the modelling of local defects in crystals: the reduced Hartree-Fock case
• Tensor Numerical Methods in Quantum Chemistry: from Hartree-Fock Energy to Excited States
• A reduced basis approach for calculation of the Bethe-Salpeter excitation energies using low-rank tensor factorizations
• Fast iterative solution of the Bethe-Salpeter eigenvalue problem using low-rank and QTT tensor approximation
• Grid-based electronic structure calculations: the tensor decomposition approach