Search

MiS Preprint Repository

We have decided to discontinue the publication of preprints on our preprint server as of 1 March 2024. The publication culture within mathematics has changed so much due to the rise of repositories such as ArXiV (www.arxiv.org) that we are encouraging all institute members to make their preprints available there. An institute's repository in its previous form is, therefore, unnecessary. The preprints published to date will remain available here, but we will not add any new preprints here.

MiS Preprint
12/2017

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

Venera Khoromskaia and Boris N. Khoromskij

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$.

Received:
Jan 30, 2017
Published:
Feb 2, 2017
MSC Codes:
65F30, 65F50, 65N35, 65F10
Keywords:
Tensor structured numerical methods for PDEs, 3D grid-based tensor approximation, Hartree-Fock equation, linearized Fock operator, periodic systems, block circulant/Toeplitz matrix, fast Fourier transform

Related publications

inJournal
2017 Repository Open Access
Venera Khoromskaia and Boris N. Khoromskij

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

In: Computational methods in applied mathematics, 17 (2017) 3, pp. 431-455