Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1504.06289 · Physical Chemistry Chemical Physics

Tensor numerical methods in quantum chemistry: from Hartree–Fock to excitation energies

Venera Khoromskaia, Boris N. Khoromskij

Published Jan 1, 2015• 68 citations
arXivHTMLPDF

Abstract

We resume the recent successes of the grid-based tensor numerical methods and discuss their prospects in real-space electronic structure calculations. These methods, based on the low-rank representation of the multidimensional functions and integral operators, first appeared as an accurate tensor calculus for the 3D Hartree potential using 1D complexity operations, and have evolved to entirely grid-based tensor-structured 3D Hartree-Fock eigenvalue solver. It benefits from tensor calculation of the core Hamiltonian and two-electron integrals (TEI) in O(n log n) complexity using the rank-structured approximation of basis functions, electron densities and convolution integral operators all represented on 3D n × n × n Cartesian grids. The algorithm for calculating TEI tensor in a form of the Cholesky decomposition is based on multiple factorizations using algebraic 1D "density fitting" scheme, which yield an almost irreducible number of product basis functions involved in the 3D convolution integrals, depending on a threshold ε > 0. The basis functions are not restricted to separable Gaussians, since the analytical integration is substituted by high-precision tensor-structured numerical quadratures. The tensor approaches to post-Hartree-Fock calculations for the MP2 energy correction and for the Bethe-Salpeter excitation energies, based on using low-rank factorizations and the reduced basis method, were recently introduced. Another direction is towards the tensor-based Hartree-Fock numerical scheme for finite lattices, where one of the numerical challenges is the summation of electrostatic potentials of a large number of nuclei. The 3D grid-based tensor method for calculation of a potential sum on a L × L × L lattice manifests the linear in L computational work, O(L), instead of the usual O(L(3) log L) scaling by the Ewald-type approaches.

DOI
10.1039/c5cp01215e

Related papers

Tensor Numerical Methods in Quantum Chemistry

Venera Khoromskaia, Boris N. Khoromskij

The conventional numerical methods when applied to multidimensional problems suffer from the so-called "curse of dimensionality", that cannot be eliminated by using parallel architectures and high performance computing. The novel tensor numerical methods are based on a "smart" rank-structured tensor representation of the multivariate functions and operators discretized on Cartesian grids thus reducing solution of the multidimensional integral-differential equations to 1D calculations. We explain basic tensor formats and algorithms and show how the orthogonal Tucker tensor decomposition originating from chemometrics made a revolution in numerical analysis, relying on rigorous results from approximation theory. Benefits of tensor approach are demonstrated in ab-initio electronic structure calculations. Computation of the 3D convolution integrals for functions with multiple singularities is replaced by a sequence of 1D operations, thus enabling accurate MATLAB calculations on a laptop using 3D uniform tensor grids of the size up to 10 15 . Fast tensor-based Hartree-Fock solver, incorporating the grid-based low-rank factorization of the two-electron integrals, serves as a prerequisite for economical calculation of the excitation energies of molecules. Tensor approach suggests efficient grid-based numerical treatment of the long-range electrostatic potentials on large 3D finite lattices with defects.The novel range-separated tensor format applies to interaction potentials of multi-particle systems of general type opening the new prospects for tensor methods in scientific computing. This research monograph presenting the modern tensor techniques applied to problems in quantum chemistry may be interesting for a wide audience of students and scientists working in computational chemistry, material science and scientific computing.

De Gruyter eBooks•Jun 8, 2018•45 citations

Tensor numerical methods for multidimensional PDES: theoretical analysis and initial applications

Boris N. Khoromskij

