Efficient High-Level Quantum Chemical Exploration of Clathrate Hydrates via Fragmentation

We present fragmentation-based MP2 and CCSD(T)-level energetic calculations on naturally occurring clathrates encapsulating gas molecules CO ${}_2$ , CH ${}_4$ , and H ${}_2$ S, within a 20-water cage. These calculations are followed by the computation and analysis of their vibrational infrared (IR) spectra. High-accuracy single-point energy evaluations at the MP2 and CCSD(T) levels using the aug-cc-pVNZ basis sets (N = T, Q, 5), involving $\sim$ 6200 basis functions are performed on a desktop workstation with 16 cores, demonstrating the practical feasibility of such demanding computations. To the best of our knowledge, this work reports, for the first time, the complete basis set (CBS) limit at the CCSD(T) level for a system of this size. Our efficient and scalable in-house developed fragment-based algorithms, REAlgo and CIC, enable high-level correlated calculations with large basis sets, opening new avenues for the exploration of intermolecular interactions in complex molecular clusters.