Hybrid Ab Initio and Empirical Machine Learning Models for the Potential Energy Surface

We propose a methodology to generate hybrid machine learning models for the potential energy surface trained simultaneously on data from ab initio electronic structure calculations and on thermodynamic and/or structural observables from experiment. The approach is based on the use of a loss function that includes the mean square error of observables with respect to their experimental values, in addition to the usual terms involving the mean square error of the energies and forces with respect to ab initio data. We employ a reweighting procedure that allows for the calculation of ensemble averages of observables during training for arbitrary values of the model parameters and on the fly. The method is general and can be applied to any set of static observables. We illustrate the usefulness of this approach by applying it to the generation of hybrid models for water that reproduce accurately the experimental density maximum and the radial distribution function of liquid water in molecular dynamics simulations.