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
62/2014

Tensor Approach to Linearized Hartree-Fock Equation for Lattice-type and Periodic Systems

Venera Khoromskaia and Boris N. Khoromskij

Abstract

This paper introduces and analyses the new grid-based tensor approach for approximate solution of the eigenvalue problem for linearized Hartree-Fock equation applied to the 3D lattice-structured and periodic systems.

The set of localized basis functions over spatial $(L_1,L_2,L_3)$ lattice in a bounding box (or supercell) is assembled by multiple replicas of those from the unit cell. All basis functions and operators are discretized on a global 3D tensor grid in the bounding box which enables rather general basis sets.

In the periodic case, the Galerkin Fock matrix is shown to have the three-level block circulant structure, that allows the FFT-based diagonalization. The proposed 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) 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.

We describe fast numerical algorithms for the block circulant representation of the core Hamiltonian in the periodic setting based on low-rank tensor representation of arising multidimensional functions.

Lattice type systems in a box with open boundary conditions are treated by our previous tensor solver for single molecules, which makes possible calculations on large $(L_1,L_2,L_3)$ lattices due to reduced numerical cost for 3D problems. The numerical simulations for box/periodic $(L,1,1)$ lattice systems in a 3D rectangular "tube" with $L$ up to several hundred confirm the theoretical complexity bounds for the tensor-structured eigenvalue solvers in the limit of large $L$.

Received:
Jun 19, 2014
Published:
Jun 25, 2014
MSC Codes:
65F30, 65F50, 65N35, 65F10
Keywords:
Hartree-Fock equation, tensor-structured numerical methods, 3D grid-based tensor approximation, Fock operator, Fourier transform, periodic systems, lattice summation, block circulant matrix

Related publications

Preprint
2014 Repository Open Access
Venera Khoromskaia and Boris N. Khoromskij

Tensor numerical approach to linearized Hartree-Fock equation for lattice-type and periodic systems