|
This article is cited in 3 scientific papers (total in 3 papers)
Mathematical physics
Prospects of tensor-based numerical modeling of the collective electrostatics in many-particle systems
V. Kh. Khoromskaiaa, B. N. Khoromskyab a D-04103 Leipzig, Inselstr. 22–26, Max Planck Institute for Mathematics in the Sciences, Germany
b Magdeburg, Max Planck Institute for Dynamics of Complex Technical Systems, Germany
Abstract:
Recently the rank-structured tensor approach suggested a progress in the numerical treatment of the long-range electrostatics in many-particle systems and the respective interaction energy and forces. 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. Our approach, initially introduced for the rank-structured grid-based calculation of the interaction potentials on 3D lattices is generalized here 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 a result, the interaction potential is represented in a parametric low-rank canonical format in $O(dLn)$ complexity. The total interaction energy can be then calculated in $O(dL)$ operations. Electrostatics in large bio-molecular systems is discretized on a fine $n^{\otimes 3}$ grid by using the novel range-separated $(\mathrm{RS})$ tensor format, which maintains the long-range part of the 3D collective potential of a many-body system in a parametric low-rank form in $O(n)$-complexity. We show how the energy and force field can be easily recovered by using the already precomputed electric field in the low-rank $\mathrm{RS}$ format. The $\mathrm{RS}$ tensor representation of the discretized Dirac delta enables the construction of the efficient energy preserving (conservative) regularization scheme for solving the 3D elliptic partial differential equations with strongly singular right-hand side arising in scientific computing. We conclude that the rank-structured tensor-based approximation techniques provide the promising numerical tools for applications to many-body dynamics in bio-sciences, protein docking and classification problems, for low-parametric interpolation of scattered data in data science, as well as in machine learning in many dimensions.
Key words:
Coulomb potential, Slater potential, long-range many-particle interactions, low-rank tensor decomposition, range-separated tensor formats, summation of electrostatic potentials, energy and force calculations.
Received: 24.12.2020 Revised: 24.12.2020 Accepted: 14.01.2021
Citation:
V. Kh. Khoromskaia, B. N. Khoromsky, “Prospects of tensor-based numerical modeling of the collective electrostatics in many-particle systems”, Zh. Vychisl. Mat. Mat. Fiz., 61:5 (2021), 895; Comput. Math. Math. Phys., 61:5 (2021), 864–886
Linking options:
https://www.mathnet.ru/eng/zvmmf11246 https://www.mathnet.ru/eng/zvmmf/v61/i5/p895
|
|