Prospects of tensor-based numerical modeling of the collective electrostatic potential in many-particle systems

Abstract

Recently the rank-structured tensor approach suggested a progress in the numerical treatment of the long-range electrostatic potentials in many-particle systems and the respective interaction energy and forces [39,40,2]. In this paper, we outline the prospects for tensor-based numerical modeling of the collective electrostatic potential on lattices and in many-particle systems of general type. We generalize the approach initially introduced for the rank-structured grid-based calculation of the collective potentials on 3D lattices [39] to the case of many-particle systems with variable charges placed on $L^{\otimes d}$ lattices and discretized on fine $n^{\otimes d}$ Cartesian grids for arbitrary dimension $d$. As result, the interaction potential is represented in a parametric low-rank canonical format in $O(dLn)$ complexity. The energy is then calculated in $O(dL)$ operations. Electrostatics in large biomolecules is modeled by using the novel range-separated (RS) tensor format [2], which maintains the long-range part of the 3D collective potential of the many-body system represented on $n\times n\times n$ grid in a parametric low-rank form in $O(n)$-complexity. We show that the force field can be easily recovered by using the already precomputed electric field in the low-rank RS format. The RS tensor representation of the discretized Dirac delta [45] enables the efficient energy preserving regularization scheme for solving the 3D elliptic PDEs with strongly singular right-hand side arising in bio-sciences. We conclude that the rank-structured tensor-based approximation techniques provide the promising numerical tools for applications to many-body dynamics, protein docking and classification problems and for low-parametric interpolation of scattered data in data science.

arXiv comment30 pages, 23 figures

References in corpus (13)

• The density-matrix renormalization group in the age of matrix product states
• Low-Rank Tensor Networks for Dimensionality Reduction and Large-Scale Optimization Problems: Perspectives and Challenges PART 1
• Comments on PPPM, FMM, and the Ewald Method for Large Periodic Coulombic Systems
• How to mesh up Ewald sums (II): An accurate error estimate for the P3M algorithm
• Time integration of tensor trains
• On degeneration of tensors and algebras
• Grid-based lattice summation of electrostatic potentials by assembled rank-structured tensor approximation
• Grid-based electronic structure calculations: the tensor decomposition approach
• Computation of the Response Surface in the Tensor Train data format
• An inexact Newton-Krylov method for stochastic eigenvalue problems
• 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
• Tensor Method for Optimal Control Problems Constrained by Fractional 3D Elliptic Operator with Variable Coefficients

Cited by in corpus (1)

• Solving intractable chemical problems by tensor decomposition