A critical analysis of least-squares tensor hypercontraction applied to MP3.

The least-squares tensor hypercontraction (LS-THC) approach is a promising method of reducing the high polynomial scaling of wavefunction methods, for example, those based on many-body perturbation theory or coupled cluster. Here, we focus on LS-THC-MP3 and identify four variants with differing errors and efficiency characteristics. The performance of LS-THC-MP3 is analyzed for regular test systems with up to 40 first-row atoms. We also analyze the size-extensivity/size-consistency and grid- and basis set dependence of LS-THC-MP3. Overall, the errors observed are favorably small in comparison with standard density fitting, and a more streamlined method of generating grids via pruning is suggested. A practical crossover (the point at which LS-THC-MP3 is cheaper than the canonical method) is achieved around 240 correlated electrons. Despite several drawbacks of LS-THC that have been identified: an initial non-linearity of error when increasing system size, poor description of angular correlation, and a potentially large increase in error with the basis set size, the results show that LS-THC has significant potential for practical application to MP3 and other wavefunction methods.

[1]  Thomas Bondo Pedersen,et al.  Reduced scaling in electronic structure calculations using Cholesky decompositions , 2003 .

[2]  Rick A. Kendall,et al.  The impact of the resolution of the identity approximate integral method on modern ab initio algorithm development , 1997 .

[3]  Roman Schutski,et al.  Tensor-structured coupled cluster theory. , 2017, The Journal of chemical physics.

[4]  F. Weigend,et al.  Efficient use of the correlation consistent basis sets in resolution of the identity MP2 calculations , 2002 .

[5]  Joonho Lee,et al.  Systematically Improvable Tensor Hypercontraction: Interpolative Separable Density-Fitting for Molecules Applied to Exact Exchange, Second- and Third-Order Møller-Plesset Perturbation Theory. , 2019, Journal of chemical theory and computation.

[6]  F. Neese,et al.  Communication: An improved linear scaling perturbative triples correction for the domain based local pair-natural orbital based singles and doubles coupled cluster method [DLPNO-CCSD(T)]. , 2018, The Journal of chemical physics.

[7]  Zoltán Rolik,et al.  A general-order local coupled-cluster method based on the cluster-in-molecule approach. , 2011, The Journal of chemical physics.

[8]  Wei Li,et al.  Cluster-in-molecule local correlation method for post-Hartree–Fock calculations of large systems , 2016 .

[9]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[10]  W. Kutzelnigg,et al.  Comparison of CEPA and CP-MET methods , 1980 .

[11]  A. Becke.  A multicenter numerical integration scheme for polyatomic molecules , 1988 .

[12]  D. Matthews.  Improved Grid Optimization and Fitting in Least Squares Tensor Hypercontraction. , 2019, Journal of chemical theory and computation.

[13]  J. Kussmann,et al.  Linear‐Scaling Methods in Quantum Chemistry , 2007 .

[14]  Robert M Parrish,et al.  Communication: Tensor hypercontraction. III. Least-squares tensor hypercontraction for the determination of correlated wavefunctions. , 2012, The Journal of chemical physics.

[15]  D. Mazziotti,et al.  Low-rank spectral expansions of two electron excitations for the acceleration of quantum chemistry calculations. , 2012, The Journal of chemical physics.

[16]  Kohn,et al.  Density functional and density matrix method scaling linearly with the number of atoms. , 1996, Physical review letters.

[17]  Benny G. Johnson,et al.  A standard grid for density functional calculations , 1993 .

[18]  H. Kiers.  Towards a standardized notation and terminology in multiway analysis , 2000 .

[19]  Tamara G. Kolda,et al.  Tensor Decompositions and Applications , 2009, SIAM Rev..

[20]  M. Schütz,et al.  On the use of the Laplace transform in local correlation methods. , 2008, Physical chemistry chemical physics : PCCP.

[21]  F. Neese,et al.  Communication: Exact analytical derivatives for the domain-based local pair natural orbital MP2 method (DLPNO-MP2). , 2018, The Journal of chemical physics.

[22]  Yang Yang,et al.  Tensor hypercontracted ppRPA: reducing the cost of the particle-particle random phase approximation from O(r(6)) to O(r(4)). , 2014, The Journal of chemical physics.

