Vibrational energy and phase relaxation at surfaces

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.

Isotopic effects in vibrational relaxation dynamics of H on a Si(100) surface

In a recent paper [U. Lorenz and P. Saalfrank, Chem. Phys. 482, 69 (2017)], we proposed a robust scheme to set up a system-bath model Hamiltonian, describing the coupling of adsorbate vibrations (system) to surface phonons (bath), from first principles. The method is based on an embedded cluster approach, using orthogonal coordinates for system and bath modes, and an anharmonic phononic expansion of the system-bath interaction up to second order. In this contribution, we use this model Hamiltonian to calculate vibrational relaxation rates of H-Si and D-Si bending modes, coupled to a fully H(D)-covered Si(100)-(2×1) surface, at zero temperature. The D-Si bending mode has an anharmonic frequency lying inside the bath frequency spectrum, whereas the H-Si bending mode frequency is outside the bath Debye band. Therefore, in the present calculations, we only take into account one-phonon system-bath couplings for the D-Si system and both one- and two-phonon interaction terms in the case of H-Si. The computation of vibrational lifetimes is performed with two different approaches, namely, Fermi's golden rule, and a generalized Bixon-Jortner model built in a restricted vibrational space of the adsorbate-surface zeroth-order Hamiltonian. For D-Si, the Bixon-Jortner Hamiltonian can be solved by exact diagonalization, serving as a benchmark, whereas for H-Si, an iterative scheme based on the recursive residue generation method is applied, with excellent convergence properties. We found that the lifetimes obtained with perturbation theory, albeit having almost the same order of magnitude-a few hundred fs for D-Si and a couple of ps for H-Si-, are strongly dependent on the discretized numerical representation of the bath spectral density. On the other hand, the Bixon-Jortner model is free of such numerical deficiencies, therefore providing better estimates of vibrational relaxation rates, at a very low computational cost. The results obtained with this model clearly show a net exponential decay of the time-dependent survival probability for the H-Si initial vibrational state, allowing an easy extraction of the bending mode "lifetime." This is in contrast with the D-Si system, whose survival probability exhibits a non-monotonic decay, making it difficult to define such a lifetime. This different behavior of the vibrational decay is rationalized in terms of the power spectrum of the adsorbate-surface system. In the case of D-Si, it consists of several, non-uniformly distributed peaks around the bending mode frequency, whereas the H-Si spectrum exhibits a single Lorentzian lineshape, whose width corresponds to the calculated lifetime. The present work gives some insight into mechanisms of vibration-phonon coupling at surfaces. It also serves as a benchmark for multidimensional system-bath quantum dynamics, for comparison with approximate schemes such as reduced, open-system density matrix theory (where the bath is traced out and a Liouville-von Neumann equation is solved) or approximate wavefunction methods to solve the combined system-bath Schrödinger equation.

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.

Optimal control of open quantum systems applied to the photochemistry of surfaces

A quantum system in a condensed phase undergoes strong dissipative processes. The last decades have seen the rise of experimental and theoretical approaches for gaining control over dissipative phenomena. From a theoretical viewpoint it is important to model such processes in a rigorous way. An efficient and accurate method to find control fields is optimal control theory (OCT). In this Letter, a control scheme relying on OCT with time-dependent targets is employed to minimize dissipation, modeled within the surrogate Hamiltonian approach, on adsorbate-surface systems.

Optimal control in a dissipative system: Vibrational excitation of CO∕Cu(100) by IR pulses

The question as to whether state-selective population of molecular vibrational levels by shaped infrared laser pulses is possible in a condensed phase environment is of central importance for such diverse fields as time-resolved spectroscopy, quantum computing, or "vibrationally mediated chemistry." This question is addressed here for a model system, representing carbon monoxide adsorbed on a Cu(100) surface. Three of the six vibrational modes are considered explicitly, namely, the CO stretch vibration, the CO-surface vibration, and a frustrated translation. Optimized infrared pulses for state-selective excitation of "bright" and "dark" vibrational levels are designed by optimal control theory in the framework of a Markovian open-system density matrix approach, with energy flow to substrate electrons and phonons, phase relaxation, and finite temperature accounted for. The pulses are analyzed by their Husimi "quasiprobability" distribution in time-energy space.