Coupled cluster with singles, doubles and perturbative triples (CCSD(T)) often provides ground state correlation energies within "chemical accuracy," but suffers from high computational cost and steep scaling with system size. We present a GPU-accelerated implementation of CCSD(T) in the TeraChem software package. The new implementation achieves state-of-the-art performance, enabling the calculation of the (T) correction for a system with 63 atoms and more than 1000 basis functions in a little under 8 hours on a single node. Additionally, we demonstrate the utility of our optimized implementation for the rapid calculation of full CCSD(T)/CBS stacking energies for all ten unique DNA base pair stacked tetramers. We expect that the TeraChem CCSD(T) implementation will enable the rapid calculation of high-level data that was not previously accessible in a reasonable timeframe.
@article{arxiv.2512.01055,
title = {Accelerating CCSD(T) on Graphical Processing Units (GPUs)},
author = {O. Jonathan Fajen and Joseph E. Kelly and Edward G. Hohenstein and Todd J. Martínez},
journal= {arXiv preprint arXiv:2512.01055},
year = {2025}
}