Related Experiment Video
Updated: Oct 13, 2026

Computation of Atmospheric Concentrations of Molecular Clusters from ab initio Thermochemistry
Published on: April 8, 2020
High-performance parallel implementation of spin-free high-order coupled-cluster theories
Yu Jin1, Christopher Hillenbrand2, Timothy C Berkelbach1,3
1Initiative for Computational Catalysis, Flatiron Institute, New York, New York 10010, USA.
Abstract:
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(N8), O(N9), and O(N10), 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), and RCCSDTQ within the pyscf package and a Message Passing Interface (MPI)-based extension, distributed-cc. 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 central processing unit (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 ∼3000 CPU cores. These developments substantially extend the practical reach of canonical high-order coupled-cluster theory, enabling CCSDT and CCSDT(Q) calculations on systems with ∼100 correlated electrons in 450 orbitals and CCSDTQ calculations on systems with ∼50 correlated electrons in 115 orbitals. Applications to π-stacked noncovalent dimers, the CO dissociation energy of Cr(CO)6, and the Cope rearrangement of semibullvalene demonstrate that canonical high-order coupled-cluster benchmarks are now feasible for chemically realistic molecular systems.
Related Concept Videos
Spin–Spin Coupling Constant: Overview
Qualitatively, any spin plus-half nucleus polarizes the spins of its electrons to the minus-half state. Consequently, the paired electron in the hydrogen–carbon bond must have a...
Spin–Spin Coupling: Two-Bond Coupling (Geminal Coupling)
The central atom need not be NMR-active because its electrons are affected by the electron polarization of the spin-active atoms. However, spin information is transmitted less effectively than in one-bond coupling, and 2J values are usually weaker than 1J values. The energy of...
Spin–Spin Coupling: One-Bond Coupling
Spin–Spin Coupling: Three-Bond Coupling (Vicinal Coupling)
The extent of coupling depends on the C‑C bond length, the two H‑C‑C angles, any electron-withdrawing substituents, and the dihedral angle between the involved orbitals. The...
¹H NMR: Interpreting Distorted and Overlapping Signals
As Δν decreases and the signals move closer, the doublets appear increasingly distorted. The intensities of the inner lines increase at the cost of those of the outer lines as the signals are slanted or...
Hybridization of Atomic Orbitals II
