skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Systematically Improvable Tensor Hypercontraction: Interpolative Separable Density-Fitting for Molecules Applied to Exact Exchange, Second- and Third-Order Møller–Plesset Perturbation Theory

Abstract

We present a systematically improvable tensor hypercontraction (THC) factorization based on interpolative separable density fitting (ISDF). We illustrate algorithmic details to achieve this within the framework of Becke's atom-centered quadrature grid. A single ISDF parameter cISDF controls the trade-off between accuracy and cost. In particular, cISDF sets the number of interpolation points used in THC, NIP = cISDF × NX with NX being the number of auxiliary basis functions. In conjunction with the resolution-of-the-identity (RI) technique, we develop and investigate the THC-RI algorithms for cubic-scaling exact exchange for Hartree-Fock and range-separated hybrids (e.g., ωB97X-V) and quartic-scaling second- and third-order Møller-Plesset theory (MP2 and MP3). These algorithms were evaluated over the W4-11 thermochemistry (atomization energy) set and A24 noncovalent interaction benchmark set with standard Dunning basis sets (cc-pVDZ, cc-pVTZ, aug-cc-pVDZ, and aug-cc-pVTZ). We demonstrate the convergence of THC-RI algorithms to numerically exact RI results using ISDF points. Based on these, we make recommendations on cISDF for each basis set and method. We also demonstrate the utility of THC-RI exact exchange and MP2 for larger systems such as water clusters and C20. We stress that more challenges await in obtaining accurate and numerically stable THC factorization for wave function amplitudes as wellmore » as for the space spanned by virtual orbitals in large basis sets and implementing sparsity-aware THC-RI algorithms.« less

Authors:
ORCiD logo [1]; ORCiD logo [2]; ORCiD logo [1]
  1. Univ. of California, Berkeley, CA (United States). Dept. of Chemistry; Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). Chemical Sciences Div.
  2. Univ. of California, Berkeley, CA (United States). Dept. of Mathematics; Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). Computational Research Div.
Publication Date:
Research Org.:
Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR). Scientific Discovery through Advanced Computing (SciDAC)
OSTI Identifier:
1779232
Grant/Contract Number:  
AC02-05CH11231
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Journal of Chemical Theory and Computation
Additional Journal Information:
Journal Volume: 16; Journal Issue: 1; Journal ID: ISSN 1549-9618
Publisher:
American Chemical Society
Country of Publication:
United States
Language:
English
Subject:
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY

Citation Formats

