Understanding the Photoinduced Desorption and Oxidation of CO on Ru(0001) Using a Neural Network Potential Energy Surface.

The study of ultrafast photoinduced dynamics of adsorbates on metal surfaces requires thorough investigation of laser-excited electrons and, in many cases, the highly excited surface lattice. While ab initio molecular dynamics with electronic friction and thermostats (Te, Tl)-AIMDEF addresses such complex modeling, it imposes severe computational costs, hindering quantitative comparison with experimental desorption probabilities. In order to bypass this limitation, we utilize the embedded atom neural network method to construct a potential energy surface (PES) for the coadsorption of CO and O on Ru(0001). Our results demonstrate that this PES not only reproduces the short-time ab initio dynamics but is also able to yield statistically significant data for long lasting trajectories that correlate well with experimental findings. Furthermore, the analysis of the laser-induced dynamics reveals the existence of a dynamic trapping state that acts as a precursor for CO desorption, and it is not observed under thermal conditions. Altogether, our results validate the underlying theoretical framework, providing robust support for the description of not only the photoinduced desorption but also the oxidation of CO in terms of nonequilibrated but thermal hot electrons and phonons.

Multicoverage Study of Femtosecond Laser-Induced Desorption of CO from Pd(111).

We study the strong coverage dependence of the femtosecond laser-induced desorption of CO from Pd(111) using molecular dynamics simulations that consistently include the effect of the laser-induced hot electrons on both the adsorbates and surface atoms. Adiabatic forces are obtained from a multicoverage neural network potential energy surface that we construct using data from density functional theory calculations for 0.33 and 0.75 monolayer (ML). Our molecular dynamics simulations performed for these two trained coverages and an additional intermediate coverage of 0.60 ML reproduce well the peculiarities of the experimental findings. The performed simulations also permit us to disentangle the relative role played by the excited electrons and phonons on the desorption process and discover interesting properties of the reaction dynamics as the relevance that the precursor physisorption well acquires during the dynamics as coverage increases.

Disentangling the role of electrons and phonons in the photoinduced CO desorption and CO oxidation on (O,CO)-Ru(0001).

The role played by electronic and phononic excitations in the femtosecond laser induced desorption and oxidation of CO coadsorbed with O on Ru(0001) is investigated using ab initio molecular dynamics with electronic friction. To this aim, simulations that account for both kind of excitations and that only consider electronic excitations are performed. Results for three different surface coverages are obtained. We unequivocally demonstrate that CO desorption is governed by phononic excitations. In the case of oxidation the low statistics does not allow to give a categorical answer. However, the analysis of the adsorbates kinetic energy gain and displacements strongly suggest that phononic excitations and surface distortion also play an important role in the oxidation process.

Why Ultrafast Photoinduced CO Desorption Dominates over Oxidation on Ru(0001).

CO oxidation on Ru(0001) is a long-standing example of a reaction that, being thermally forbidden in ultrahigh vacuum, can be activated by femtosecond laser pulses. In spite of its relevance, the precise dynamics of the photoinduced oxidation process as well as the reasons behind the dominant role of the competing CO photodesorption remain unclear. Here we use ab initio molecular dynamics with electronic friction that account for the highly excited and nonequilibrated system created by the laser to investigate both reactions. Our simulations successfully reproduce the main experimental findings: the existence of photoinduced oxidation and desorption, the large desorption to oxidation branching ratio, and the changes in the O K-edge X-ray absorption spectra attributed to the initial stage of the oxidation process. Now, we are able to monitor in detail the ultrafast CO desorption and CO oxidation occurring in the highly excited system and to disentangle what causes the unexpected inertness to the otherwise energetically favored oxidation.

Photoinduced Desorption Dynamics of CO from Pd(111): A Neural Network Approach.

Modeling the ultrafast photoinduced dynamics and reactivity of adsorbates on metals requires including the effect of the laser-excited electrons and, in many cases, also the effect of the highly excited surface lattice. Although the recent ab initio molecular dynamics with electronic friction and thermostats, (Te,Tl)-AIMDEF [Alducin, M.; Phys. Rev. Lett. 2019, 123, 246802], enables such complex modeling, its computational cost may limit its applicability. Here, we use the new embedded atom neural network (EANN) method [Zhang, Y.; J. Phys. Chem. Lett. 2019, 10, 4962] to develop an accurate and extremely complex potential energy surface (PES) that allows us a detailed and reliable description of the photoinduced desorption of CO from the Pd(111) surface with a coverage of 0.75 monolayer. Molecular dynamics simulations performed on this EANN-PES reproduce the (Te,Tl)-AIMDEF results with a remarkable level of accuracy. This demonstrates the outstanding performance of the obtained EANN-PES that is able to reproduce available density functional theory (DFT) data for an extensive range of surface temperatures (90-1000 K); a large number of degrees of freedom, those corresponding to six CO adsorbates and 24 moving surface atoms; and the varying CO coverage caused by the abundant desorption events.