We present a brief survey on the modern tensor numerical methods for multidimensional stationary and time-dependent partial differential equations (PDEs). The guiding principle of the tensor approach is the rank-structured separable approximation of multivariate functions and operators represented on a grid. Recently, the traditional Tucker, canonical, and matrix product states (tensor train) tensor models have been applied to the grid-based electronic structure calculations, to parametric PDEs, and to dynamical equations arising in scientific computing. The essential progress is based on the quantics tensor approximation method proved to be capable to represent (approximate) function related d-dimensional data arrays of size Nd with log-volume complexity, O(dlog N). Combined with the traditional numerical schemes, these novel tools establish a new promising approach for solving multidimensional integral and differential equations using low-parametric rank-structured tensor formats. As the main example, we describe the grid-based tensor numerical approach for solving the 3D nonlinear Hartree-Fock eigenvalue problem, that was the starting point for the developments of tensor-structured numerical methods for large-scale computations in solving real-life multidimensional problems. We also discuss a new method for the fast 3D lattice summation of electrostatic potentials by assembled low-rank tensor approximation capable to treat the potential sum over millions of atoms in few seconds. We address new results on tensor approximation of the dynamical Fokker-Planck and master equations in many dimensions up to d = 20. Numerical tests demonstrate the benefits of the rank-structured tensor approximation on the aforementioned examples of multidimensional PDEs. In particular, the use of grid-based tensor representations in the reduced basis of atomics orbitals yields an accurate solution of the Hartree-Fock equation on large N × N × N grids with a grid size of up to N = 105.

ESAIM Proceedings and Surveys•Jan 1, 2015•58 citations

QTT Representation of the Hartree and Exchange Operators in Electronic Structure Calculations

Venera Khoromskaia, Boris N. Khoromskij, Reinhold Schneider

Abstract In this paper, the tensor-structured numerical evaluation of the Coulomb and exchange operators in the Hartree-Fock equation is supplemented by the usage of recent quantized-TT (QTT) formats. It leads to O(log n) complexity at computationally extensive stages in the rank-structured calculation with the respective 3D Hartree and exchange potentials discretized on large n×n×n Cartesian grids. The numerical examples for some volumetric organic molecules confirm that the QTT ranks of these potentials are nearly independent of the one-dimension grid size n. Thus, paradoxically, the complexity of the grid-based evaluation of the Coulumb and exchange matrices becomes almost independent of the grid size, being regulated only by the structure of a molecular system. As a result, the grid approximation of the Hartree-Fock equation allows to gain the high resolution with a guaranteed accuracy.

Computational Methods in Applied Mathematics•Jan 1, 2011•26 citations

RESEARCH ARTICLE Tensor Representation Techniques in post-Hartree Fock Methods: Matrix Product State Tensor Format

Udo Benedikt, Henry Auer, Mike Espig +2 more

;In this proof-of-principle study, we discuss the application of various tensor representation formats and their implications on memory requirements and computational eort for tensor manipulations as they occur in typical post Hartree-Fock methods. A successive tensor decomposition / rank reduction scheme in the matrix product state (MPS) format for the two electron integrals in the AO and MO basis and an estimate of the t2 amplitudes as obtained from MP2 are described. Furthermore, the AO-MO integral transformation, the calculation of the MP2 energy and the potential usage of tensors in low rank MPS representation for the tensor contractions in CC theory are discussed in detail. We are able to show that the overall scaling of the memory requirements is reduced from the conventional N 4 scaling to approximately N 3 and the scaling of computational eort for tensor contractions in post-HF methods can be reduced to roughly N4 while the decomposition itself scales as N5. While ecient algorithms with low prefactor for the tensor decomposition have yet to be devised, this Ansatz oers the possibility to nd a robust approximation with low scaling behaviour with system and basis set size for post Hartree Fock ab initio methods.

•Jan 1, 2013•0 citations

A Tensor-product approach for large scale electronic structure calculations using Kohn–Sham density functional theory

Phani Motamarri, Vikram Gavini, Thomas Blesgen

