We implemented eq. ( ) in ASE that we use for all dynamics calculations. In the LDFA, Since the neural network potential is by orders of magnitude faster than DFT, it is easy to run many more trajectories than in ab-initio molecular dynamics. In Figure It can be seen that desorption events start as early as 0.2 ps after laser excitation, while the highest number of atoms in the desorbing channel takes place when the kinetic energy is the highest, between 0.3 and 0.4 ps after pulse excitation.