High-Dimensional Atomistic Neural Network Potential to Study the Alignment-Resolved O2 Scattering from Highly Oriented Pyrolytic Graphite.

A high dimensional and accurate atomistic neural network potential energy surface (ANN-PES) that describes the interaction between one O2 molecule and a highly oriented pyrolytic graphite (HOPG) surface has been constructed using the open-source package (aenet). The validation of the PES is performed by paying attention to static characteristics as well as by testing its performance in reproducing previous ab initio molecular dynamics simulation results. Subsequently, the ANN-PES is used to perform quasi-classical molecular dynamics calculations of the alignment-dependent scattering of O2 from HOPG. The results are obtained for 200 meV O2 molecules with different initial alignments impinging with a polar incidence angle with respect to the surface normal of 22.5° on a thermalized (110 and 300 K) graphite surface. The choice of these initial conditions in our simulations is made to perform comparisons to recent experimental results on this system. Our results show that the scattering of O2 from the HOPG surface is a rather direct process, that the angular distributions are alignment dependent, and that the final translational energy of end-on molecules is around 20% lower than that of side-on molecules. Upon collision with the surface, the molecules that are initially aligned perpendicular to the surface become highly rotationally excited, whereas a very small change in the rotational state of the scattered molecules is observed for the initial parallel alignments. The latter confirms the energy transfer dependence on the stereodynamics for the present system. The results of our simulations are in overall agreement with the experimental observations regarding the shape of the angular distributions and the alignment dependence of the in-plane reflected molecules.

Ultrafast Transient Dynamics of Adsorbates on Surfaces Deciphered: The Case of CO on Cu(100).

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.

CO Stretch Vibration Lives Long on Au(111).

Measured lifetimes of the CO internal stretch mode on various metal surfaces routinely lie in the picosecond regime. These short vibrational lifetimes, which are actually reproduced by current first-principles nonadiabatic calculations, are attributed to the rapid vibrational energy loss that is caused by the facile excitation of electron-hole pairs in metals. However, this explanation was recently questioned by the huge discrepancy that exists for CO on Au(111) between the experimental vibrational lifetime that is larger than 100 ps and the previous theoretical predictions of 4.8 and 1.6 ps. Here, we show that the state-of-the-art nonadiabatic theory does reproduce the long CO lifetime measured in Au(111) provided the molecule-surface interaction is properly described. Importantly, our new results confirm that the current understanding of the adsorbates' vibrational relaxation at metal surfaces is indeed valid.

Reactive and Nonreactive Scattering of HCl from Au(111): An Ab Initio Molecular Dynamics Study.

The HCl + Au(111) system has recently become a benchmark for highly activated dissociative chemisorption, which presumably is strongly affected by electron-hole pair excitation. Previous dynamics calculations, which were based on density functional theory at the generalized gradient approximation level (GGA-DFT) for the molecule-surface interaction, have all overestimated measured reaction probabilities by at least an order of magnitude. Here, we perform ab initio molecular dynamics (AIMD) and AIMD with electronic friction (AIMDEF) calculations employing a density functional that includes the attractive van der Waals interaction. Our calculations model the simultaneous and possibly synergistic effects of surface temperature, surface atom motion, electron-hole pair excitation, the molecular beam conditions of the experiments, and the van der Waals interaction on the reactivity. We find that reaction probabilities computed with AIMDEF and the SRP32-vdW functional still overestimate the measured reaction probabilities, by a factor 18 for the highest incidence energy at which measurements were performed (≈2.5 eV). Even granting that the experiment could have underestimated the sticking probability by about a factor three, this still translates into a considerable overestimation of the reactivity by the current theory. Likewise, scaled transition probabilities for vibrational excitation from ν = 1, j = 1 to ν = 2 are overestimated by the AIMDEF theory, by factors 3-8 depending on the initial conditions modeled. Energy losses to the surface and translational energy losses are, however, in good agreement with experimental values.

Electrons and Phonons Cooperate in the Laser-Induced Desorption of CO from Pd(111).

Femtosecond laser induced desorption of CO from a CO-covered Pd(111) surface is investigated with ab initio molecular dynamics with electronic friction that incorporates effects due to the excited electronic and phononic systems, as well as out-of-phase coadsorbate interactions. Our simulations show evidence of an important electron-phonon synergy in promoting CO desorption that has largely been neglected in other similar systems. At the saturated coverage of 0.75 ML, effects due to CO-CO interadsorbate energy exchange are also important. Our dynamics simulations, in concert with site-specific desorption energy calculations, allow us to understand the large coverage dependence of the desorption yields observed in experiments.

