Cluster-in-Molecule Method Combined with the Domain-Based Local Pair Natural Orbital Approach for Electron Correlation Calculations of Periodic Systems

The cluster-in-molecule (CIM) method was extended to systems with periodic boundary conditions (PBCs) in a previous work (PBC-CIM) [J. Chem. Theory Comput.2019, 15, 2933], which is able to compute the electronic structures of periodic systems at second-order Møller-Plesset perturbation theory (MP2) and coupled cluster singles and doubles (CCSD) levels. However, the high computational costs of CCSD with respect to the size of clusters limit the usage of PBC-CIM to crystals with small or medium unit cells. In this work, we further develop the PBC-CIM method by employing the domain-based local pair natural orbital (DLPNO) methods for the electron correlation calculations of clusters to reduce the computational costs. The combined approach allows CCSD with perturbative triples, denoted as CCSD(T), to be computationally available for accurate descriptions of periodic systems. The distant-pair correction is also implemented to improve the accuracy of PBC-CIM. As in the molecular cases, the distant pair correction significantly improves the accuracy of various PBC-CIM methods with few additional costs. The PBC-CIM-DLPNO-CCSD(T) approach has been applied to investigate the optimized lattice parameter of the cubic LiCl crystal and two adsorption problems (CO on the NaCl(100) surface and H2O on the h-BN surface). The results show that the CIM-DLPNO-CCSD(T) method offers accurate and efficient descriptions for the studied systems. Another application to the cohesive energy of the acetic acid crystal reveals that large basis sets are necessary for reliable calculations on the cohesive energies of molecular crystals.