Modeling Many-Body Interactions in Water with Gaussian Process Regression

We report a first-principles water dimer potential that captures many-body interactions through Gaussian process regression (GPR). Modeling is upgraded from previous work by using a custom kernel function implemented through the KeOps library, allowing for much larger GPR models to be constructed and interfaced with the next-generation machine learning force field FFLUX. A new synthetic water dimer data set, called WD24, is used for model training. The resulting models can predict 90% of dimer geometries within chemical accuracy for a test set and in a simulation. The curvature of the potential energy surface is captured by the models, and a successful geometry optimization is completed with a total energy error of just 2.6 kJ mol^-1, from a starting structure where water molecules are separated by nearly 4.3 Å. Dimeric modeling of a flexible, noncrystalline system with FFLUX is shown for the first time.