Dynamics of N2 sticking on W(100): the decisive role of van der Waals interactions.

The reactive dynamics of N2 on W(100) has been investigated by means of quasi-classical trajectory calculations using an interpolated six-dimensional potential energy surface (PES) based on density functional theory energies obtained employing the vdW-DF2 functional. The dynamics are compared to those obtained using the PW91 functional and to experimental data. The results show that the new PES provides a significant improvement in the description of the reactivity in this system. We show that the long standing problem that constituted the large qualitative disagreement between the simulations performed with the PW91-PES and the experiments was due to the presence of energy barriers in the entrance channel that disappear when vdW forces are accounted for.

Electron-Mediated Phonon-Phonon Coupling Drives the Vibrational Relaxation of CO on Cu(100).

We bring forth a consistent theory for the electron-mediated vibrational intermode coupling that clarifies the microscopic mechanism behind the vibrational relaxation of adsorbates on metal surfaces. Our analysis points out the inability of state-of-the-art nonadiabatic theories to quantitatively reproduce the experimental linewidth of the CO internal stretch mode on Cu(100) and it emphasizes the crucial role of the electron-mediated phonon-phonon coupling in this regard. The results demonstrate a strong electron-mediated coupling between the internal stretch and low-energy CO modes, but also a significant role of surface motion. Our nonadiabatic theory is also able to explain the temperature dependence of the internal stretch phonon linewidth, thus far considered a sign of the direct anharmonic coupling.

Electronic Stopping of Slow Protons in Oxides: Scaling Properties.

Electronic stopping of slow protons in ZnO, VO_{2} (metal and semiconductor phases), HfO_{2}, and Ta_{2}O_{5} was investigated experimentally. As a comparison of the resulting stopping cross sections (SCS) to data for Al_{2}O_{3} and SiO_{2} reveals, electronic stopping of slow protons does not correlate with electronic properties of the specific material such as band gap energies. Instead, the oxygen 2p states are decisive, as corroborated by density functional theory calculations of the electronic densities of states. Hence, at low ion velocities the SCS of an oxide primarily scales with its oxygen density.

Strong Anisotropic Interaction Controls Unusual Sticking and Scattering of CO at Ru(0001).

Complete sticking at low incidence energies and broad angular scattering distributions at higher energies are often observed in molecular beam experiments on gas-surface systems which feature a deep chemisorption well and lack early reaction barriers. Although CO binds strongly on Ru(0001), scattering is characterized by rather narrow angular distributions and sticking is incomplete even at low incidence energies. We perform molecular dynamics simulations, accounting for phononic (and electronic) energy loss channels, on a potential energy surface based on first-principles electronic structure calculations that reproduce the molecular beam experiments. We demonstrate that the mentioned unusual behavior is a consequence of a very strong rotational anisotropy in the molecule-surface interaction potential. Beyond the interpretation of scattering phenomena, we also discuss implications of our results for the recently proposed role of a precursor state for the desorption and scattering of CO from ruthenium.

Vibrational Excitation of H2 Scattering from Cu(111): Effects of Surface Temperature and of Allowing Energy Exchange with the Surface.

