How adiabatic is activated adsorption/associative desorption?

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.

The effects of electron-hole pair coupling on the infrared laser-controlled vibrational excitation of NO on Au(111)

In this work, we present theoretical simulations of laser-driven vibrational control of NO adsorbed on a gold surface. Our goal is to tailor laser pulses to selectively excite specific modes and vibrational eigenstates, as well as to favor photodesorption of the adsorbed molecule. To this end, various control schemes and algorithms are applied. For adsorbates at metallic surfaces, the creation of electron-hole pairs in the substrate is known to play a dominant role in the transfer of energy from the system to the surroundings. These nonadiabatic couplings are included perturbatively in our reduced density matrix simulations using a generalization of the state-resolved position-dependent anharmonic rate model we recently introduced. An extension of the reduced density matrix is also proposed to provide a sound model for photodesorption in dissipative systems.

Electronic excitations induced by hydrogen surface chemical reactions on gold

Associated with chemical reactions at surfaces energy may be dissipated exciting surface electronic degrees of freedom. These excitations are detected using metal-insulator-metal (MIM) heterostructures (Ta-TaOx-Au) and the reactions of H with and on a Au surface are probed. A current corresponding to 5×10(-5) electrons per adsorbing H atom and a marked isotope effect are observed under steady-state conditions. Analysis of the current trace when the H atom flux is intermitted suggests that predominantly the recombination reaction creates electronic excitations. Biasing the front versus the back electrode of the MIM structure provides insights into the spectrum of electronic excitations. The observed spectra differ for the two isotopes H and D and are asymmetric when comparing negative and positive bias voltages. Modeling indicates that the excited electrons and the concurrently created holes differ in their energy distributions.

Apparent failure of the Born-Oppenheimer static surface model for vibrational excitation of molecular hydrogen on copper

The accuracy of dynamical models for reactive scattering of molecular hydrogen, H(2), from metal surfaces is relevant to the validation of first principles electronic structure methods for molecules interacting with metal surfaces. The ability to validate such methods is important to progress in modeling heterogeneous catalysis. Here, we study vibrational excitation of H(2) on Cu(111) using the Born-Oppenheimer static surface model. The potential energy surface (PES) used was validated previously by calculations that reproduced experimental data on reaction and rotationally inelastic scattering in this system with chemical accuracy to within errors ≤ 1 kcal/mol ≈ 4.2 kJ/mol [Díaz C, et al. (2009) Science 326:832-834]. Using the same PES and model, our dynamics calculations underestimate the contribution of vibrational excitation to previously measured time-of-flight spectra of H(2) scattered from Cu(111) by a factor 3. Given the accuracy of the PES for the experiments for which the Born-Oppenheimer static surface model is expected to hold, we argue that modeling the effect of the surface degrees of freedom will be necessary to describe vibrational excitation with similar high accuracy.

Quantum corrected Langevin dynamics for adsorbates on metal surfaces interacting with hot electrons

We investigate the importance of including quantized initial conditions in Langevin dynamics for adsorbates interacting with a thermal reservoir of electrons. For quadratic potentials the time evolution is exactly described by a classical Langevin equation and it is shown how to rigorously obtain quantum mechanical probabilities from the classical phase space distributions resulting from the dynamics. At short time scales, classical and quasiclassical initial conditions lead to wrong results and only correctly quantized initial conditions give a close agreement with an inherently quantum mechanical master equation approach. With CO on Cu(100) as an example, we demonstrate the effect for a system with ab initio frictional tensor and potential energy surfaces and show that quantizing the initial conditions can have a large impact on both the desorption probability and the distribution of molecular vibrational states.

Role of electron-hole pair excitations in the dissociative adsorption of diatomic molecules on metal surfaces

We quantitatively evaluate the contribution of electron-hole pair excitations to the reactive dynamics of H2 on Cu(110) and N2 on W(110), including the six dimensionality of the process in the entire calculation. The interaction energy between molecule and surface is represented by an ab initio six-dimensional potential energy surface. Electron friction coefficients are calculated with density functional theory in a local density approximation. Contrary to previous claims, only minor differences between the adiabatic and nonadiabatic results for dissociative adsorption are found. Our calculations demonstrate the validity of the adiabatic approximation to analyze adsorption dynamics in these two representative systems.

