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

Title: Highly Efficient and Scalable Compound Decomposition of Two-Electron Integral Tensor and Its Application in Coupled Cluster Calculations

Journal Article · · Journal of Chemical Theory and Computation
ORCiD logo [1]; ORCiD logo [1]
  1. William R. Wiley Environmental Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory, K8-91, P. O. Box 999, Richland, Washington 99352, United States

The representation and storage of two-electron integral tensors are vital in large- scale applications of accurate electronic structure methods. Low-rank representation and efficient storage strategy of integral tensors can significantly reduce the numerical overhead and consequently time-to-solution of these methods. In this paper, by combining pivoted incomplete Cholesky decomposition (CD) with a follow-up truncated singular vector decomposition (SVD), we develop a decomposition strategy to approximately represent the two-electron integral tensor in terms of low-rank vectors. A systematic benchmark test on a series of 1-D, 2-D, and 3-D carbon-hydrogen systems demonstrates high efficiency and scalability of the compound two-step decomposition of the two-electron integral tensor in our implementation. For the size of atomic basis set N_b ranging from ~ 100 up to ~ 2, 000, the observed numerical scaling of our implementation shows O(N_b^{2.5~3}) versus O(N_b^{3~4}) of single CD in most of other implementations. More importantly, this decomposition strategy can significantly reduce the storage requirement of the atomic-orbital (AO) two-electron integral tensor from O(N_b^4) to O(N_b^2 log_{10}(N_b)) with moderate decomposition thresholds. The accuracy tests have been performed using ground- and excited-state formulations of coupled- cluster formalism employing single and double excitations (CCSD) on several bench- mark systems including the C_{60} molecule described by nearly 1,400 basis functions. The results show that the decomposition thresholds can be generally set to 10^{-4} to 10^{-3} to give acceptable compromise between efficiency and accuracy.

Research Organization:
Pacific Northwest National Lab. (PNNL), Richland, WA (United States). Environmental Molecular Sciences Lab. (EMSL)
Sponsoring Organization:
USDOE
DOE Contract Number:
AC05-76RL01830
OSTI ID:
1414535
Report Number(s):
PNNL-SA-126829; 49220
Journal Information:
Journal of Chemical Theory and Computation, Vol. 13, Issue 9; ISSN 1549-9618
Publisher:
American Chemical Society
Country of Publication:
United States
Language:
English

Similar Records

Low-rank factorization of electron integral tensors and its application in electronic structure theory
Journal Article · Wed Mar 01 00:00:00 EST 2017 · Chemical Physics Letters · OSTI ID:1414535

Low-rank factorization of electron integral tensors and its application in electronic structure theory
Journal Article · Wed Jan 25 00:00:00 EST 2017 · Chemical Physics Letters · OSTI ID:1414535

Tensor decomposition in electronic structure calculations on 3D Cartesian grids
Journal Article · Tue Sep 01 00:00:00 EDT 2009 · Journal of Computational Physics · OSTI ID:1414535