Towards the ‘Gold Standard’ in Bootstrap Embedding

We present a scalable implementation of Coupled Cluster Singles, Doubles, and Perturbative Triples [CCSD(T)] density matrix evaluation, integrated within the Bootstrap Embedding (BE2) framework and accelerated through the use of orbital permutation symmetry. Our approach achieves substantial reduction in computational cost compared to canonical all-electron CCSD(T) calculations, while retaining high accuracy in relative energies. Benchmark calculations on annulene systems show that BE2-CCSD(T) achieves accurate energetics with near-linear scaling in system size, enabling routine application to systems that are otherwise intractable with conventional post-Hartree–Fock methods using both double and triple-zeta basis sets. We then benchmark the BE2-CCSD(T) scheme with Intrinsic Atomic Orbital and Projected Atomic Orbital (IAO+PAO) partitioning versus canonical CCSD(T) in isomermisations. This work lays the groundwork for extending fragment-based correlated wavefunction methods to broader classes of molecular systems.