Abstract:The mathematical representation of large data sets of electronic energies has seen substantial progress in the past 10 years. The so-called Permutationally Invariant Polynomial (PIP) representation is one established approach. This approach dates from 2003, when a global potential energy surface (PES) for CH5+ was reported using a basis of polynomials that are invariant with respect to the 120 permutations of the five equivalent H atoms. More recently, several approaches from "machine learning" have been applied to fit these large data sets. Gaussian Process (GP) regression is such an approach. Here, we consider the implementation of the (full) GP due to Krems and co-workers, with a modification that renders it permutationally invariant, which we denote by PIP-GP. This modification uses the approach of Guo and co-workers and later extended by Zhang and co-workers, to achieve permutational invariance for neural-network fits. The PIP, GP, and PIP-GP approaches are applied to four case studies for fitting data sets of electronic energies: H3O+, OCHCO+, and H2CO/ cis-HCOH/ trans-HCOH with the goal of assessing precision, accuracy in normal-mode analysis and barrier heights, and timings. We also report an application to (HCOOH)2, where the full PIP approach is possible but where the PIP-GP one is not feasible. However, by replicating data, which is feasible in this case, the GP approach is able to represent the data with precision comparable to that of the PIP approach. We examine these assessments for varying sizes of data sets in each case to determine the dependence of properties of the fits on the training data size. We conclude with some comments on the different aspects of computational effort of the PIP, GP, and PIP-GP approaches and also challenges these methods face for more "rugged" PESs, exemplified here by H2CO/ cis-HCOH/ trans-HCOH.

Gaussian Processes for Finite Size Extrapolation of Many-Body Simulations

Gaussian Processes for Finite Size Extrapolation of Many-Body Simulations

Residual Gaussian Process: A Tractable Nonparametric Bayesian Emulator for Multi-Fidelity Simulations

Interpolation and extrapolation of global potential energy surfaces for polyatomic systems by Gaussian processes with composite kernels

Optimal Allocation of Computational Resources Based on Gaussian Process: Application to Molecular Dynamics Simulations

Gaussian Process Regression Adaptive Density-Guided Approach: Towards Calculations of Potential Energy Surfaces for Larger Molecules

Gaussian Processes with Spectral Delta kernel for higher accurate Potential Energy surfaces for large molecules

Gaussian Processes for Surrogate Modeling of Discharged Fuel Nuclide Compositions

Solving High Frequency and Multi-Scale PDEs with Gaussian Processes

Modeling of thermophysical properties and vapor-liquid equilibrium using Gaussian process regression

Multifidelity Gaussian Processes for Predicting Shear Viscosity over Wide Ranges of Liquid State Points Based on Molecular Dynamics Simulations

Robust Gaussian Process Regression method for efficient reaction pathway optimization: application to surface processes

Efficient force field and energy emulation through partition of permutationally equivalent atoms.

Gaussian Process Emulation for High-Dimensional Coupled Systems

Interpolating many-body wave functions for accelerated molecular dynamics on the near-exact electronic surface

Accelerated Scale Bridging with Sparsely Approximated Gaussian Learning.

Accelerated Bayesian Inference for Molecular Simulations using Local Gaussian Process Surrogate Models

A Graphical Multi-Fidelity Gaussian Process Model, with Application to Emulation of Heavy-Ion Collisions

Building nonparametric $n$-body force fields using Gaussian process regression

Assessing Gaussian Process Regression and Permutationally Invariant Polynomial Approaches To Represent High-Dimensional Potential Energy Surfaces.

Space-local memory in generalized master equations: Reaching the thermodynamic limit for the cost of a small lattice simulation