Quantum-mechanical calculations based on Kohn–Sham density functional theory (DFT) played a significant role in accurately predicting various aspects of materials behavior over the past decade. The Kohn–Sham approach to DFT reduces the many-body Schrodinger (eigen value) problem of interacting electrons into an equivalent problem of noninteracting electrons in an effective mean field that is governed by electron-density. Despite the reduced computational complexity of Kohn–Sham DFT, large-scale DFT calculations are still computationally very demanding with the resulting computational complexity scaling cubically with number of atoms in a given materials system. Numerical algorithms with reduced scaling behavior which are robust, computationally efficient and scalable on parallel computing architectures are always desirable to enable simulations at larger scales and on more complex systems. Following this line of thought, this study explores the use of tensor structured methods for ab-initio numerical solution of Kohn–Sham equations arising in DFT calculations. Earlier studies on tensor-structured methods have been quite successful in the accurate calculation of Hartree and the nonlocal exchange operators arising in the Hartree–Fock equations. A recent investigation of low-rank Tucker-type decomposition of the electron-density of large aluminum clusters (obtained from the finite-element discretization of orbital free DFT) shows the exponential decay of approximation error with respect to Tucker rank (number of tensor-basis functions in Tucker type representation). The results also indicate a smaller Tucker rank for the accurate representation of the electron density and is only weakly dependent on the system sizes studied. The promising success of tensor-structured techniques in resolving the electronic structure of material systems has enabled us to take a step further. In this study, we propose a systematic way of computing a globally adapted Tucker-type basis for solving the Kohn–Sham DFT problem by using a separable approximation of the Kohn–Sham Hamiltonian. Further, the resulting Kohn–Sham eigenvalue problem is projected into the aforementioned Tucker basis and is solved for ground-state energy using a self-consistent field iteration. The rank of the resulting Tucker representation and the computational complexity of the calculation are examined on representative benchmark examples involving metallic and insulating systems.

•Jan 1, 2014•0 citations

Tensor product approximation with optimal rank in quantum chemistry

Sambasiva Rao Chinnamsetty, Mike Espig, Boris N. Khoromskij +2 more

Tensor product decompositions with optimal separation rank provide an interesting alternative to traditional Gaussian-type basis functions in electronic structure calculations. We discuss various applications for a new compression algorithm, based on the Newton method, which provides for a given tensor the optimal tensor product or so-called best separable approximation for fixed Kronecker rank. In combination with a stable quadrature scheme for the Coulomb interaction, tensor product formats enable an efficient evaluation of Coulomb integrals. This is demonstrated by means of best separable approximations for the electron density and Hartree potential of small molecules, where individual components of the tensor product can be efficiently represented in a wavelet basis. We present a fairly detailed numerical analysis, which provides the basis for further improvements of this novel approach. Our results suggest a broad range of applications within density fitting schemes, which have been recently successfully applied in quantum chemistry.

