HF trimer: A new full-dimensional potential energy surface and rigorous 12D quantum calculations of vibrational states

Abstract

HF trimer, as the lightest cyclic hydrogen-bonded (HB) trimer, has long been a favorite prototype system for spectroscopic and theoretical investigations of the structure, energetics, and dynamics of hydrogen-bond networks, and the role that nonadditive, three-body interactions play in shaping these properties. Recently, rigorous quantum 12D calculations of the coupled intra- and intermolecular vibrations of this fundamental HB trimer [P.M. Felker and Z. Bačić, J. Chem. Phys. 2023, 158, 234109] were performed, employing an older ab initio-based many-body potential energy surface (PES). While the theoretical results were found to be in reasonably good agreement with the available spectroscopic data, it was also evident that it is highly desirable to develop a more accurate 12D PES of HF trimer. Motivated by this, here we report a new 12D PES of this paradigmatic system. Approximately 42,540 geometries were sampled and calculated at the level of CCSD(T)-F12a/AVTZ. The permutationally invariant polynomial-neural network based Δ-machine learning approach [Y. Liu and J. Li, J. Phys. Chem. Lett. 2022, 13, 4729] was employed to perform cost-efficient calculations of the basis-set-superposition error (BSSE) correction. By strategically selecting data points, this approach facilitated the construction of a high-precision PES with BSSE correction, while requiring only a minimal number of BSSE value computations. The fitting error of the final PES is only 0.035 kcal/mol. To assess its performance, the 12D fully coupled quantum calculations of excited intra- and intermolecular vibrational states of HF trimer are carried out using the rigorous methodology developed by us earlier. The results are found to be in a significantly better agreement with the available spectroscopic data than those obtained with the previously existing 12D PES.