[23]  Wolfgang Hackbusch,et al.  F ¨ Ur Mathematik in Den Naturwissenschaften Leipzig Tensor Decomposition in Post-hartree Fock Methods. I. Two-electron Integrals and Mp2 , 2022 .

[24]  Frank Neese,et al.  An efficient and near linear scaling pair natural orbital based local coupled cluster method. , 2013, The Journal of chemical physics.

[25]  A. Auer,et al.  Tensor decomposition in post-Hartree-Fock methods. II. CCD implementation. , 2013, The Journal of chemical physics.

[26]  Shuhua Li,et al.  Cluster-in-molecule local correlation method for large systems , 2013, Science China Chemistry.

[27]  Edward F. Valeev,et al.  State-Averaged Pair Natural Orbitals for Excited States: A Route toward Efficient Equation of Motion Coupled-Cluster. , 2018, Journal of chemical theory and computation.

[28]  Robert M Parrish,et al.  Tensor hypercontraction. II. Least-squares renormalization. , 2012, The Journal of chemical physics.

[29]  Hans-Joachim Werner,et al.  Scalable Electron Correlation Methods. 3. Efficient and Accurate Parallel Local Coupled Cluster with Pair Natural Orbitals (PNO-LCCSD). , 2017, Journal of chemical theory and computation.

[30]  Zoltán Rolik,et al.  An efficient linear-scaling CCSD(T) method based on local natural orbitals. , 2013, The Journal of chemical physics.

[31]  Thomas Kjærgaard,et al.  Linear-Scaling Coupled Cluster with Perturbative Triple Excitations: The Divide-Expand-Consolidate CCSD(T) Model. , 2015, Journal of chemical theory and computation.

[32]  R. Lindh,et al.  Low-cost evaluation of the exchange Fock matrix from Cholesky and density fitting representations of the electron repulsion integrals. , 2007, The Journal of chemical physics.

[33]  T. Martínez,et al.  Tensor Hypercontraction Second-Order Møller-Plesset Perturbation Theory: Grid Optimization and Reaction Energies. , 2015, Journal of chemical theory and computation.

[34]  Robert M Parrish,et al.  Tensor hypercontraction density fitting. I. Quartic scaling second- and third-order Møller-Plesset perturbation theory. , 2012, The Journal of chemical physics.

[35]  F. Neese,et al.  Exploring the Accuracy of a Low Scaling Similarity Transformed Equation of Motion Method for Vertical Excitation Energies. , 2018, Journal of chemical theory and computation.

[36]  Christof Hättig,et al.  Local explicitly correlated second- and third-order Møller-Plesset perturbation theory with pair natural orbitals. , 2011, The Journal of chemical physics.

[37]  Frank Neese,et al.  Towards a pair natural orbital coupled cluster method for excited states. , 2016, The Journal of chemical physics.

[38]  Andreas Grüneis,et al.  Low rank factorization of the Coulomb integrals for periodic coupled cluster theory. , 2016, The Journal of chemical physics.

[39]  K. Kitaura,et al.  Fragment molecular orbital method: an approximate computational method for large molecules , 1999 .

[40]  T. Martínez,et al.  Rank reduced coupled cluster theory. I. Ground state energies and wavefunctions. , 2019, The Journal of chemical physics.

[41]  Yang Yang,et al.  The tensor hypercontracted parametric reduced density matrix algorithm: coupled-cluster accuracy with O(r(4)) scaling. , 2013, The Journal of chemical physics.

[42]  N. Shenvi,et al.  Comparison of low-rank tensor expansions for the acceleration of quantum chemistry computations. , 2013, The Journal of chemical physics.

[43]  Mark S. Gordon,et al.  The Effective Fragment Potential Method: A QM-Based MM Approach to Modeling Environmental Effects in Chemistry , 2001 .

[44]  Peter M. W. Gill,et al.  SG‐0: A small standard grid for DFT quadrature on large systems , 2006, J. Comput. Chem..

[45]  Roland Lindh,et al.  Unbiased auxiliary basis sets for accurate two-electron integral approximations. , 2007, The Journal of chemical physics.

[46]  Lan Cheng,et al.  Coupled-cluster techniques for computational chemistry: The CFOUR program package. , 2020, The Journal of chemical physics.