“…The main challenge in constructing multi-dimensional PESs is to represent the function between the potential energies and the molecular nuclear coordinates based on the discrete ab initio data. Fitting PESs with machine learning models has been gaining popularity in recent years, and using an artificial neural network (NN) [ 39 , 40 , 41 , 42 , 43 , 44 , 45 , 46 , 47 , 48 , 49 , 50 , 51 , 52 , 53 , 54 ] or a Gaussian process (GP) [ 55 , 56 , 57 , 58 , 59 , 60 , 61 , 62 , 63 , 64 , 65 , 66 , 67 , 68 , 69 , 70 ] are the two most common approaches. GP is a kernel-based supervised statistical learning method [ 71 ], which has been widely used to solve physical chemistry problems such as mapping high-dimensional PESs and simulating quantum scattering dynamics.…”