A high
dimensional and accurate atomistic neural network potential
energy surface (ANN-PES) that describes the interaction between one
O2 molecule and a highly oriented pyrolytic graphite (HOPG)
surface has been constructed using the open-source package (aenet).
The validation of the PES is performed by paying attention to static
characteristics as well as by testing its performance in reproducing
previous ab initio molecular dynamics simulation results. Subsequently,
the ANN-PES is used to perform quasi-classical molecular dynamics
calculations of the alignment-dependent scattering of O2 from HOPG. The results are obtained for 200 meV O2 molecules
with different initial alignments impinging with a polar incidence
angle with respect to the surface normal of 22.5° on a thermalized
(110 and 300 K) graphite surface. The choice of these initial conditions
in our simulations is made to perform comparisons to recent experimental
results on this system. Our results show that the scattering of O2 from the HOPG surface is a rather direct process, that the
angular distributions are alignment dependent, and that the final
translational energy of end-on molecules is around 20% lower than
that of side-on molecules. Upon collision with the surface, the molecules
that are initially aligned perpendicular to the surface become highly
rotationally excited, whereas a very small change in the rotational
state of the scattered molecules is observed for the initial parallel
alignments. The latter confirms the energy transfer dependence on
the stereodynamics for the present system. The results of our simulations
are in overall agreement with the experimental observations regarding
the shape of the angular distributions and the alignment dependence
of the in-plane reflected molecules.