Paper detail

Implementation of the full CCSDT electronic structure model with tensor decompositions

We report a complete implementation of the coupled-cluster method with single, double, and triple excitations (CCSDT) where tensor decompositions are used to reduce its scaling and overall computational costs. For the decomposition of the electron repulsion integrals the standard density fitting (or Cholesky decomposition) format is used. The coupled-cluster single and double amplitudes are treated conventionally, and for the triple amplitudes tensor we employ the Tucker-3 compression formula, $t_{ijk}^{abc} \approx t_{XYZ} \,U^X_{ai}\,U^Y_{bj} \,U^Z_{ck}$. The auxiliary quantities $U^X_{ai}$ come from singular value decomposition (SVD) of an approximate triple amplitudes tensor based on perturbation theory. The efficiency of the proposed method relies on an observation that the dimension of the "compressed" tensor $t_{XYZ}$ sufficient to deliver a constant relative accuracy of the correlation energy grows only linearly with the size of the system, $N$. This fact, combined with proper factorization of the coupled-cluster equations, leads to practically $N^6$ scaling of the computational costs of the proposed method, as illustrated numerically for linear alkanes with increasing chain length. This constitutes a considerable improvement over the $N^8$ scaling of the conventional (uncompressed) CCSDT theory. The accuracy of the proposed method is verified by benchmark calculations of total and relative energies for several small molecular systems and comparison with the exact CCSDT method. The accuracy levels of 1 kJ/mol are easily achievable with reasonable SVD subspace size, and even more demanding levels of accuracy can be reached with a considerable reduction of the computational costs. Extensions of the proposed method to include higher excitations are briefly discussed, along with possible strategies of reducing other residual errors.