Search papers, labs, and topics across Lattice.
This work presents high-performance parallel implementations of high-order coupled-cluster theories (CCSDT, CCSDT(Q), and CCSDTQ) within the PySCF package, addressing the computational challenges associated with their steep scaling and memory requirements. By leveraging compact triangular storage and a multithreaded tensor contraction backend, the authors achieve near-ideal thread scaling on many-core CPU architectures and extend the methods to distributed-memory systems using MPI, allowing for efficient calculations on large molecular systems. The implementations enable benchmark-quality calculations for systems with up to 100 correlated electrons, significantly broadening the applicability of high-order CC theories to chemically relevant scenarios.
Achieving benchmark-quality coupled-cluster calculations for large molecular systems is now feasible, opening new avenues for high-accuracy quantum chemistry.
High-order coupled-cluster theories with iterative triples (CCSDT), perturbative quadruples [CCSDT(Q)], and iterative quadruples (CCSDTQ) provide benchmark-quality correlation energies, but their steep computational scalings, $O(N^8), O(N^9)$, and $O(N^{10})$, together with the large memory requirements of high-order amplitude tensors, have historically limited their application to small molecules. In this work, we develop efficient open-source implementations of spin-restricted CCSDT (RCCSDT), RCCSDT(Q), RCCSDTQ, and spin-unrestricted CCSDT (UCCSDT) within the PySCF package. The shared-memory implementation combines compact triangular storage of the highest-order amplitude tensors with the multithreaded tensor contraction backend pytblis, enabling efficient use of modern many-core CPU architectures. This design delivers near-ideal thread scaling up to 90 cores and achieves wall times shorter than or comparable to existing single-node implementations for representative benchmark molecules. We further extend RCCSDT, RCCSDT(Q), and RCCSDTQ to distributed-memory architectures using MPI-based algorithms. By distributing compact high-order amplitudes across MPI ranks and overlapping communication with computation through nonblocking data transfers, the distributed implementation achieves near-ideal strong scaling on up to 32 nodes, corresponding to approximately 3,000 CPU cores. These developments substantially extend the practical reach of canonical high-order CC theory, enabling CCSDT(Q) calculations with approximately 100 correlated electrons in 450 orbitals and CCSDTQ calculations with approximately 50 correlated electrons in 115 orbitals. Applications to $\pi$-stacked noncovalent dimers, the CO dissociation energy of Cr(CO)$_6$, and the Cope rearrangement of semibullvalene demonstrate that canonical high-order CC benchmarks are now feasible for chemically realistic molecular systems.
