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.

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

Efficient non-parametric fitting of potential energy surfaces for polyatomic molecules with Gaussian processes

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

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

Neural networks vs Gaussian process regression for representing potential energy surfaces: A comparative study of fit quality and vibrational spectrum accuracy.

Permutationally invariant polynomial regression for energies and gradients, using reverse differentiation, achieves orders of magnitude speed-up with high precision compared to other machine learning methods

Permutationally Invariant Fitting of Intermolecular Potential Energy Surfaces: A Case Study of the Ne-C2H2 System.

A Gaussian Process Based Δ-Machine Learning Approach to Reactive Potential Energy Surfaces.

Construction of Reactive Potential Energy Surfaces with Gaussian Process Regression: Active Data Selection

Assessing PIP and sGDML Potential Energy Surfaces for H3O2-

Representing Globally Accurate Reactive Potential Energy Surfaces with Complex Topography by Combining Gaussian Process Regression and Neural Networks.

Neural network Gaussian processes as efficient models of potential energy surfaces for polyatomic molecules

Quantum Gaussian process model of potential energy surface for a polyatomic molecule

Comparison of permutationally invariant polynomials, neural networks, and Gaussian approximation potentials in representing water interactions through many-body expansions

Data Quality, Data Sampling and Data Fitting: A Tutorial Guide for Constructing Full-Dimensional Accurate Potential Energy Surfaces (PESs) of Molecules and Reactions

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

Gaussian Processes for Surrogate Modeling of Discharged Fuel Nuclide Compositions

Polynomial Representation of the Gaussian Process

Unified Regression Model in Fitting Potential Energy Surfaces for Quantum Dynamics

Potential Energy Surfaces from High Fidelity Fitting Ofab Initiopoints: the Permutation Invariant Polynomial - Neural Network Approach

Permutation Invariant Polynomial Neural Network Approach to Fitting Potential Energy Surfaces