On the Best Rank-1 and Rank-(R1,R2,. . .,RN) Approximation of Higher-Order Tensors
SIAM Journal on Matrix Analysis and Applications
Rank-One Approximation to High Order Tensors
SIAM Journal on Matrix Analysis and Applications
Algorithms for Numerical Analysis in High Dimensions
SIAM Journal on Scientific Computing
Tensor-product approximation to operators and functions in high dimensions
Journal of Complexity
On tensor approximation of Green iterations for Kohn-Sham equations
Computing and Visualization in Science
Tucker Dimensionality Reduction of Three-Dimensional Arrays in Linear Time
SIAM Journal on Matrix Analysis and Applications
Multigrid Accelerated Tensor Approximation of Function Related Multidimensional Arrays
SIAM Journal on Scientific Computing
Journal of Computational and Applied Mathematics
Numerical Solution of the Hartree-Fock Equation in Multilevel Tensor-Structured Format
SIAM Journal on Scientific Computing
Wedderburn Rank Reduction and Krylov Subspace Method for Tensor Approximation. Part 1: Tucker Case
SIAM Journal on Scientific Computing
Hi-index | 31.45 |
In this paper, we investigate a novel approach based on the combination of Tucker-type and canonical tensor decomposition techniques for the efficient numerical approximation of functions and operators in electronic structure calculations. In particular, we study applicability of tensor approximations for the numerical solution of Hartree-Fock and Kohn-Sham equations on 3D Cartesian grids. We show that the orthogonal Tucker-type tensor approximation of electron density and Hartree potential of simple molecules leads to low tensor rank representations. This enables an efficient tensor-product convolution scheme for the computation of the Hartree potential using a collocation-type approximation via piecewise constant basis functions on a uniform nxnxn grid. Combined with the Richardson extrapolation, our approach exhibits O(h^3) convergence in the grid-size h=O(n^-^1). Moreover, this requires O(3rn+r^3) storage, where r denotes the Tucker rank of the electron density with r=O(logn), almost uniformly in n. For example, calculations of the Coulomb matrix and the Hartree-Fock energy for the CH"4 molecule, with a pseudopotential on the C atom, achieved accuracies of the order of 10^-^6 hartree with a grid-size n of several hundreds. Since the tensor-product convolution in 3D is performed via 1D convolution transforms, our scheme markedly outperforms the 3D-FFT in both the computing time and storage requirements.