Abstract
Multilevel coupled cluster theory offers reduced scaling computation of intensive properties in systems that are too large for standard coupled cluster calculations. A significant benefit of the multilevel coupled cluster framework is the possibility of calculating intensive properties that are not tightly localized if an appropriate set of active orbitals is used. Correlated natural transition orbitals (CNTOs) are tailored to describe excitation processes. For multilevel coupled cluster singles and doubles (MLCCSD) and singles and perturbative doubles (MLCC2) calculations, the construction of CNTOs generally becomes the computational bottleneck. Here, we demonstrate how CNTOs can be obtained with O(N3) operations, eliminating the O(N5)-scaling steps involved in the original approach. This reduction in scaling moves the bottleneck of MLCC2 and MLCCSD calculations from the active orbital space preparation to the MLCC2 and MLCCSD equations with O(N4)-scaling.