The inelastic scattering of D2 and HD molecules impinging on a graphite surface in well-defined initial roto-vibrational states has been studied by using the computational setup recently developed to characterize important selectivities in the molecular dynamics occurring at the gas-surface interface. In order to make an immediate comparison of determined elastic and inelastic scattering probabilities, we considered for D2 and HD molecules the same initial states, as well as the same collision energy range, previously selected for the investigation of H2 behaviour. The analysis of the back-scattered molecules shows that, while low-lying initial vibrational states are preserved, the medium–high initial ones give rise to final states covering the complete ladder of vibrational levels, although with different probability for the various cases investigated. Moreover, propensities in the formation of the final rotational states are found to depend strongly on the initial ones, on the collision energy, and on the isotopologue species.