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.

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.

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.

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.

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.

Dissociative dynamics of O2 on Ag(110).

We study the dissociative dynamics of O2 on Ag(110) by performing classical and quasiclassical trajectory calculations on an adiabatic six-dimensional potential energy surface (PES). The PES is constructed from the interpolation of a large set of energies that are calculated using spin-polarized density functional theory. The minimum energy barrier to dissociation amounts to 0.36 eV. This value, which is considerably lower than the barriers of about 1.1 eV found in the Ag(100) and Ag(111) surfaces, is in line with the measured much higher reactivity of the (110) surface. Our classical dynamics calculations show that under normal incidence conditions no significant dissociation occurs below an initial energy of 0.9 eV (0.6 eV in the quasiclassical calculations). This result is an indication of a very much reduced configurational space leading to dissociation and also explains why direct dissociation has not been observed experimentally at low incidence energies. Our calculations also show that for off-normal incidence, most of the dissociation takes place close to the long-bridge site, a region of the configurational space where the energy barriers to dissociation are higher than 0.7 eV, resulting in still lower dissociation probabilities.

Influence of the van der Waals interaction in the dissociation dynamics of N2 on W(110) from first principles.

Using ab initio molecular dynamics (AIMD) calculations, we investigate the role of the van der Waals (vdW) interaction in the dissociative adsorption of N2 on W(110). Hitherto, existing classical dynamics calculations performed on six-dimensional potential energy surfaces based on density functional theory (DFT), and the semi-local PW91 and RPBE [Hammer et al. Phys. Rev. B 59, 7413 (1999)] exchange-correlation functionals were unable to fully describe the dependence of the initial sticking coefficient on the molecular beam incidence conditions as found in experiments. N2 dissociation on W(110) was shown to be very sensitive not only to short molecule-surface distances but also to large distances where the vdW interaction, not included in semilocal-DFT, should dominate. In this work, we perform a systematic study on the dissociative adsorption using a selection of existing non-local functionals that include the vdW interaction (vdW-functionals). Clearly, the inclusion of the non-local correlation term contributes in all cases to correct the unrealistic energy barriers that were identified in the RPBE at large molecule-surface distances. Among the tested vdW-functionals, the original vdW-DF by Dion et al. [Phys. Rev. Lett. 92, 246401 (2004)] and the ulterior vdW-DF2 give also an adequate description of the N2 adsorption energy and energy barrier at the transition state, i.e., of the properties that are commonly used to verify the quality of any exchange-correlation functional. However, the results of our AIMD calculations, which are performed at different incidence conditions and hence extensively probe the multi-configurational potential energy surface of the system, do not seem as satisfactory as the preliminary static analysis suggested. When comparing the obtained dissociation probabilities with existing experimental data, none of the used vdW-functionals seems to provide altogether an adequate description of the N2/W(110) interaction at short and large distances.

Vibrational lifetimes of hydrogen on lead films: an ab initio molecular dynamics with electronic friction (AIMDEF) study.

Using density functional theory and Ab Initio Molecular Dynamics with Electronic Friction (AIMDEF), we study the adsorption and dissipative vibrational dynamics of hydrogen atoms chemisorbed on free-standing lead films of increasing thickness. Lead films are known for their oscillatory behaviour of certain properties with increasing thickness, e.g., energy and electron spillout change in discontinuous manner, due to quantum size effects [G. Materzanini, P. Saalfrank, and P. J. D. Lindan, Phys. Rev. B 63, 235405 (2001)]. Here, we demonstrate that oscillatory features arise also for hydrogen when chemisorbed on lead films. Besides stationary properties of the adsorbate, we concentrate on finite vibrational lifetimes of H-surface vibrations. As shown by AIMDEF, the damping via vibration-electron hole pair coupling dominates clearly over the vibration-phonon channel, in particular for high-frequency modes. Vibrational relaxation times are a characteristic function of layer thickness due to the oscillating behaviour of the embedding surface electronic density. Implications derived from AIMDEF for frictional many-atom dynamics, and physisorbed species will also be given.

Surface strain improves molecular adsorption but hampers dissociation for N2 on the Fe/W(110) surface.

We compare the adsorption dynamics of N(2) on the unstrained Fe(110) and on a 10% expanded Fe monolayer grown on W(110) by performing classical molecular dynamics simulations that use potential energy surfaces calculated with density functional theory. Our results allow us to understand why, experimentally, the molecular adsorption of N(2) is observed on the strained layer but not on Fe(110). Surprisingly, we also find that while surface strain favors the molecular adsorption of N(2) it seems, on the contrary, to impede the dissociative adsorption. This result contrasts with previous examples for which strain is found to modify equally the energetics of chemisorption and dissociation.

Role of physisorption states in molecular scattering: a semilocal density-functional theory study on O2/Ag(111).

We simulate the scattering of O2 from Ag(111) with classical dynamics simulations performed on a six-dimensional potential energy surface calculated within semilocal density-functional theory. The enigmatic experimental trends that originally required the conjecture of two types of repulsive walls, arising from a physisorption and chemisorption part of the interaction potential, are fully reproduced. Given the inadequate description of the physisorption properties in semilocal density-functional theory, our work casts severe doubts on the prevalent notion to use molecular scattering data as indirect evidence for the existence of such states.

