A model for treating excitation and relaxation of adsorbates at metallic surfaces induced by non-adiabatic coupling is developed. The derivation is based on the concept of resonant electron transfer, where the adsorbate serves as a molecular bridge for the inelastic transition between an electron source and a sink. In this picture, energy relaxation and scanning tunneling microscopy (STM) at metallic surfaces are treated on an equal footing as a quasi-thermal process. The model goes beyond the local harmonic approximation and allows for an unbiased description of floppy systems with multiple potential wells. Further, the limitation of the product ansatz for the vibronic wave function to include the position-dependence of the non-adiabatic couplings is avoided by explicitly enforcing detailed balance. The theory is applied to the excitation of hydrogen on palladium, which has multiple local potential minima connected by low energy barriers. The main aspects investigated are the lifetimes of adsorbate vibrations in different adsorption sites, as well as the dependence of the excitation, response, and transfer rates on an applied potential bias. The excitation and relaxation simulations reveal intricate population dynamics that depart significantly from the simplistic tunneling model in a truncated harmonic potential. In particular, the population decay from an initially occupied local minimum induced by the contact with an STM tip is found to be better described by a double exponential. The two rates are interpreted as a response to the system perturbation and a transfer rate following the perturbation. The transfer rate is found to obey a power law, as was the case in previous experimental and theoretical work.
