Hot or not:An embedding technique for metallic systems makes it possible to model energy dissipation into substrate phonons during surface chemical reactions from first principles. Application to O2 dissociation at Pd(100) predicts translationally "hot" oxygen adsorbates as a consequence of the released adsorption energy (ca. 2.6 eV). This questions the instant thermalization of reaction enthalpies generally assumed in heterogeneous catalysis modeling.Modeling Heat Dissipation at the Nanoscale:An Abstract: We present an embedding technique for metallic systems that makes it possible to model energy dissipation into substrate phonons during surface chemical reactions from first principles. The separation of chemical and elastic contributions to the interaction potential provides a quantitative description of both electronic and phononic band structure. Application to the dissociation of O 2 at Pd(100) predicts translationally "hot" oxygen adsorbates as a consequence of the released adsorption energy (ca. 2.6 eV). This finding questions the instant thermalization of reaction enthalpies generally assumed in models of heterogeneous catalysis.Exothermic surface chemical reactions may easily release several electron volts of energy. Though sizable in view of potential microscopic dissipation channels, the prevalent picture in chemical kinetics is that this energy is quasi-instantaneously thermalized, ultimately into * joerg.meyer@ch.tum.de present address: Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands phononic degrees of freedom. This motivates theoretical treatments on the level of a local temperature and separates the continuous chemical motion into rare-event sequences of thermal reactions. The resulting Markovian state-to-state hopping underlies, for example, all presentday microkinetic formulations in heterogeneous catalysis. [1,2] The validation of this picture would require detailed insight into the energy-conversion process at the interface. To date this is limited at best, and if at all centered on reactions at well-defined single-crystal surfaces in ultrahigh vacuum. For a prototypical model reaction like the dissociative adsorption of O 2 molecules scanningtunneling microscopy experiments have suggested the formation of so-called "hot adatoms" on several metal surfaces. [3][4][5][6] As a consequence of the released chemical energy, this transient mobility thus intricately couples the elementary reaction steps of dissociation and diffusion.As the experimental quest to generate molecular movies of such reactions is still ongoing, theory has been challenged to elucidate the equilibration dynamics of this process. [7][8][9] Here, the bond breaking and making in highly corrugated surface potentials dictates computationally demanding quantum mechanical (QM) treatments, in par- ticular periodic boundary condition (PBC) supercell approaches to adequately describe the delocalized (surface) metallic band structure. [2] Complementing this with a...