Grid-based lattice summation of electrostatic potentials by assembled rank-structured tensor approximation

Abstract

We introduce and study a novel tensor approach for fast and accurate assembled summation of a large number of lattice-allocated potentials represented on 3D $N\times N\times N$ grid with the computational requirements only \emph{weakly dependent} on the number of summed potentials. It is based on the assembled low-rank canonical tensor representations of the collected potentials using pointwise sums of shifted canonical vectors representing the single generating function, say the Newton kernel. For a sum of electrostatic potentials over $L\times L\times L$ lattice embedded in a box the required storage scales linearly in the 1D grid-size, $O(N)$, while the numerical cost is estimated by $O(NL)$. For periodic boundary conditions, the storage demand remains proportional to the 1D grid-size of a unit cell, $n=N/L$, while the numerical cost reduces to $O(N)$, that outperforms the FFT-based Ewald-type summation algorithms of complexity $O(N^3\log N)$. The complexity in the grid parameter $N$ can be reduced even to the logarithmic scale $O(\log N)$ by using data-sparse representation of canonical $N$-vectors via the quantics tensor approximation. For justification, we prove an upper bound on the quantics ranks for the canonical vectors in the overall lattice sum. The presented approach is beneficial in applications which require further functional calculus with the lattice potential, say, scalar product with a function, integration or differentiation, which can be performed easily in tensor arithmetics on large 3D grids with 1D cost. Numerical tests illustrate the performance of the tensor summation method and confirm the estimated bounds on the tensor ranks.

This paper is an essentially improved version of the preprint of the Max-Planck Institute for Mathematics in the Sciences [48]. arXiv admin note: text overlap with arXiv:1411.1994

References in corpus (1)

• DMRG and periodic boundary conditions: a quantum information perspective

Cited by in corpus (13)

• Tensor Numerical Methods in Quantum Chemistry: from Hartree-Fock Energy to Excited States
• Fast iterative solution of the Bethe-Salpeter eigenvalue problem using low-rank and QTT tensor approximation
• Computation of the Response Surface in the Tensor Train data format
• Tensor Numerical Methods for High-dimensional PDEs: Basic Theory and Initial Applications
• Prospects of tensor-based numerical modeling of the collective electrostatic potential in many-particle systems
• Computing electrostatic potentials using regularization based on the range-separated tensor format
• Range-separated tensor representation of the discretized multidimensional Dirac delta and elliptic operator inverse
• Tucker tensor method for fast grid-based summation of long-range potentials on 3D lattices with defects
• Polynomial Chaos Expansion of random coefficients and the solution of stochastic partial differential equations in the Tensor Train format
• Computing the density of states for optical spectra by low-rank and QTT tensor approximation
• Tensor Numerical Approach to Linearized Hartree-Fock Equation for Lattice-type and Periodic Systems
• Block circulant and Toeplitz structures in the linearized Hartree-Fock equation on finite lattices: tensor approach
• Rank structured approximation method for quasi--periodic elliptic problems