Machine-learning interatomic potentials achieving CCSD(T) accuracy for systems with extended covalent networks and van der Waals interactions

Machine-learning interatomic potentials (MLIPs) enable large-scale atomistic simulations at moderate computational cost while retaining ab initio accuracy. MLIPs trained on coupled-cluster data, particularly CCSD(T), have emerged as a promising route to achieve chemical accuracy beyond the limits of density functional theory (DFT) and to incorporate non-empirical van der Waals (vdW) interactions. Most existing approaches are, however, still not straightforwardly applicable for systems with extended covalent networks such as covalent organic frameworks (COFs) due to the limited availability of CCSD(T) for periodic systems. Here we present a methodology to train MLIPs with CCSD(T) accuracy for these systems. The approach uses the Δ-learning method with a dispersion-corrected tight-binding baseline. This strategy enables training on compact molecular fragments while preserving transferability toward the periodic systems. Dispersion interactions are accounted for by adding vdW-bound multimers in the training set, and the combination with a vdW-aware tight-binding baseline allows the formally local MLIP to attain CCSD(T)-level accuracy even for systems dominated by long-range vdW forces. The resulting potential yields root-mean-square energy errors below 0.4 meV/atom on training and test sets and reproduces electronic total atomization energies, bond lengths, harmonic vibrational frequencies, and inter-molecular interaction energies for benchmark molecular systems. We apply the method to a prototypical quasi-two-dimensional COF composed of carbon and hydrogen. The COF structure, inter-layer binding energies, and hydrogen absorption are analyzed at CCSD(T) accuracy. The developed methodology opens a practical route to large-scale atomistic simulations for systems with extended covalent networks and vdW interactions with chemical accuracy.