Search papers, labs, and topics across Lattice.
The paper introduces a canonical polyadic decomposition (CPD)-based low-level solver to accelerate the environment-level solver within the MPCC embedding framework. By factorizing density-fitting two-electron integral (DF TEI) tensors using CPD, the method reduces storage complexity from $\mathcal{O}(N^3)$ to $\mathcal{O}(NR)$ and computational scaling from $\mathcal{O}(N^4)$ to $\mathcal{O}(NR^2)$. Benchmarks on water clusters and alkane chains demonstrate that CPD-compressed DF TEI tensors reproduce the DF reference convergence behavior with only small, rank-controlled shifts in absolute energies, preserving chemically relevant energy differences.
Squeezing quantum chemistry calculations: CPD tensor decomposition slashes storage and compute costs in MPCC embedding, opening the door to simulating larger, more complex chemical environments.
We incorporate a canonical polyadic decomposition (CPD) based low-level solver as a means to accelerate the environment-level solver for the recently developed MPCC embedding framework. Using CPD, we both factorize the three dominant order-three density-fitting two-electron integral (DF TEI) tensors and develop a novel formulation that reduces the storage complexity of the low-level solver from $\mathcal{O}(N^3)$ to $\mathcal{O}(NR)$, where $R$ is the CPD rank, and the computational scaling of the most time-consuming contractions from $\mathcal{O}(N^4)$ to $\mathcal{O}(NR^2)$. We provide benchmarks on representative chemical environments, namely water clusters $\mathrm{(H_2O)_n}$ with $n = 1$ to $6$ and linear alkane chains $\mathrm{C_nH_{2n+2}}$ with $n = 1$ to $6$. For both test sets, using the CPD-compressed DF TEI tensors reproduces the DF reference convergence behavior of the low-level solver, the subsequent high-level step, and the fully self-consistent MPCC iterations, while introducing only small, rank-controlled shifts in absolute energies. At a fixed tolerance in the absolute MPCC energy, the CP ranks required for these tensor approximations increase linearly with system size. Chemically relevant energy differences are likewise preserved, as demonstrated for water-cluster dissociation energies and in a proof-of-concept embedding calculation of methane in a four-water cluster.