Dissociative chemisorption of H2 on the Cu(110) surface: a quantum and quasiclassical dynamical study

Six-dimensional quantum dynamical and quasiclassical trajectory (QCT) calculations are reported for the reaction and vibrationally inelastic scattering of (v = 0,1,j = 0) H(2) scattering from Cu(110), and for the reaction and rovibrationally elastic and inelastic scattering of (v = 1,j = 1) H(2) scattering from Cu(110). The dynamics results were obtained using a potential energy surface obtained with density functional theory using the PW91 functional. The reaction probabilities computed with quantum dynamics for (v = 0,1,j = 0) were in excellent agreement with the QCT results obtained earlier for these states, thereby validating the QCT approach to sticking of hydrogen on Cu(110). The vibrational de-excitation probability P(v=1,j = 0 --> v = 0) computed with the QCT method is in remarkably good agreement with the quantum dynamical results for normal incidence energies E(n) between 0.2 and 0.6 eV. The QCT result for the vibrational excitation probability P(v = 0,j = 0 --> v = 1) is likewise accurate for E(n) between 0.8 and 1 eV, but the QCT method overestimates vibrational excitation for lower E(n). The QCT method gives probabilities for rovibrationally (in)elastic scattering, P(v = 1,j = 1 --> v('),j(')), which are in remarkably good agreement with quantum dynamical results. The rotationally averaged, initial vibrational state-selective reaction probability obtained with QCT agrees well with the initial vibrational state-selective reaction probability extracted from molecular beam experiments for v = 1, for the range of collision energies for which the v=1 contribution to the measured total sticking probability dominates. The quantum dynamical probabilities for rovibrationally elastic scattering of (v = 1,j = 1) H(2) from Cu(110) are in good agreement with experiment for E(n) between 0.08 and 0.25 eV.

Reactive and nonreactive scattering of N2 from Ru(0001): a six-dimensional adiabatic study

We have studied the dissociative chemisorption and scattering of N(2) on and from Ru(0001), using a six-dimensional quasiclassical trajectory method. The potential energy surface, which depends on all the molecular degrees of freedom, has been built applying a modified Shepard interpolation method to a data set of results from density functional theory, employing the RPBE generalized gradient approximation. The frozen surface and Born-Oppenheimer [Ann. Phys. (Leipzig) 84, 457 (1927)] approximations were used, neglecting phonons and electron-hole pair excitations. Dissociative chemisorption probabilities are found to be very small even for translational energies much higher than the minimum reaction barrier, in good agreement with experiment. A comparison to previous low dimensional calculations shows the importance of taking into account the multidimensional effects of N(2) rotation and translation parallel to the surface. The new calculations strongly suggest a much smaller role of nonadiabatic effects than previously assumed on the basis of a comparison between low dimensional results and experiments [J. Chem. Phys. 115, 9028 (2001)]. Also in agreement with experiment, our theoretical results show a strong dependence of reaction on the initial vibrational state. Computed angular scattering distributions and parallel translation energy distributions are in good agreement with experiments on scattering, but the theory overestimates vibrational and rotational excitations in scattering.

Spectrum of electronic excitations due to the adsorption of atoms on metal surfaces

The time-dependent, mean-field Newns-Anderson model for a spin-polarized adsorbate approaching a metallic surface is solved in the wide-band limit. Equations for the time evolution of the electronic structure of the adsorbate-metal system are derived and the spectrum of electronic excitations is found. The behavior of the model is demonstrated for a set of physically reasonable parameters.

Ab initio simulation of the spin transition during chemisorption: H/Al(111)

Electronically nonadiabatic effects during the chemisorption of hydrogen atoms on an Al(111) surface are simulated ab initio using time-dependent density-functional theory for the electrons in combination with Ehrenfest dynamics for the nuclei. Strongly nonadiabatic effects close to the spin transition of the H atom are identified, and the dissipated energy as well as the electron-hole pair excitation spectra are calculated. The recent Newns-Anderson-model approach by Mizielinski et al. is confirmed. The simulations illustrate the physical processes that contribute to internal exoelectron emission.

Calculations of scattering of N2 molecules from Ru(0001)

Recent experimental measurements of state resolved scattering of nitrogen molecules from a Ru(0001) surface are discussed in comparison with a mixed quantum-classical theory that has given reasonable explanations for similar data on other systems. Acceptable agreement between data and calculations is obtained, but only upon assuming an effective mass of the surface equal to 2.3 times the mass of a single Ru atom.

