Several studies have shown that 2 H NMR, among many other techniques, is a versatile tool for the study of hydrogen coordination modes in transition metal complexes [1] or clusters.[2]Thanks to 1 H gas-phase and 2 H solid-state MAS-NMR spectroscopy, the coexistence of ancillary organic ligands, with mobile and reactive hydride ligands coordinated to ruthenium nanoparticles (NPs) has been observed.[3] Similar results were also found in the case of Ru NPs embedded in the cavities of metal-organic frameworks (MOFs).[4] However, secure assignment of the experimental information is hardly straightforward, as few reliable reference data are available for different bonding situations and chemical environments. As a matter of fact, solid-state NMR spectra of molecular compounds with nonequivalent deuterons generally consist of a complex superposition of subspectra with different shapes, which depend on the motion of the corresponding deuterons.[5] Here, the need for systematic theoretical studies for better interpretation of NMR spectra is clear.Recently, some of the authors have proposed the foundation of this work, with molecular density functional theory (DFT) calculations of quadrupole coupling constants Q cc and the asymmetry parameter h Q , (see the Supporting Information for definitions) of deuteron in ruthenium complexes [6] directly comparable with the experimental data.[1g] There is no doubt now that joint theoretical/experimental studies are a very powerful approach to complete the partial understanding of deuteron solid-state NMR spectra. Herein, we propose to explore the stability and to estimate quadrupole coupling parameters of deuterium atoms, on and below the surface of Ru NPs at low coverage, by means of DFT calculations on an infinite Ru(0001) slab surface model. Moreover, new quantitative insights into the kinetic behavior of chemically adsorbed hydrogen atoms are brought with the help of diffusion barrier estimations for surface hopping and subsurface paths. These investigations should form a first basis for the assignment of deuteron MAS-NMR spectra of ruthenium NPs.At a coverage value of 1/4 monolayer (ML), high-symmetry adsorption sites on the Ru surface, as depicted in Figure 1, were considered, yielding energetic and structural parameters in good agreement with the literature.[7] The Fcc site is slightly more stable than the Hcp site, while sites with less local symmetry order, namely On-top and Bridge, respectively lie 3 and 10 kcal mol À1 higher. Note that the terminal (On-top) and the two-fold (Bridge) sites are distinguished from the other sites by an imaginary frequency, which indicates that these states are saddle points at this particular low coverage value. The subsurface sites, Ts and Os, are significantly higher in energy, 24.2 and 8.7 kcal mol À1 respectively. Consequently, unless coverage increases significantly or co-adsorption species can reasonably invert their energy sign, there is a very low probability that hydrogen atoms can be found in these four-and six-fold sites at low co...