-
Learning tensor networks with tensor cross interpolation: new algorithms and libraries
Authors:
Yuriel Núñez Fernández,
Marc K. Ritter,
Matthieu Jeannin,
Jheng-Wei Li,
Thomas Kloss,
Thibaud Louvet,
Satoshi Terasaki,
Olivier Parcollet,
Jan von Delft,
Hiroshi Shinaoka,
Xavier Waintal
Abstract:
The tensor cross interpolation (TCI) algorithm is a rank-revealing algorithm for decomposing low-rank, high-dimensional tensors into tensor trains/matrix product states (MPS). TCI learns a compact MPS representation of the entire object from a tiny training data set. Once obtained, the large existing MPS toolbox provides exponentially fast algorithms for performing a large set of operations. We di…
▽ More
The tensor cross interpolation (TCI) algorithm is a rank-revealing algorithm for decomposing low-rank, high-dimensional tensors into tensor trains/matrix product states (MPS). TCI learns a compact MPS representation of the entire object from a tiny training data set. Once obtained, the large existing MPS toolbox provides exponentially fast algorithms for performing a large set of operations. We discuss several improvements and variants of TCI. In particular, we show that replacing the cross interpolation by the partially rank-revealing LU decomposition yields a more stable and more flexible algorithm than the original algorithm. We also present two open source libraries, xfac in Python/C++ and TensorCrossInterpolation.jl in Julia, that implement these improved algorithms, and illustrate them on several applications. These include sign-problem-free integration in large dimension, the superhigh-resolution quantics representation of functions, the solution of partial differential equations, the superfast Fourier transform, the computation of partition functions, and the construction of matrix product operators.
△ Less
Submitted 4 July, 2024; v1 submitted 2 July, 2024;
originally announced July 2024.
-
Tensorized orbitals for computational chemistry
Authors:
Nicolas Jolly,
Yuriel Núñez Fernández,
Xavier Waintal
Abstract:
Choosing a basis set is the first step of a quantum chemistry calculation and it sets its maximum accuracy. This choice of orbitals is limited by strong technical constraints as one must be able to compute a large number of high dimensional integrals from these orbitals. Here we use tensor network techniques to construct representations of orbitals that essentially lift these technical constraints…
▽ More
Choosing a basis set is the first step of a quantum chemistry calculation and it sets its maximum accuracy. This choice of orbitals is limited by strong technical constraints as one must be able to compute a large number of high dimensional integrals from these orbitals. Here we use tensor network techniques to construct representations of orbitals that essentially lift these technical constraints. We show that a large class of orbitals can be put into ``tensorized'' form including the Gaussian orbitals, the exact hydrogenoid atomic orbitals and orbitals for chemical bonds. Our method provides a path for building more accurate and more compact basis sets beyond what has been accessible with previous technology.
△ Less
Submitted 7 August, 2023;
originally announced August 2023.
-
Quantics Tensor Cross Interpolation for High-Resolution, Parsimonious Representations of Multivariate Functions in Physics and Beyond
Authors:
Marc K. Ritter,
Yuriel Núñez Fernández,
Markus Wallerberger,
Jan von Delft,
Hiroshi Shinaoka,
Xavier Waintal
Abstract:
Multivariate functions of continuous variables arise in countless branches of science. Numerical computations with such functions typically involve a compromise between two contrary desiderata: accurate resolution of the functional dependence, versus parsimonious memory usage. Recently, two promising strategies have emerged for satisfying both requirements: (i) The quantics representation, which e…
▽ More
Multivariate functions of continuous variables arise in countless branches of science. Numerical computations with such functions typically involve a compromise between two contrary desiderata: accurate resolution of the functional dependence, versus parsimonious memory usage. Recently, two promising strategies have emerged for satisfying both requirements: (i) The quantics representation, which expresses functions as multi-index tensors, with each index representing one bit of a binary encoding of one of the variables; and (ii) tensor cross interpolation (TCI), which, if applicable, yields parsimonious interpolations for multi-index tensors. Here, we present a strategy, quantics TCI (QTCI), which combines the advantages of both schemes. We illustrate its potential with an application from condensed matter physics: the computation of Brillouin zone integrals.
△ Less
Submitted 25 March, 2024; v1 submitted 21 March, 2023;
originally announced March 2023.