Quasiclassical molecular dynamics simulations are performed to investigate energy dissipation to the (100) and (110) tungsten surfaces upon Eley−Rideal (ER) recombination of H 2 and N 2 . Calculations are carried out within the single adsorbate limit under normal incidence. A generalized Langevin surface oscillator (GLO) scheme is used to simulate the coupling to phonons, whereas electron−hole (e-h) pair excitations are implemented using the local density friction approximation (LDFA). Phonon excitations are found to reduce the ER reactivity for N 2 recombination, but do not affect H abstraction. In contrast, the effect of e-h pair excitations on the ER recombination cross section is small for N 2 , but can be important for H 2 . The analysis of the energy lost by the recombined species shows that most of the energy is dissipated into phonon excitations in the N 2 recombination and into electronic excitations in the H 2 recombination. In all cases, the energy dissipated into e-h pairs is taken away from the translational kinetic energy of the formed molecules, whereas dissipation to phonons, only significant for N 2 , also affects vibration. Interestingly, the electron mediated energy losses are found to be smaller in the case of N 2 when surface motion is allowed.