Abstract
Grazing incidence fast molecule diffraction (GIFMD) has been recently used to study a number of surfaces, but this experimental effort has not been followed, to present, by a subsequent theoretical endeavor. Aiming at filling this gap, in this work, we have carried out GIFMD simulations for the benchmark system H 2 /LiF(001). To perform our study, we have built a six-dimensional potential energy surface (6D-PES) by applying a modified version of the corrugation reducing procedure (CRP) to a set of density functional theory (DFT) energies. Based on this CRP interpolated PES, we have conducted quantum dynamics calculations using both the multiconfiguration time-dependent Hartree and the time-dependent wave packet propagation methods. We have compared the results of our GIFMD simulations with available experimental spectra. From this comparison, we have uncovered a prominent role of the interaction between the quadrupole moment of H 2 and the electric field associated with LiF(001) for specific incidence crystallographic directions. We show that, on the one hand, the molecule's initial rotation strongly affects its diffractive scattering and, on the other hand, the scattering is predominantly rotationally elastic over a wide range of incidence conditions typical for GIFMD experiments.