Hot carriers at metal surfaces can drive nonthermal reactions of adsorbates. Characterizing nonequilibrium statistics among various degrees of freedom in an ultrafast time scale is crucial to understand and develop hot carrier-driven chemistry. Here we demonstrate multidimensional vibrational dynamics of carbon monoxide (CO) on Cu(100) along hot-carrier induced desorption studied by using time-resolved vibrational sum-frequency generation with phase-sensitive detection. Instantaneous frequency and amplitude of the CO internal stretching mode are tracked with a subpicosecond time resolution that is shorter than the vibrational dephasing time. These experimental results in combination with numerical analysis based on Langevin simulations enable us to extract nonequilibrium distributions of external vibrational modes of desorbing molecules. Superstatistical distributions are generated with modedependent frictional couplings in a few hundred femtoseconds after hot-electron excitation, and energy flow from hot electrons and intermode anharmonic coupling play crucial roles in the subsequent evolution of the non-Boltzman distributions. DOI: 10.1103/PhysRevLett.117.186101 Hot-carrier induced surface reaction of adsorbates on metal is triggered nonthermally by high-temperature electrons, which are generated by irradiation of metal with an intense femtosecond laser pulse [1][2][3][4][5][6][7][8][9][10][11][12][13] or by excitation of the localized surface plasmon of metallic nanoparticles [14][15][16][17]. Ultrafast energy transfer from electrons in metal to adsorbate nuclear degrees of freedom (DOF) via nonadiabatic coupling [18][19][20][21][22][23][24][25] is of fundamental importance because it induces novel surface reactions [14,26,27].Despite extensive studies at flat metal surfaces over decades [26], our understanding of the ultrafast energy flow during hot-carrier induced reactions is still in its infancy. Earlier two-pulse correlation studies [1,[3][4][5][8][9][10][11]28] have specified the time scale of reactions and provided rich information on energy partitioning in desorption by product state analysis but were incapable of probing vibrational dynamics. Ultrafast time-resolved vibrational spectroscopy, particularly probing the internal stretching mode (S mode) of diatomic adsorbates [12,29,30], has been used to track excitation of external vibrational modes: frustrated translation (FT), frustrated rotation (FR), and external stretch (ES) modes. While the frequency shift of the S mode is dominated by coupling with the FT mode under moderate excitation [31][32][33][34][35], intense fs pulse excitation leads to larger S-mode frequency shifts due to hot-carrier driven FR mode excitation; thus, the FR mode was considered to trigger lateral hopping [12] and/or desorption of adsorbates [36,37]. Because the time resolution of the time-resolved vibrational sum-frequency generation (TR-VSFG) studies is inherently limited by an S-mode dephasing time (∼1 ps), it has not been possible to clarify the dynamics in the sub-p...