Fast tensor method for summation of long-range potentials on 3D lattices with defects.

DOI10.1002/NLA.2023zbMATH Open1413.65141arXiv1411.1994OpenAlexW2148442881MaRDI QIDQ5739760FDOQ5739760

Authors: Boris N. Khoromskij, V. Khoromskaia 

Publication date: 19 July 2016

Published in: Numerical Linear Algebra with Applications (Search for Journal in Brave)

Abstract: In this paper, we present a method for fast summation of long-range potentials on 3D lattices with multiple defects and having non-rectangular geometries, based on rank-structured tensor representations. This is a significant generalization of our recent technique for the grid-based summation of electrostatic potentials on the rectangular $LimesLimesL$ lattices by using the canonical tensor decompositions and yielding the $O(L)$ computational complexity instead of $O(L^3)$ by traditional approaches. The resulting lattice sum is calculated as a Tucker or canonical representation whose directional vectors are assembled by the 1D summation of the generating vectors for the shifted reference tensor, once precomputed on large $NimesNimesN$ representation grid in a 3D bounding box. The tensor numerical treatment of defects is performed in an algebraic way by simple summation of tensors in the canonical or Tucker formats. To diminish the considerable increase in the tensor rank of the resulting potential sum the $varepsilon$-rank reduction procedure is applied based on the generalized reduced higher-order SVD scheme. For the reduced higher-order SVD approximation to a sum of canonical/Tucker tensors, we prove the stable error bounds in the relative norm in terms of discarded singular values of the side matrices. The required storage scales linearly in the 1D grid-size, $O(N)$, while the numerical cost is estimated by $O(NL)$. Numerical tests confirm the efficiency of the presented tensor summation method: we demonstrate that a sum of millions of Newton kernels on a 3D lattice with defects/impurities can be computed in seconds in Matlab implementation.

Full work available at URL: https://arxiv.org/abs/1411.1994

Recommendations

• Grid-based lattice summation of electrostatic potentials by assembled rank-structured tensor approximation
• Prospects of tensor-based numerical modeling of the collective electrostatics in many-particle systems
• Range-separated tensor format for many-particle modeling
• Tensor numerical methods in quantum chemistry
• Tensor numerical methods for multidimensional PDEs: theoretical analysis and initial applications

zbMATH Keywords

lattice sumselectronic structure calculationstensor numerical methodscanonical and Tucker tensor formatsreduced higher-order SVDdefected lattice in a boxlong-range interaction potentials

Mathematics Subject Classification ID

Multilinear algebra, tensor calculus (15A69) Numerical computation of eigenvalues and eigenvectors of matrices (65F15) Lattice systems (Ising, dimer, Potts, etc.) and systems on graphs arising in equilibrium statistical mechanics (82B20)

Cites Work

• Orthogonal tensor decompositions
• A fast algorithm for particle simulations
• Tensor Decompositions and Applications
• Tensor-train decomposition
• A Multilinear Singular Value Decomposition
• Title not available (Why is that?)
• A literature survey of low-rank tensor approximation techniques
• The density-matrix renormalization group in the age of matrix product states
• Low-rank quadrature-based tensor approximation of the Galerkin projected Newton/Yukawa kernels
• Low rank Tucker-type tensor approximation to classical potentials
• Multigrid accelerated tensor approximation of function related multidimensional arrays
• Formulation and optimization of the energy-based blended quasicontinuum method
• Fast and accurate 3D tensor calculation of the Fock operator in a general basis
• Structured Rank-(r1, . . . , rd) Decomposition of Function-related Tensors in R_D
• \(O(d \log N)\)-quantics approximation of \(N\)-\(d\) tensors in high-dimensional numerical modeling
• Low-rank Kronecker-product approximation to multi-dimensional nonlocal operators I. Separable approximation of multi-variate functions
• Mathematical modeling of point defects in materials science
• Data-sparse approximation to a class of operator-valued functions
• Comments on P\(^3\)M, FMM, and the Ewald method for large periodic Coulombic systems
• Periodic Schrödinger operators with local defects and spectral pollution
• Ewald summation techniques in perspective: A survey
• Asymptotics for the approximation of wave functions by exponential sums
• Grid-based lattice summation of electrostatic potentials by assembled rank-structured tensor approximation
• Black-box Hartree-Fock solver by tensor numerical methods

Cited In (10)

• Range-separated tensor decomposition of the discretized Dirac delta and elliptic operator inverse
• Low-rank tensor methods for partial differential equations
• Tucker tensor analysis of Matérn functions in spatial statistics
• Prospects of tensor-based numerical modeling of the collective electrostatics in many-particle systems
• NFFT based Ewald summation for electrostatic systems with charges and dipoles
• On the efficient computation of large scale singular sums with applications to long-range forces in crystal lattices
• Block circulant and Toeplitz structures in the linearized Hartree-Fock equation on finite lattices: tensor approach
• Range-separated tensor format for many-particle modeling
• Mathematical methods in quantum chemistry. Abstracts from the workshop held March 18--24, 2018
• Grid-based lattice summation of electrostatic potentials by assembled rank-structured tensor approximation

Uses Software

• Matlab
• CRYSCOR

This page was built for publication: Fast tensor method for summation of long-range potentials on 3D lattices with defects.

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5739760)