Electronic friction dominates hydrogen hot-atom relaxation on Pd(100).

We study the dynamics of transient hot H atoms on Pd(100) that originated from dissociative adsorption of H2. The methodology developed here, denoted AIMDEF, consists of ab initio molecular dynamics simulations that include a friction force to account for the energy transfer to the electronic system. We find that the excitation of electron-hole pairs is the main channel for energy dissipation, which happens at a rate that is five times faster than energy transfer into Pd lattice motion. Our results show that electronic excitations may constitute the dominant dissipation channel in the relaxation of hot atoms on surfaces.

Vibrational deexcitation and rotational excitation of H2 and D2 scattered from Cu(111): adiabatic versus non-adiabatic dynamics.

We have studied survival and rotational excitation probabilities of H(2)(v(i) = 1, J(i) = 1) and D(2)(v(i) = 1, J(i) = 2) upon scattering from Cu(111) using six-dimensional (6D) adiabatic (quantum and quasi-classical) and non-adiabatic (quasi-classical) dynamics. Non-adiabatic dynamics, based on a friction model, has been used to analyze the role of electron-hole pair excitations. Comparison between adiabatic and non-adiabatic calculations reveals a smaller influence of non-adiabatic effects on the energy dependence of the vibrational deexcitation mechanism than previously suggested by low-dimensional dynamics calculations. Specifically, we show that 6D adiabatic dynamics can account for the increase of vibrational deexcitation as a function of the incidence energy, as well as for the isotope effect observed experimentally in the energy dependence for H(2)(D(2))/Cu(100). Furthermore, a detailed analysis, based on classical trajectories, reveals that in trajectories leading to vibrational deexcitation, the minimum classical turning point is close to the top site, reflecting the multidimensionally of this mechanism. On this site, the reaction path curvature favors vibrational inelastic scattering. Finally, we show that the probability for a molecule to get close to the top site is higher for H(2) than for D(2), which explains the isotope effect found experimentally.

Diffusion of hydrogen in Pd assisted by inelastic ballistic hot electrons.

Sykes et al. [Proc. Natl. Acad. Sci. U.S.A. 102, 17907 (2005)] have reported how electrons injected from a scanning tunneling microscope modify the diffusion rates of H buried beneath Pd(111). A key point in that experiment is the symmetry between positive and negative voltages for H extraction, which is difficult to explain in view of the large asymmetry in Pd between the electron and hole densities of states. Combining concepts from the theory of ballistic electron microscopy and electron-phonon scattering we show that H diffusion is driven by the s-band electrons only, which explains the observed symmetry.

Dissociative and non-dissociative adsorption dynamics of N2 on Fe(110).

We study the adsorption dynamics of N(2) on the Fe(110) surface. Classical molecular dynamics calculations are performed on top of a six-dimensional potential energy surface calculated within density functional theory. Our results show that N(2) dissociation on this surface is a highly activated process that takes place along a very narrow reaction path with an energy barrier of around 1.1 eV, which explains the measured low reactivity of this system. By incorporating energy exchange with the lattice in the dynamics, we also study the non-dissociative molecular adsorption process. From the analysis of the potential energy surface, we observe the presence of two distinct N(2) adsorption wells. Our dynamics calculations show that the relative population of these adsorption sites varies with the incident energy of the molecule and the surface temperature. We find an activation energy of around 150 meV that prevents molecular adsorption under thermal and hypothermal N(2) gas exposure of the surface. This finding is also consistent with the available experimental information.

Competition between electron and phonon excitations in the scattering of nitrogen atoms and molecules off tungsten and silver metal surfaces.

We investigate the role played by electron-hole pair and phonon excitations in the interaction of reactive gas molecules and atoms with metal surfaces. We present a theoretical framework that allows us to evaluate within a full-dimensional dynamics the combined contribution of both excitation mechanisms while the gas particle-surface interaction is described by an ab initio potential energy surface. The model is applied to study energy dissipation in the scattering of N(2) on W(110) and N on Ag(111). Our results show that phonon excitation is the dominant energy loss channel, whereas electron-hole pair excitations represent a minor contribution. We substantiate that, even when the energy dissipated is quantitatively significant, important aspects of the scattering dynamics are well captured by the adiabatic approximation.

Non-reactive scattering of N2 from the W(110) surface studied with different exchange-correlation functionals.

The non-reactive scattering of N(2) from the W(110) surface is studied with six dimensional (6D) classical dynamics and two distinct potential energy surfaces (PES). Here, we use the PESs calculated with density functional theory and two different exchange-correlation functionals, the PW91 [J. E. Perdew et al., Phys. Rev. B, 1992, 46, 6671] and the RPBE [B. Hammer et al., Phys. Rev. B, 1999, 59, 7413]. By analyzing the final rotational state and angular distributions, we extract information on the characteristics of the two PESs in the 6D configurational space. Comparison of the theoretical results with the available experimental data provides detailed information on the validity of each functional. In general, the PW91 PES is more corrugated than the RPBE one in all the configurational space, meaning that there is a stronger dependence of the potential energy on the molecular orientation and position over the surface unit cell. Furthermore, we find that the larger corrugation and the less repulsive character exhibited by the PW91 PES seems to be realistic at distances above the chemisorption well. In contrast, the less corrugated RPBE PES performs better in the region below the chemisorption well.