Global properties of the energy landscape: a testing and training arena for machine learned potentials

Machine learning interatomic potentials (MLIPs) have achieved remarkable accuracy on standard benchmarks, yet their ability to reproduce molecular kinetics -- critical for reaction rate calculations -- remains largely unexplored. We introduce Landscape17, a dataset of complete kinetic transition networks (KTNs) for the molecules of the MD17 dataset, computed using hybrid-level density functional theory. Each KTN contains minima, transition states, and approximate steepest-descent paths, along with energies, forces, and Hessian eigenspectra at stationary points. We develop a comprehensive test suite to evaluate the MLIP ability to reproduce these reference landscapes and apply it to a number of state-of-the-art architectures. Our results reveal limitations in current MLIPs: all the models considered miss over half of the DFT transition states and generate stable unphysical structures throughout the potential energy surface. Data augmentation with pathway configurations improves reproduction of DFT potential energy surfaces, resulting in significant improvement in the global kinetics. However, these models still produce many spurious stable structures, indicating that current MLIP architectures face underlying challenges in capturing the topology of molecular potential energy surfaces. The Landscape17 benchmark provides a straightforward but demanding test of MLIPs for kinetic applications, requiring only up to a few hours of compute time. We propose this test for validation of next-generation MLIPs targeting reaction discovery and rate prediction.