We investigate the efficient combination of the canonical polyadic decomposition (CPD) and tensor hyper-contraction (THC) approaches. We first present a novel low-cost CPD solver that leverages a precomputed THC factorization of an order-4 tensor to efficiently optimize the order-4 CPD with O ( N R 2 ) scaling. With the matrix-free THC-based optimization strategy in hand, we can efficiently generate CPD factorizations of the order-4 two-electron integral tensors and develop novel electronic structure methods that take advantage of both the THC and CPD approximations. Next, we investigate the application of a combined CPD and THC approximation of the Laplace transform (LT) second-order Mo̷ller-Plesset (MP2) method. We exploit the ability to switch efficiently between the THC and CPD factorizations of the two-electron integrals to reduce the computational complexity of the LT MP2 method while preserving the accuracy of the approach. Furthermore, we take advantage of the robust fitting approximation to eliminate the leading-order error in the CPD approximated tensor networks. Finally, we show that modest values of THC and CPD rank preserve the accuracy of the LT MP2 method and that this CPD + THC LT MP2 strategy realizes a performance advantage over canonical LT MP2 in both computational wall-times and memory resource requirements.
We investigate the efficient combination of the canonical polyadic decomposition (CPD) and tensor hyper-contraction (THC) approaches. We first present a novel low-cost CPD solver which leverages a precomputed THC factorization of an order-
4 tensor to efficiently optimize the order-
4 CPD with
O(NR2) scaling. With the matrix-free THC-based optimization strategy in hand we can: efficiently generate CPD factorizations of the order-4 two-electron integral tensors; and develop novel electronic structure methods which take advantage of both the THC and CPD approximations. Next, we investigate the application of a combined CPD and THC approximation of the Laplace transform (LT) second-order Møller-Plesset (MP2) method. We exploit the ability to switch efficiently between the THC and CPD factorizations of the two electron integrals to reduce the computational complexity of the LT MP2 method while preserving the accuracy of the approach. Furthermore we take advantage of the robust fitting approximation to eliminate leading order error in the CPD approximated tensor networks. Finally, we show that modest values of THC and CPD rank preserve the accuracy of the LT MP2 method and that this CPD+THC LT MP2 strategy realizes a performance advantage over canonical LT MP2 in both computational wall-times and memory resource requirements.
Abstract An efficient parallel algorithm is developed for second‐order Møller–Plesset perturbation theory with the resolution‐of‐identity approximation of two‐electron repulsion integrals (RI‐MP2) to perform MP2 energy calculations of large molecules on distributed memory processors. Benchmark calculations are carried out for taxol (C47H51NO14), valinomycin (C54H90N6O18), and two‐layer nanographene sheets (C96H24)2, which show the high parallel efficiency of the developed algorithm. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009