MP2 Is Not Good Enough: Transfer Learning ML Models for Accurate VPT2 Frequencies

The calculation of the anharmonic modes of small to medium sized molecules for assigning experimentally measured frequencies to the corresponding type of molecular motions is computationally challenging at sufficiently high levels of quantum chemical theory. Here, a practical and affordable way to calculate coupled-cluster quality anharmonic frequencies using second order vibrational perturbation theory (VPT2) from machine-learned models is presented. The approach, referred to as "NN + VPT2", uses a high-dimensional neural network (PhysNet) to learn potential energy surfaces (PESs) at different levels of theory from which harmonic and VPT2 frequencies can be efficiently determined. The NN + VPT2 approach is applied to eight small to medium sized molecules (H$_2$CO, trans-HONO, HCOOH, CH$_3$OH, CH$_3$CHO, CH$_3$NO$_2$, CH$_3$COOH and CH$_3$CONH$_2$) and frequencies are reported from NN-learned models at the MP2/aug-cc-pVTZ, CCSD(T)/aug-cc-pVTZ and CCSD(T)-F12/aug-cc-pVTZ-F12 levels of theory. For the largest molecules and at the highest levels of theory, transfer learning (TL) is used to determine the necessary full-dimensional, near-equilibrium PESs. Overall, NN + VPT2 yields anharmonic frequencies to within 20 cm$^{-1}$ of experimentally determined frequencies for close to 90 % of the modes for the highest quality PES available and to within 10 cm$^{-1}$ for more than 60 % of the modes. For the MP2 PESs only around 60 % of the NN + VPT2 frequencies were within 20~cm$^{-1}$ of the experiment, with outliers up to 150 cm$^{-1}$ compared with experiment. It is also demonstrated that the approach allows to provide correct assignments for strongly interacting modes such as the OH bending and the OH torsional modes in formic acid monomer and the CO-stretch and OH-bend mode in acetic acid.

PDF Abstract