In scattering of H2 from Cu(111), vibrational excitation has so far defied an accurate theoretical description. To expose the causes of the large discrepancies with experiment, we investigate how the feature due to vibrational excitation (the "gain peak") in the simulated time-of-flight spectrum of (v = 1, j = 3) H2 scattering from Cu(111) depends on the surface temperature (Ts) and the possibility of energy exchange with surface phonons and electron-hole pairs (ehp's). Quasi-classical dynamics calculations are performed on the basis of accurate semiempirical density functionals for the interaction with H2 + Cu(111). The methods used include the quasi-classical trajectory method within the Born-Oppenheimer static surface model, the generalized Langevin oscillator (GLO) method incorporating energy transfer to surface phonons, the GLO + friction (GLO+F) method also incorporating energy exchange with ehp's, and ab initio molecular dynamics with electronic friction (AIMDEF). Of the quasi-classical methods tested, comparison with AIMDEF suggests that the GLO+F method is accurate enough to describe vibrational excitation as measured in the experiments. The GLO+F calculations also suggest that the promoting effect of raising Ts on the measured vibrational excitation is due to an electronically nonadiabatic mechanism. However, by itself, enabling energy exchange with the surface by modeling surface phonons and ehp's leads to reduced vibrational excitation, further decreasing the agreement with experiment. The simulated gain peak is quite sensitive to energy shifts in calculated vibrational excitation probabilities and to shifts in a specific experimental parameter (the chopper opening time). While the GLO+F calculations allow important qualitative conclusions, comparison to quantum dynamics results suggests that, with the quasi-classical way of describing nuclear motion and the present box quantization method for assigning the final vibrational state, the gain peak is not yet described with quantitative accuracy. Ways in which this problem might be resolved in the future are discussed.

Electronic Stopping of Slow Protons in Transition and Rare Earth Metals: Breakdown of the Free Electron Gas Concept.

The electronic stopping cross sections (SCS) of Ta and Gd for slow protons have been investigated experimentally. The data are compared to the results for Pt and Au to learn how electronic stopping in transition and rare earth metals correlates with features of the electronic band structures. The extraordinarily high SCS observed for protons in Ta and Gd cannot be understood in terms of a free electron gas model, but are related to the high densities of both occupied and unoccupied electronic states in these metals.

Molecular dynamics simulation of O2 adsorption on Ag(110) from first principles electronic structure calculations.

We perform a detailed study of the static and dynamical properties of molecular oxygen adsorption on Ag(110) based on semi-local density functional theory (DFT) calculations and compare the results to experimental studies. For the classical dynamics calculations we use two complementary approaches, ab initio molecular dynamics and dynamics on a precalculated potential energy surface. In contrast to the molecular beam experiments, at low beam incidence energies we obtain high molecular adsorption probabilities that are related to the physisorption-like adsorption wells at the bridge sites of Ag(110). Semi-local DFT seems to overbind O2 in these wells. Based on our dynamics calculations we propose a model for adsorption in the chemisorption wells via initial adsorption in the bridge wells. In this model the measured low adsorption probabilities at low incidence energies are explained by the existence of energy barriers between the physisorption-like and chemisorption wells.

Electron-hole pair effects in methane dissociative chemisorption on Ni(111).

The dissociative chemisorption of methane on metal surfaces has attracted much attention in recent years as a prototype of gas-surface reactions in understanding the mode specific and bond selective chemistry. In this work, we systematically investigate the influence of electron-hole pair excitations on the dissociative chemisorption of CH4/CH3D/CHD3 on Ni(111). The energy dissipation induced by surface electron-hole pair excitations is modeled as a friction force introduced in the generalized Langevin equation, in which the independent atomic friction coefficients are determined within the local-density friction approximation. Quasi-classical trajectory calculations for CH4/CH3D/CHD3 have been carried out on a recently developed twelve-dimensional potential energy surface. Comparing the dissociation probabilities obtained with and without friction, our results clearly indicate that the electron-hole pair effects are generally small, both on absolute reactivity of each vibrational state and on the mode specificity and bond selectivity. Given similar observations in both water and methane dissociation processes, we conclude that electron-hole pair excitations would not play an important role as long as the reaction is direct and the interaction time between the molecule and metal electrons is relatively short.

Electronic Friction-Based Vibrational Lifetimes of Molecular Adsorbates: Beyond the Independent-Atom Approximation.

We assess the accuracy of vibrational damping rates of diatomic adsorbates on metal surfaces as calculated within the local-density friction approximation (LDFA). An atoms-in-molecules (AIM) type charge partitioning scheme accounts for intramolecular contributions and overcomes the systematic underestimation of the nonadiabatic losses obtained within the prevalent independent-atom approximation. The quantitative agreement obtained with theoretical and experimental benchmark data suggests the LDFA-AIM scheme as an efficient and reliable approach to account for electronic dissipation in ab initio molecular dynamics simulations of surface chemical reactions.

The dynamics of adsorption and dissociation of N2 in a monolayer of iron on W(110).

We study the adsorption dynamics of N2 on an expanded monolayer of Fe grown pseudomorphically on W(110). To this aim we have performed molecular dynamics simulations in a six-dimensional potential energy surface calculated within density functional theory. Our results show that N2 dissociation on this surface is a highly activated process with an energy barrier of around 1.25 eV. Regarding molecular adsorption, we find that the energetically most favorable adsorption well corresponds to a parallel orientation of the molecule with an adsorption energy of around 520 meV. However, at low molecular energies and surface temperatures, the molecules preferentially adsorb vertically to the surface with an adsorption energy of around 480 meV. A comparative analysis with the results previously obtained on a clean Fe(110) surface shows that while surface strain favors molecular adsorption of N2 in this system, it impedes dissociative adsorption. The former is consistent with the experimental observations showing that the inertness of Fe(110) towards N2 uptake is reduced in the strained surface. The latter leads us to suggest that the experimental observation of dissociated atomic N in the strained surface when increasing surface temperature must be related to the presence of step/defects at the surface.