Quantum Dynamical Approach to Ultrafast Molecular Desorption from Surfaces

Communication: importance sampling including path correlation in semiclassical initial value representation calculations for time correlation functions

Full semiclassical (SC) initial value representation (IVR) for time correlation functions involves a double phase space average over a set of two phase points, each of which evolves along a classical path. Conventionally, the two initial phase points are sampled independently for all degrees of freedom (DOF) in the Monte Carlo procedure. Here, we present an efficient importance sampling scheme by including the path correlation between the two initial phase points for the bath DOF, which greatly improves the performance of the SC-IVR calculations for large molecular systems. Satisfactory convergence in the study of quantum coherence in vibrational relaxation has been achieved for a benchmark system-bath model with up to 21 DOF.

Photoexcitation of adsorbates on metal surfaces: one-step or three-step

In this essay we discuss the light-matter interactions at molecule-covered metal surfaces that initiate surface photochemistry. The hot-electron mechanism for surface photochemistry, whereby the absorption of light by a metal surface creates an electron-hole pair, and the hot electron scatters through an unoccupied resonance of adsorbate to initiate nuclear dynamics leading to photochemistry, has become widely accepted. Yet, ultrafast spectroscopic measurements of molecule-surface electronic structure and photoexcitation dynamics provide scant support for the hot electron mechanism. Instead, in most cases the adsorbate resonances are excited through photoinduced substrate-to-adsorbate charge transfer. Based on recent studies of the role of coherence in adsorbate photoexcitation, as measured by the optical phase and momentum resolved two-photon photoemission measurements, we examine critically the hot electron mechanism, and propose an alternative description based on direct charge transfer of electrons from the substrate to adsorbate. The advantage of this more quantum mechanically rigorous description is that it informs how material properties of the substrate and adsorbate, as well as their interaction, influence the frequency dependent probability of photoexcitation and ultimately how light can be used to probe and control surface femtochemistry.

Desorption of H2, HD, and D2 through temporary negative ion formation

Physically adsorbed H(2), HD, and D(2) desorbs via the low-energy 2Σ(u)+ electron scattering resonance with characteristic electron energy dependent cross sections as we show by desorption and electron energy-loss measurements. Only H(2) and HD desorb, for energetic reasons, via excitation of the j = 0 → 2 rotational mode and its subsequent decay into free translational motion, while all three isotopes desorb by direct transitions to unbound continuum states, with a probability close to unity and an unexpectedly weak mass dependence ~m(-1/2).

Nonequilibrium thermodynamics of interacting tunneling transport: variational grand potential, density functional formulation and nature of steady-state forces

The standard formulation of tunneling transport rests on an open-boundary modeling. There, conserving approximations to nonequilibrium Green function or quantum statistical mechanics provide consistent but computational costly approaches; alternatively, the use of density-dependent ballistic-transport calculations (e.g., Lang 1995 Phys. Rev. B 52 5335), here denoted 'DBT', provides computationally efficient (approximate) atomistic characterizations of the electron behavior but has until now lacked a formal justification. This paper presents an exact, variational nonequilibrium thermodynamic theory for fully interacting tunneling and provides a rigorous foundation for frozen-nuclei DBT calculations as a lowest-order approximation to an exact nonequilibrium thermodynamic density functional evaluation. The theory starts from the complete electron nonequilibrium quantum statistical mechanics and I identify the operator for the nonequilibrium Gibbs free energy which, generally, must be treated as an implicit solution of the fully interacting many-body dynamics. I demonstrate a minimal property of a functional for the nonequilibrium thermodynamic grand potential which thus uniquely identifies the solution as the exact nonequilibrium density matrix. I also show that the uniqueness-of-density proof from a closely related Lippmann-Schwinger collision density functional theory (Hyldgaard 2008 Phys. Rev. B 78 165109) makes it possible to express the variational nonequilibrium thermodynamic description as a single-particle formulation based on universal electron-density functionals; the full nonequilibrium single-particle formulation improves the DBT method, for example, by a more refined account of Gibbs free energy effects. I illustrate a formal evaluation of the zero-temperature thermodynamic grand potential value which I find is closely related to the variation in the scattering phase shifts and hence to Friedel density oscillations. This paper also discusses the difference between the here-presented exact thermodynamic forces and the often-used electrostatic forces. Finally the paper documents an inherent adiabatic nature of the thermodynamic forces and observes that these are suited for a nonequilibrium implementation of the Born-Oppenheimer approximation.

Stochastic approach to laser-induced ultrafast dynamics: the desorption of H(2)/D(2) from Ru(0001)

