The relaxation dynamics of hot H, N, and N 2 on Pd(100), Ag(111), and Fe(110), respectively, is studied by means of ab initio molecular dynamics with electronic friction. This method is adapted here to account for the electron density changes caused by lattice vibrations, thus treating on an equal footing both electron-hole (e-h) pair and phonon excitations. We find that even if the latter increasingly dominate the heavier is the hot species, the contribution of e-h pairs is by no means negligible in these cases because it gains relevance at the last stage of the relaxation process. The quantitative details of energy dissipation depend on the interplay of the potential energy surface, electronic structure, and kinetic factors. DOI: 10.1103/PhysRevB.92.201411 PACS number(s): 82.65.+r, 34.35.+a, 34.50.Bw, 68.43.−h In dynamic gas-surface environments, where gas-phase atomic and molecular species impinge on the surface at energies of the order of up to a few eV, energy dissipation occurs by the excitation of electron-hole (e-h) pairs and the excitation of lattice vibrations, i.e., phonons [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17]. In the adsorption processes of atomic and molecular species, dissociative as well as nondissociative, the species trapped by the surface gradually lose their energy until they become thermalized on the surface. The competition between the e-h pairs and phonon channels governs the relaxation dynamics of the transient hot species, and thus it plays a decisive role in the system reactivity properties. The reason is that it rules the traveled length and relaxation time of a hot atom or molecule on the surface and, consequently, the probability to undergo a recombination reaction with another adsorbate [18][19][20][21][22][23].Recent ab initio molecular dynamics (AIMD) simulations with electronic friction (AIMDEF) have shown that e-h pair excitations are the dominant relaxation mechanism for hot H atoms on Pd(100) that originate from the dissociative adsorption of H 2 [16]. More particularly, this channel dissipates energy at a five times faster rate than the phonon channel [10]. The two main reasons behind this behavior are the long H-Pd interaction time, of hundreds of fs, and the low adsorbate-tosurface atom mass ratio, γ = m H /m Pd = 0.0094. The case of H on Pd(100) represents a limiting case. For heavier adsorbates, the relative weight of e-h pairs and phonons in the energy loss is expected to vary. The energy transfer to the substrate will be determined not only by kinetic factors, such as the value of γ and the incidence conditions, but also by the topography of the multidimensional potential energy surface (PES) and the electronic structure details of the configurations probed along the relaxation trajectory.In this Rapid Communication, we investigate the relaxation dynamics of hot species in three adsorption scenarios that are representative of different energy loss regimes. We have chosen atomic N on Ag(111) (γ = 0.13), N 2 on Fe(110) (γ = 0.5), and the aforemen...