Search papers, labs, and topics across Lattice.
This paper introduces a block tensor decomposition (BTD) and canonical polyadic decomposition (CPD) framework to reduce the computational scaling of second-order perturbation theory (PT2) to O(N^3). The method uses a dual-grid scheme for BTD and a block-based two-stage alternating least squares (ALS) for CPD to factorize the exchange channel. Applied to MP2 and renormalized PT2 (rPT2), the BTD-CPD approach achieves accuracy comparable to canonical RI-MP2 and rPT2 while significantly reducing storage requirements to O(N^2).
Forget O(N^5) scaling: this new tensor decomposition approach slashes the computational cost of second-order perturbation theory down to O(N^3) without sacrificing accuracy.
Block tensor decomposition (BTD) and canonical polyadic decomposition (CPD) are combined into a unified $O(N^3)$-scaling framework for second-order perturbation theory (PT2), demonstrated on MP2 and renormalized PT2 (rPT2). BTD constructs the tensor hyper-contraction kernel at $O(N^3)$ via a dual-grid scheme; CPD factorizes the exchange channel through a block-based two-stage ALS. An asymmetric half-kernel design applies bare Coulomb to one vertex and coupling-constant-averaged screening to the other, capturing the SOSEX component of rPT2 without a frequency-dependent CPD. For MP2, BTD-CPD reproduces canonical RI-MP2 to 0.058~kcal/mol per heavy atom. For rPT2@PBE0 on the S66x8 benchmark, the mean absolute error is 0.36~kcal/mol (ME $-$0.19, RMSE 0.46) over 528 data points. The CPD-compressed intermediates yield $O(N^2)$ storage alongside $O(N^3)$ scaling.