The Journal of Chemical Physics•Aug 28, 2007•61 citations
We resume the recent successes of the grid-based tensor numerical methods and discuss their prospects in real-space electronic structure calculations. These methods, based on the low-rank representation of the multidimensional functions and integral operators, first appeared as an accurate tensor calculus for the 3D Hartree potential using 1D complexity operations, and have evolved to entirely grid-based tensor-structured 3D Hartree-Fock eigenvalue solver. It benefits from tensor calculation of the core Hamiltonian and two-electron integrals (TEI) in O(n log n) complexity using the rank-structured approximation of basis functions, electron densities and convolution integral operators all represented on 3D n × n × n Cartesian grids. The algorithm for calculating TEI tensor in a form of the Cholesky decomposition is based on multiple factorizations using algebraic 1D "density fitting" scheme, which yield an almost irreducible number of product basis functions involved in the 3D convolution integrals, depending on a threshold ε > 0. The basis functions are not restricted to separable Gaussians, since the analytical integration is substituted by high-precision tensor-structured numerical quadratures. The tensor approaches to post-Hartree-Fock calculations for the MP2 energy correction and for the Bethe-Salpeter excitation energies, based on using low-rank factorizations and the reduced basis method, were recently introduced. Another direction is towards the tensor-based Hartree-Fock numerical scheme for finite lattices, where one of the numerical challenges is the summation of electrostatic potentials of a large number of nuclei. The 3D grid-based tensor method for calculation of a potential sum on a L × L × L lattice manifests the linear in L computational work, O(L), instead of the usual O(L(3) log L) scaling by the Ewald-type approaches.
The conventional numerical methods when applied to multidimensional problems suffer from the so-called "curse of dimensionality", that cannot be eliminated by using parallel architectures and high performance computing. The novel tensor numerical methods are based on a "smart" rank-structured tensor representation of the multivariate functions and operators discretized on Cartesian grids thus reducing solution of the multidimensional integral-differential equations to 1D calculations. We explain basic tensor formats and algorithms and show how the orthogonal Tucker tensor decomposition originating from chemometrics made a revolution in numerical analysis, relying on rigorous results from approximation theory. Benefits of tensor approach are demonstrated in ab-initio electronic structure calculations. Computation of the 3D convolution integrals for functions with multiple singularities is replaced by a sequence of 1D operations, thus enabling accurate MATLAB calculations on a laptop using 3D uniform tensor grids of the size up to 10 15 . Fast tensor-based Hartree-Fock solver, incorporating the grid-based low-rank factorization of the two-electron integrals, serves as a prerequisite for economical calculation of the excitation energies of molecules. Tensor approach suggests efficient grid-based numerical treatment of the long-range electrostatic potentials on large 3D finite lattices with defects.The novel range-separated tensor format applies to interaction potentials of multi-particle systems of general type opening the new prospects for tensor methods in scientific computing. This research monograph presenting the modern tensor techniques applied to problems in quantum chemistry may be interesting for a wide audience of students and scientists working in computational chemistry, material science and scientific computing.
We present a brief survey on the modern tensor numerical methods for multidimensional stationary and time-dependent partial differential equations (PDEs). The guiding principle of the tensor approach is the rank-structured separable approximation of multivariate functions and operators represented on a grid. Recently, the traditional Tucker, canonical, and matrix product states (tensor train) tensor models have been applied to the grid-based electronic structure calculations, to parametric PDEs, and to dynamical equations arising in scientific computing. The essential progress is based on the quantics tensor approximation method proved to be capable to represent (approximate) function related d-dimensional data arrays of size Nd with log-volume complexity, O(dlog N). Combined with the traditional numerical schemes, these novel tools establish a new promising approach for solving multidimensional integral and differential equations using low-parametric rank-structured tensor formats. As the main example, we describe the grid-based tensor numerical approach for solving the 3D nonlinear Hartree-Fock eigenvalue problem, that was the starting point for the developments of tensor-structured numerical methods for large-scale computations in solving real-life multidimensional problems. We also discuss a new method for the fast 3D lattice summation of electrostatic potentials by assembled low-rank tensor approximation capable to treat the potential sum over millions of atoms in few seconds. We address new results on tensor approximation of the dynamical Fokker-Planck and master equations in many dimensions up to d = 20. Numerical tests demonstrate the benefits of the rank-structured tensor approximation on the aforementioned examples of multidimensional PDEs. In particular, the use of grid-based tensor representations in the reduced basis of atomics orbitals yields an accurate solution of the Hartree-Fock equation on large N × N × N grids with a grid size of up to N = 105.
Abstract In this paper, the tensor-structured numerical evaluation of the Coulomb and exchange operators in the Hartree-Fock equation is supplemented by the usage of recent quantized-TT (QTT) formats. It leads to O(log n) complexity at computationally extensive stages in the rank-structured calculation with the respective 3D Hartree and exchange potentials discretized on large n×n×n Cartesian grids. The numerical examples for some volumetric organic molecules confirm that the QTT ranks of these potentials are nearly independent of the one-dimension grid size n. Thus, paradoxically, the complexity of the grid-based evaluation of the Coulumb and exchange matrices becomes almost independent of the grid size, being regulated only by the structure of a molecular system. As a result, the grid approximation of the Hartree-Fock equation allows to gain the high resolution with a guaranteed accuracy.
;In this proof-of-principle study, we discuss the application of various tensor representation formats and their implications on memory requirements and computational eort for tensor manipulations as they occur in typical post Hartree-Fock methods. A successive tensor decomposition / rank reduction scheme in the matrix product state (MPS) format for the two electron integrals in the AO and MO basis and an estimate of the t2 amplitudes as obtained from MP2 are described. Furthermore, the AO-MO integral transformation, the calculation of the MP2 energy and the potential usage of tensors in low rank MPS representation for the tensor contractions in CC theory are discussed in detail. We are able to show that the overall scaling of the memory requirements is reduced from the conventional N 4 scaling to approximately N 3 and the scaling of computational eort for tensor contractions in post-HF methods can be reduced to roughly N4 while the decomposition itself scales as N5. While ecient algorithms with low prefactor for the tensor decomposition have yet to be devised, this Ansatz oers the possibility to nd a robust approximation with low scaling behaviour with system and basis set size for post Hartree Fock ab initio methods.
Quantum-mechanical calculations based on Kohn–Sham density functional theory (DFT) played a significant role in accurately predicting various aspects of materials behavior over the past decade. The Kohn–Sham approach to DFT reduces the many-body Schrodinger (eigen value) problem of interacting electrons into an equivalent problem of noninteracting electrons in an effective mean field that is governed by electron-density. Despite the reduced computational complexity of Kohn–Sham DFT, large-scale DFT calculations are still computationally very demanding with the resulting computational complexity scaling cubically with number of atoms in a given materials system. Numerical algorithms with reduced scaling behavior which are robust, computationally efficient and scalable on parallel computing architectures are always desirable to enable simulations at larger scales and on more complex systems. Following this line of thought, this study explores the use of tensor structured methods for ab-initio numerical solution of Kohn–Sham equations arising in DFT calculations. Earlier studies on tensor-structured methods have been quite successful in the accurate calculation of Hartree and the nonlocal exchange operators arising in the Hartree–Fock equations. A recent investigation of low-rank Tucker-type decomposition of the electron-density of large aluminum clusters (obtained from the finite-element discretization of orbital free DFT) shows the exponential decay of approximation error with respect to Tucker rank (number of tensor-basis functions in Tucker type representation). The results also indicate a smaller Tucker rank for the accurate representation of the electron density and is only weakly dependent on the system sizes studied. The promising success of tensor-structured techniques in resolving the electronic structure of material systems has enabled us to take a step further. In this study, we propose a systematic way of computing a globally adapted Tucker-type basis for solving the Kohn–Sham DFT problem by using a separable approximation of the Kohn–Sham Hamiltonian. Further, the resulting Kohn–Sham eigenvalue problem is projected into the aforementioned Tucker basis and is solved for ground-state energy using a self-consistent field iteration. The rank of the resulting Tucker representation and the computational complexity of the calculation are examined on representative benchmark examples involving metallic and insulating systems.
Tensor product decompositions with optimal separation rank provide an interesting alternative to traditional Gaussian-type basis functions in electronic structure calculations. We discuss various applications for a new compression algorithm, based on the Newton method, which provides for a given tensor the optimal tensor product or so-called best separable approximation for fixed Kronecker rank. In combination with a stable quadrature scheme for the Coulomb interaction, tensor product formats enable an efficient evaluation of Coulomb integrals. This is demonstrated by means of best separable approximations for the electron density and Hartree potential of small molecules, where individual components of the tensor product can be efficiently represented in a wavelet basis. We present a fairly detailed numerical analysis, which provides the basis for further improvements of this novel approach. Our results suggest a broad range of applications within density fitting schemes, which have been recently successfully applied in quantum chemistry.