The desorption of molecular hydrogen and deuterium induced by femtosecond-laser pulses is studied theoretically for the so-called DIMET (Desorption Induced by Multiple Electronic Transitions) process. These investigations are based on nonadiabatic classical Monte Carlo trajectory (CMCT) simulations on a ground and an excited state potential energy surface, including up to all six adsorbate degrees of freedom. The focus is on the hot-electron mediated energy transfer from the surface to the molecule and back, and the energy partitioning between the different degrees of freedom of the desorbing molecules. We first validate for a two-mode model comprising the desorption mode and the internal vibrational coordinate, the classical Monte Carlo trajectory method by comparing with Monte Carlo wavepacket (MCWP) calculations arising from a fully quantum mechanical open-system density matrix treatment. We then proceed by extending the CMCT calculations to include all six nuclear degrees of freedom of the desorbing molecule. This allows for a detailed comparison between theory and experiment concerning isotope effects, energy partitioning (translational, vibrational, and rotational energies and their distributions), and the dependence of these properties on the laser fluence. The most important findings are as follows. (i) CMCT agrees qualitative with the MCWP scheme. (ii) The basic experimental features such as the large isotope effect, the non-linear increase of yield with laser fluence, translationally hot products (in the order of several 1000 K) and non-equipartitioning of translational and internal energies (E(trans) > E(vib) > E(rot)) are well reproduced. (iii) Predictions concerning a strong angular dependence of translational energies at large observation angles are also made.

Pi resonance of chemisorbed alkali atoms on noble metals

We have performed a joint experimental and theoretical study of the unoccupied electronic structure of alkali adsorbates on the (111) surfaces of Cu and Ag. Combining angle- and time-resolved two-photon photoemission spectroscopy with wave packet propagation calculations we show that, along with the well known sigma resonance oriented along the surface normal, there exist long-lived alkali-localized resonances oriented parallel to the surface (pi symmetry). These new resonances are stabilized by the projected band gap of the substrate and emerge primarily from the mixing of the p and d Rydberg orbitals of the free alkali atom modified by the interaction with the surface.

Temperature-dependent femtosecond photoinduced desorption in CO/Pd(111)

The desorption of CO from a Pd(111) surface following absorption of 120 fs pulses of 780 nm light occurs on two distinct and well-separated time scales. Two-pulse correlation measurements show a fast subpicosecond decay followed by a slower, approximately 40 ps decay. Simulations based on the two-temperature model of electron and phonon heat baths within the substrate, and an empirical friction model to treat coupling to the adsorbate, support the assignment of the desorption mechanism as an electron-mediated process. The photodesorption yield and overall width of the temporal response exhibit a marked dependence on the initial surface temperature in the 100-375 K range despite the much higher transient electronic temperatures (approximately 7000 K) achieved. The observed temperature dependences can be attributed directly to variations in the initial temperature within the frictional coupling picture. Simulations of this extended data set imply that the activation barrier to photoinduced desorption is equal in magnitude to that derived from thermal desorption experiments for this system within the limits of a one-dimensional Arrhenius desorption model. The simulations also imply that the slower decay is not the result of phonon-driven desorption. Though we cannot unambiguously determine the strength of the adsorbate-phonon coupling, our results suggest that its role is to moderate the degree of the adsorbate excitation.

Adsorption-state-dependent subpicosecond photoinduced desorption dynamics

Femtosecond laser excitation has been used to initiate desorption of molecular oxygen from the (111) surface of Pd and to study the adsorption-state dependence of the substrate-adsorbate coupling. The relative populations of the two chemical states, peroxo (O2(2-)) and superoxo (O2-), were varied by changing the total coverage. Two-pulse correlation measurements exhibit a dominant 400 fs response and a slower 10 ps decay that are relatively independent of the initial O2 coverage. In contrast, the photodesorption yield and the nonlinearity of the fluence dependence show a systematic coverage dependence. The coverage-independent subpicosecond response indicates that the photoinduced desorption from the two states is driven primarily by the same electron-mediated mechanism, while the coverage dependence of the yield indicates that the desorption efficiency from the superoxo state is greater than that from the peroxo state. These results are discussed in the context of the electron-phonon two-temperature model with an empirical adsorbate-electron frictional coupling that depends on both the electronic temperature and the activation energy for desorption. With a coupling strength that decreases as the activation energy decreases, the trends with varying coverage, absorbed fluence, and time delay can all be reproduced. The model is consistent with a transition from a resonantly enhanced (diabatic) regime to an adiabatic regime as the system relaxes, accounting for the biexponential correlation behavior.

Ontology of temperature in nonequilibrium systems

The laws of thermodynamics provide a clear concept of the temperature for an equilibrium system in the continuum limit. Meanwhile, the equipartition theorem allows one to make a connection between the ensemble average of the kinetic energy and the uniform temperature. When a system or its environment is far from equilibrium, however, such an association does not necessarily apply. In small systems, the regression hypothesis may not even apply. Herein, it is shown that in small nonequilibrium systems, the regression hypothesis still holds, though with a generalized definition of the temperature. The latter must now be defined for each such manifestation.