High translational energy release in H2 (D2) associative desorption from H (D) chemisorbed on C(0001)

Highly energetic translational energy distributions are reported for hydrogen and deuterium molecules desorbing associatively from the atomic chemisorption states on highly oriented pyrolytic graphite (HOPG). Laser assisted associative desorption is used to measure the time of flight of molecules desorbing from a hydrogen (deuterium) saturated HOPG surface produced by atomic exposure from a thermal atom source at around 2100 K. The translational energy distributions normal to the surface are very broad, from approximately 0.5 to approximately 3 eV, with a peak at approximately 1.3 eV. The highest translational energy measured is close to the theoretically predicted barrier height. The angular distribution of the desorbing molecules is sharply peaked along the surface normal and is consistent with thermal broadening contributing to energy release parallel to the surface. All results are in qualitative agreement with recent density functional theory calculations suggesting a lowest energy para-type dimer recombination path.

Femtosecond laser induced associative desorption of H2 from Ru(0001): Comparison of “first principles” theory with experiment

A three dimensional model based on molecular dynamics with electronic frictions is developed to describe the femtosecond laser induced associative desorption of H2 from Ru(0001)(1 x 1)H. Two molecular coordinates (internuclear separation d and center of mass distance to surface z) and a single phonon coordinate are included in the dynamics. Both the potential energy surface and the electronic friction tensor are calculated by density functional theory so that there are no adjustable parameters in the comparison of this model with the wide range of experiments available for this system. This "first principles" dynamic model gives results in semiquantitative agreement with all experimental results; nonlinear fluence dependence of the yield, isotope effect, two pulse correlation, and energy partitioning. The good agreement of theory with experiment supports a description of this surface femtochemistry in terms of thermalized hot electron induced chemistry with coupling to nuclear coordinates through electronic frictions. By comparing the dynamics with the analytical one dimensional frictional model used previously to fit the experiments for this system, we show that the success of the one dimensional model is based on the rapid intermixing of the z and d coordinates as the H-H climbs out of the adsorption well. However, projecting the three dimensional dynamics onto one dimension introduces a fluence (adsorbate temperature) dependent "entropic" barrier in addition to the potential barrier for the chemistry. This implies that some caution must be used in interpreting activation energies obtained in fitting experiments to the one dimensional model.

Reactive and Nonreactive Scattering of H2 from a Metal Surface Is Electronically Adiabatic

The Born-Oppenheimer approximation of uncoupled electronic and nuclear motion is a standard tool of the computational chemist. However, its validity for molecule-metal surface reactions, which are important to heterogeneous catalysis, has been questioned because of the possibility of electron-hole pair excitations. We have performed experiments and calculations on the scattering of molecular hydrogen from a catalytically relevant metal surface, obtaining absolute probabilities for changes in the molecule's velocity parallel to the representative Pt(111) surface. The comparison for in-plane and out-of-plane scattering and results for dissociative chemisorption in the same system show that for hydrogen-metal systems, reaction and diffractive scattering can be accurately described using the Born-Oppenheimer approximation.

Multidimensional effects on dissociation of N2 on Ru(0001)

The applicability of the Born-Oppenheimer approximation to molecule-metal surface reactions is presently a topic of intense debate. We have performed classical trajectory calculations on a prototype activated dissociation reaction, of N2 on Ru(0001), using a potential energy surface based on density functional theory. The computed reaction probabilities are in good agreement with molecular beam experiments. Comparison to previous calculations shows that the rotation of N2 and its motion along the surface affect the reactivity of N2 much more than nonadiabatic effects.

Theoretical evidence for nonadiabatic vibrational deexcitation in H2(D2) state-to-state scattering from Cu(100)

Dynamical calculations are presented for electronically nonadiabatic vibrational deexcitation of H2 and D2 in scattering from Cu(111). Both the potential energy surface and the nonadiabatic coupling strength were obtained from density functional calculations. The theoretically predicted magnitude of the deexcitation and its dependence on incident energy and isotope are all in agreement with state-to-state scattering experiments [on Cu(100)], and this gives indirect evidence for a nonadiabatic mechanism of the observed deexcitation. Direct evidence could be obtained by measuring the chemicurrent associated with the deexcitation, and its properties have been predicted.