Time-resolved vibrational spectroscopy constitutes an invaluable experimental tool for monitoring hot-carrier-induced surface reactions. However, the absence of a full understanding of the precise microscopic mechanisms causing the transient spectral changes has limited its applicability. Here we introduce a robust first-principles theoretical framework that successfully explains both the nonthermal frequency and linewidth changes of the CO internal stretch mode on Cu(100) induced by femtosecond laser pulses. Two distinct processes engender the changes: electron-hole pair excitations underlie the nonthermal frequency shifts, while electron-mediated vibrational mode coupling gives rise to linewidth changes. Furthermore, the origin and precise sequence of coupling events are finally identified.One of the ultimate goals in surface science is to comprehend the fundamental processes that bring about the specific timescales of surface reactions [1][2][3]. To acquire such a time-resolved insight, numerous experiments have studied the ultrafast elementary motions of adsorbates on metal surfaces by means of time-dependent techniques, including vibrational motion [4][5][6], molecular desorption [7][8][9][10][11], diffusion [12,13], and dissociation [14]. In these experiments, the transient condition achieved by the use of intense femtosecond laser pulses initiates the energy exchange mechanisms between the laser-excited surface electrons and the vibrational modes of the adsorbates and surface lattice.In time-resolved infrared (IR) spectroscopy experiments, the initial adsorbate dynamics commenced by the pump pulse is directly probed by tracking the frequency shift and linewidth changes of the IR-active internal stretch (IS) mode. In all adsorbate-surface systems investigated thus far [e.g., CO/Ru(001) [10], NO/Ir (111) [14], CO/Pt(111) [6,15], CO/Cu(100) [11]], the IS frequency mode exhibits an initial redshift followed by a rapid blueshift. However, the origin of such a characteristic behavior is still not fully understood. In the early works it was ascribed to anharmonic coupling with other low-energy (LE) modes [10,16,17], while later energy transfer from the laser-excited hot electrons to the adsorbate motion via nonadiabatic coupling (NC) was also considered [6,11,14,18]. Nonetheless, the lack of a general quantitative theory effectively prevents us from harnessing the full potential of time-resolved vibrational spectroscopy, as well as from extracting the information about subpicosecond dynamics of surface reactions buried within.In this Letter, we introduce a general first-principles theoretical framework that allows us to calculate directly the vibrational spectra changes and identify the spe-cific mechanisms behind them. Our theory relies on a recently developed approach based on many-body and density functional perturbation theories [19] that we extend here to treat nonequilibrium conditions. Our firstprinciples formalism accounts for electron-phonon coupling processes up to second order, including electronhole p...