“…Such an interpolation of high-dimensional PESs is a tedious task, but in recent years, several techniques, like analytical fits 29,31,32,33,34 , tight binding representations 35,36,37,38 , genetic programming 39 and the modified Shephard method 40,41 , possibly combined with a corrugation-reduction scheme 42,43 , have been developed. In the present work we employ a very general neural network fitting scheme, which has already proven to be a powerful tool for the accurate representation of multi-dimensional PESs in similar applications 44,45,46,47 . Since the evaluation of the energy and forces from the neural network representation is about 5 to 6 orders of magnitude faster than direct ab initio calculations, a large number of MD trajectories can be calculated in the last step of the "divide and conquer" approach to obtain the sticking probabilities at various molecular kinetic energies.…”