Lee, Joonho, Lin, Lin, and Head-Gordon, Martin. Systematically Improvable Tensor Hypercontraction: Interpolative Separable Density-Fitting for Molecules Applied to Exact Exchange, Second- and Third-Order Møller–Plesset Perturbation Theory. United States: N. p., 2019. Web. doi:10.1021/acs.jctc.9b00820.
Lee, Joonho, Lin, Lin, & Head-Gordon, Martin. Systematically Improvable Tensor Hypercontraction: Interpolative Separable Density-Fitting for Molecules Applied to Exact Exchange, Second- and Third-Order Møller–Plesset Perturbation Theory. United States. https://doi.org/10.1021/acs.jctc.9b00820
Lee, Joonho, Lin, Lin, and Head-Gordon, Martin. 2019. "Systematically Improvable Tensor Hypercontraction: Interpolative Separable Density-Fitting for Molecules Applied to Exact Exchange, Second- and Third-Order Møller–Plesset Perturbation Theory". United States. https://doi.org/10.1021/acs.jctc.9b00820. https://www.osti.gov/servlets/purl/1779232.
@article{osti_1779232,
title = {Systematically Improvable Tensor Hypercontraction: Interpolative Separable Density-Fitting for Molecules Applied to Exact Exchange, Second- and Third-Order Møller–Plesset Perturbation Theory},
author = {Lee, Joonho and Lin, Lin and Head-Gordon, Martin},
abstractNote = {We present a systematically improvable tensor hypercontraction (THC) factorization based on interpolative separable density fitting (ISDF). We illustrate algorithmic details to achieve this within the framework of Becke's atom-centered quadrature grid. A single ISDF parameter cISDF controls the trade-off between accuracy and cost. In particular, cISDF sets the number of interpolation points used in THC, NIP = cISDF × NX with NX being the number of auxiliary basis functions. In conjunction with the resolution-of-the-identity (RI) technique, we develop and investigate the THC-RI algorithms for cubic-scaling exact exchange for Hartree-Fock and range-separated hybrids (e.g., ωB97X-V) and quartic-scaling second- and third-order Møller-Plesset theory (MP2 and MP3). These algorithms were evaluated over the W4-11 thermochemistry (atomization energy) set and A24 noncovalent interaction benchmark set with standard Dunning basis sets (cc-pVDZ, cc-pVTZ, aug-cc-pVDZ, and aug-cc-pVTZ). We demonstrate the convergence of THC-RI algorithms to numerically exact RI results using ISDF points. Based on these, we make recommendations on cISDF for each basis set and method. We also demonstrate the utility of THC-RI exact exchange and MP2 for larger systems such as water clusters and C20. We stress that more challenges await in obtaining accurate and numerically stable THC factorization for wave function amplitudes as well as for the space spanned by virtual orbitals in large basis sets and implementing sparsity-aware THC-RI algorithms.},
doi = {10.1021/acs.jctc.9b00820},
url = {https://www.osti.gov/biblio/1779232}, journal = {Journal of Chemical Theory and Computation},
issn = {1549-9618},
number = 1,
volume = 16,
place = {United States},
year = {Tue Dec 03 00:00:00 EST 2019},
month = {Tue Dec 03 00:00:00 EST 2019}
}

Works referenced in this record:

Semi-exact concentric atomic density fitting: Reduced cost and increased accuracy compared to standard density fitting
journal, February 2014


Overcoming the Memory Bottleneck in Auxiliary Field Quantum Monte Carlo Simulations with Interpolative Separable Density Fitting
journal, November 2018


An implementation of RI-SCF on parallel computers
journal, January 1997


Tensor Hypercontraction Second-Order Møller–Plesset Perturbation Theory: Grid Optimization and Reaction Energies
journal, June 2015


Principles for a direct SCF approach to LICAO - MO ab - initio calculations : Direct
journal, September 1982


Pseudospectral Hartree–Fock theory: Applications and algorithmic improvements
journal, September 1990


Fitting the Coulomb potential variationally in Xα molecular calculations
journal, March 1983


W4-11: A high-confidence benchmark dataset for computational thermochemistry derived from first-principles W4 data
journal, July 2011


Auxiliary basis sets to approximate Coulomb potentials
journal, June 1995


Tensor decomposition in post-Hartree–Fock methods. II. CCD implementation
journal, December 2013


The prism algorithm for two-electron integrals
journal, December 1991


A method for two‐electron Gaussian integral and integral derivative evaluation using recurrence relations
journal, November 1988


Two-electron repulsion integrals over Gaussians functions
journal, December 1991


Fast, accurate evaluation of exact exchange: The occ-RI-K algorithm
journal, July 2015


Generalized Unitary Coupled Cluster Wave functions for Quantum Computation
journal, November 2018


Gaussian basis sets for use in correlated molecular calculations. III. The atoms aluminum through argon
journal, January 1993


Linear scaling density fitting
journal, November 2006


Solution of the Hartree–Fock equations by a pseudospectral method: Application to diatomic molecules
journal, August 1986


Laplace transform techniques in Mo/ller–Plesset perturbation theory
journal, January 1992


Pseudospectral full configuration interaction
journal, August 1992


Self‐consistent molecular orbital methods. XVIII. Constraints and stability in Hartree–Fock theory
journal, April 1977


What Are the Ground State Structures of C 20 and C 24 ? An Explicitly Correlated Ab Initio Approach
journal, December 2015


Improvements on the direct SCF method: Improved Direct SCF Method
journal, January 1989


Fast Hartree–Fock theory using local density fitting approximations
journal, November 2004


Advances in molecular quantum chemistry contained in the Q-Chem 4 program package
journal, September 2014


A direct atomic orbital driven implementation of the coupled cluster singles and doubles (CCSD) model
journal, September 1994


Linear-scaling atomic orbital-based second-order Møller–Plesset perturbation theory by rigorous integral screening criteria
journal, February 2009


New pseudospectral algorithms for electronic structure calculations: Length scale separation and analytical two‐electron integral corrections
journal, September 1994


Robust and variational fitting: Removing the four-center integrals from center stage in quantum chemistry
journal, September 2000


Linear scaling second-order Moller–Plesset theory in the atomic orbital basis for large molecular systems
journal, February 1999


Efficient use of the correlation consistent basis sets in resolution of the identity MP2 calculations
journal, February 2002


Linear scaling density functional calculations via the continuous fast multipole method
journal, May 1996


M�ller-Plesset (MP2) perturbation theory for large molecules
journal, November 1993


Cubic scaling algorithms for RPA correlation using interpolative separable density fitting
journal, December 2017


Linear Scaling Density Functional Calculations with Gaussian Orbitals
journal, June 1999


Large-scale correlated electronic structure calculations: the RI-MP2 method on parallel computers
journal, March 1996


Adaptively Compressed Polarizability Operator for Accelerating Large Scale Ab Initio Phonon Calculations
journal, January 2017


The continuous fast multipole method
journal, November 1994


Tensor-structured coupled cluster theory
journal, November 2017


Quartic scaling second-order approximate coupled cluster singles and doubles via tensor hypercontraction: THC-CC2
journal, March 2013


Computation of electron repulsion integrals involving contracted Gaussian basis functions
journal, May 1978


Scaled opposite-spin second order Møller–Plesset correlation energy: An economical electronic structure method
journal, November 2004


Making the random phase approximation to electronic correlation accurate
journal, October 2009


Efficient Implementation of the Pair Atomic Resolution of the Identity Approximation for Exact Exchange for Hybrid and Range-Separated Density Functionals
journal, January 2015


An efficient algorithm for the generation of two-electron repulsion integrals over gaussian basis functions
journal, April 1989


Tensor hypercontracted ppRPA: Reducing the cost of the particle-particle random phase approximation from O ( r6 ) to O ( r4 )
journal, July 2014


Tensor hypercontraction density fitting. I. Quartic scaling second- and third-order Møller-Plesset perturbation theory
journal, July 2012


Simplifications in the generation and transformation of two-electron integrals in molecular calculations
journal, October 1977


Low rank factorization of the Coulomb integrals for periodic coupled cluster theory
journal, March 2017


Pseudospectral Hartree–Fock gradient calculations
journal, June 1991


Tensor hypercontraction. II. Least-squares renormalization
journal, December 2012


On some approximations in applications of X α theory
journal, October 1979


Coupled-Cluster Valence-Bond Singles and Doubles for Strongly Correlated Systems: Block-Tensor Based Implementation and Application to Oligoacenes
journal, January 2017


Attractive electron-electron interactions within robust local fitting approximations
journal, April 2013


Communication: Acceleration of coupled cluster singles and doubles via orbital-weighted least-squares tensor hypercontraction
journal, May 2014


Coupled‐cluster theory employing approximate integrals: An approach to avoid the input/output and storage bottlenecks
journal, July 1994


Separable resolution-of-the-identity with all-electron Gaussian bases: Application to cubic-scaling RPA
journal, May 2019


A geometric approach to direct minimization
journal, June 2002


Reduced scaling CASPT2 using supporting subspaces and tensor hyper-contraction
journal, July 2018


Fast linear scaling second-order Møller-Plesset perturbation theory (MP2) using local and density fitting approximations
journal, May 2003


Electron affinities of the first‐row atoms revisited. Systematic basis sets and wave functions
journal, May 1992


Quartic scaling MP2 for solids: A highly parallelized algorithm in the plane wave basis
journal, March 2017


Coulombic potential energy integrals and approximations
journal, May 1973


Auxiliary basis sets for main row atoms and transition metals and their use to approximate Coulomb potentials
journal, October 1997

  • Eichkorn, Karin; Weigend, Florian; Treutler, Oliver
  • Theoretical Chemistry Accounts: Theory, Computation, and Modeling (Theoretica Chimica Acta), Vol. 97, Issue 1-4
  • https://doi.org/10.1007/s002140050244

Solution of the Hartree–Fock equations for polyatomic molecules by a pseudospectral method
journal, March 1987


Pseudospectral Mo/ller–Plesset perturbation theory through third order
journal, March 1994


Efficient recursive computation of molecular integrals over Cartesian Gaussian functions
journal, April 1986


On Resolution-of-the-Identity Electron Repulsion Integral Approximations and Variational Stability
journal, September 2017


Use of approximate integrals in ab initio theory. An application in MP2 energy calculations
journal, June 1993


Computation of electron repulsion integrals using the rys quadrature method
journal, July 1983


A multicenter numerical integration scheme for polyatomic molecules
journal, February 1988


Convergence acceleration of iterative sequences. the case of scf iteration
journal, July 1980


Pseudospectral double excitation configuration interaction
journal, May 1993


ImprovedSCF convergence acceleration
journal, January 1982


Discrete variable representation in electronic structure theory: Quadrature grids for least-squares tensor hypercontraction
journal, May 2013


Tensor Hypercontraction Equation-of-Motion Second-Order Approximate Coupled Cluster: Electronic Excitation Energies in O(N4) Time
journal, July 2013


Self-consistent molecular Hartree—Fock—Slater calculations I. The computational procedure
journal, September 1973


Linear and sublinear scaling formation of Hartree–Fock-type exchange matrices
journal, August 1998


Atomic orbital-based SOS-MP2 with tensor hypercontraction. II. Local tensor hypercontraction
journal, January 2017


Analytical gradients for tensor hyper-contracted MP2 and SOS-MP2 on graphical processing units
journal, October 2017


Electron repulsion integral approximations and error bounds: Molecular applications
journal, September 1974


Communication: Tensor hypercontraction. III. Least-squares tensor hypercontraction for the determination of correlated wavefunctions
journal, December 2012


A standard grid for density functional calculations
journal, July 1993


Integral approximations for LCAO-SCF calculations
journal, October 1993


Atomic orbital-based SOS-MP2 with tensor hypercontraction. I. GPU-based tensor construction and exploiting sparsity
journal, May 2016


Tensor Hypercontraction Equation-of-Motion Second-Order Approximate Coupled Cluster: Electronic Excitation Energies in O(N4) Time
journal, July 2013


Communication: Acceleration of coupled cluster singles and doubles via orbital-weighted least-squares tensor hypercontraction
journal, May 2014


Tensor-Structured Coupled Cluster Theory
text, January 2017