We propose a reformulation of the traditional (T) triples correction to the coupled cluster singles and doubles (CCSD) energy in terms of local Hartree-Fock (HF) orbitals such that its structural form aligns with our recently developed linear-scaling divide-expand-consolidate (DEC) coupled cluster family of local correlation methods. In a DEC-CCSD(T) calculation, a basis of local occupied and virtual HF orbitals is used to partition the correlated calculation on the full system into a number of independent atomic fragment and pair fragment calculations, each performed within a truncated set of the complete orbital space. In return, this leads to a massively parallel algorithm for the evaluation of the DEC-CCSD(T) correlation energy, which formally scales linearly with the size of the full system and has a tunable precision with respect to a conventional CCSD(T) calculation via a single energy-based input threshold. The theoretical developments are supported by proof of concept DEC-CCSD(T) calculations on a series of medium-sized molecular systems.

The Resolution of the Identity second-order Møller-Plesset perturbation theory (RI-MP2) method is implemented within the linear-scaling Divide-Expand-Consolidate (DEC) framework. In a DEC calculation, the full molecular correlated calculation is replaced by a set of independent fragment calculations each using a subset of the total orbital space. The number of independent fragment calculations scales linearly with the system size, rendering the method linear-scaling and massively parallel. The DEC-RI-MP2 method can be viewed as an approximation to the DEC-MP2 method where the RI approximation is utilized in each fragment calculation. The individual fragment calculations scale with the fifth power of the fragment size for both methods. However, the DEC-RI-MP2 method has a reduced prefactor compared to DEC-MP2 and is well-suited for implementation on massively parallel supercomputers, as demonstrated by test calculations on a set of medium-sized molecules. The DEC error control ensures that the standard RI-MP2 energy can be obtained to the predefined precision. The errors associated with the RI and DEC approximations are compared, and it is shown that the DEC-RI-MP2 method can be applied to systems far beyond the ones that can be treated with a conventional RI-MP2 implementation.

We present a local framework for the calculation of coupled cluster excitation energies of large molecules (LoFEx). The method utilizes time-dependent Hartree-Fock information about the transitions of interest through the concept of natural transition orbitals (NTOs). The NTOs are used in combination with localized occupied and virtual Hartree-Fock orbitals to generate a reduced excitation orbital space (XOS) specific to each transition where a standard coupled cluster calculation is carried out. Each XOS is optimized to ensure that the excitation energies are determined to a predefined precision. We apply LoFEx in combination with the RI-CC2 model to calculate the lowest excitation energies of a set of medium-sized organic molecules. The results demonstrate the black-box nature of the LoFEx approach and show that significant computational savings can be gained without affecting the accuracy of CC2 excitation energies.

The Divide‐Expand‐Consolidate (DEC) scheme is a linear‐scaling and massively parallel framework for high accuracy coupled cluster (CC) calculations on large molecular systems. It is designed as a black‐box method, which ensures error control in the correlation energy and molecular properties. DEC is combined with a massively parallel implementation to fully utilize modern manycore architectures providing a fast time to solution. The implementation ensures performance portability and will straightforwardly benefit from new hardware developments. The DEC scheme has been applied to several levels of CC theory and extended the range of application of those methods. WIREs Comput Mol Sci 2017, 7:e1319. doi: 10.1002/wcms.1319 This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Software > Quantum Chemistry

The theoretical foundation for solving coupled cluster singles and doubles (CCSD) amplitude equations to a desired precision in terms of independent fragment calculations using restricted local orbital spaces is reinvestigated with focus on the individual error sources. Four different error sources are identified theoretically and numerically and it is demonstrated that, for practical purposes, local orbital spaces for CCSD calculations can be identified from calculations at the MP2 level. The development establishes a solid theoretical foundation for local CCSD calculations for the independent fragments, and thus for divide-expand-consolidate coupled cluster calculations for large molecular systems with rigorous error control. Based on this theoretical foundation, we have developed an algorithm for determining the orbital spaces needed for obtaining the single fragment energies to a requested precision and numerically demonstrated the robustness and precision of this algorithm.

We present a new framework for calculating coupled cluster (CC) excitation energies at a reduced computational cost. It relies on correlated natural transition orbitals (NTOs), denoted CIS(D')-NTOs, which are obtained by diagonalizing generalized hole and particle density matrices determined from configuration interaction singles (CIS) information and additional terms that represent correlation effects. A transition-specific reduced orbital space is determined based on the eigenvalues of the CIS(D')-NTOs, and a standard CC excitation energy calculation is then performed in that reduced orbital space. The new method is denoted CorNFLEx (Correlated Natural transition orbital Framework for Low-scaling Excitation energy calculations). We calculate second-order approximate CC singles and doubles (CC2) excitation energies for a test set of organic molecules and demonstrate that CorNFLEx yields excitation energies of CC2 quality at a significantly reduced computational cost, even for relatively small systems and delocalized electronic transitions. In order to illustrate the potential of the method for large molecules, we also apply CorNFLEx to calculate CC2 excitation energies for a series of solvated formamide clusters (up to 4836 basis functions).

The cluster perturbation series, CPS(D), for coupled cluster singles and doubles excitation energies is considered. It is demonstrated that the second-order model CPS(D-2) is identical to the configuration interaction singles with perturbative doubles, CIS(D) model. The thirdorder model, CPS(D-3), provides excitation energies of coupled cluster singles and doubles (CCSD) quality in the sense that the difference between CPS(D-3) and CCSD excitation energies is of the same size or smaller than the effect of adding triples corrections to CCSD excitation energies. We further show that the third-order corrections can be efficiently implemented, in particular, when the resolution of the identity approximation is used for integrals. We also show that the CPS(D-3) excitation energies can be determined for system sizes that are far beyond what can be considered in conventional CCSD excitation energy calculations.

We present a scalable cross-platform hybrid MPI / OpenMP / OpenACC implementation of the Divide-Expand-Consolidate (DEC) formalism with portable performance on heterogeneous HPC architectures. The Divide-Expand-Consolidate formalism is designed to reduce the steep computational scaling of conventional many-body methods employed in electronic structure theory to linear scaling, while providing a simple mechanism for controlling the error introduced by this approximation. Our massively parallel implementation of this general scheme has three levels of parallelism, being a hybrid of the loosely coupled task-based parallelization approach and the conventional MPI+X programming model, where X is either OpenMP or OpenACC. We demonstrate strong and weak scalability of this implementation on heterogeneous HPC systems, namely on the GPU-based Cray XK7 Titan supercomputer at the Oak Ridge National Laboratory. Using the "resolution of the identity second-order Møller-Plesset perturbation theory" (RI-MP2) as the physical model for simulating correlated electron motion, the linear-scaling DEC implementation is applied to 1-aza-adamantane-trione (AAT) supramolecular wires containing up to 40 monomers (2440 atoms, 6800 correlated electrons, 24440 basis functions and 91280 auxiliary functions). This represents the largest molecular system treated at the MP2 level of theory, demonstrating an efficient removal of the scaling wall pertinent to conventional quantum many-body methods.