Surrogate Hamiltonian study of electronic relaxation in the femtosecond laser induced desorption of NO/NiO(100)

Wave Function Realization of a Thermal Collision Model

An efficient algorithm to simulate dynamics of open quantum system is presented. The method describes the dynamics by unraveling stochastic wave functions converging to a density operator description. The stochastic techniques are based on the quantum collision model. Modeling systems dynamics with wave functions and modeling the interaction with the environment with a collision sequence reduces the scale of the complexity significantly. The algorithm developed can be implemented on quantum computers. We introduce stochastic methods that exploit statistical characteristics of the model such as Markovianity, Brownian motion, and binary distribution. The central limit theorem is employed to study the convergence of distributions of stochastic dynamics of pure quantum states represented by wave vectors. By averaging a sample of functions in the distribution we prove and demonstrate the convergence of the dynamics to the mixed quantum state described by a density operator.

Dissipative quantum dynamics using the stochastic surrogate Hamiltonian approach

Based on the original Stochastic Surrogate Hamiltonian approach for the simulation of open quantum system dynamics, we present a modified algorithm for the swap procedure. The implementation is tested with nuclear relaxation dynamics of model systems known from literature, i.e., the harmonic oscillator and the Morse oscillator. Finally, the stochastic surrogate Hamiltonian is applied to simulate the dynamical electronic excitation and relaxation of a photodesorption process. This is the first application of the stochastic surrogate Hamiltonian in an ab initio context. A comparison to a surrogate Hamiltonian benchmark allows us to evaluate the results obtained. For this purpose, the well-studied laser-induced desorption of NO from NiO(100) is chosen.

Non-stochastic matrix Schrödinger equation for open systems

We propose an extension of the Schrödinger equation for a quantum system interacting with environment. This extension describes dynamics of a collection of auxiliary wavefunctions organized as a matrix m, from which the system density matrix can be reconstructed as ρ̂=mm(†). We formulate a compatibility condition, which ensures that the reconstructed density satisfies a given quantum master equation for the system density. The resulting non-stochastic evolution equation preserves positive-definiteness of the system density and is applicable to both Markovian and non-Markovian system-bath treatments. Our formalism also resolves a long-standing problem of energy loss in the time-dependent variational principle applied to mixed states of closed systems.

Quantum dynamical study of femtosecond photodesorption of CO from TiO2(110)

The photodesorption of CO from TiO2(110) by femtosecond pulses is investigated with the Surrogate Hamiltonian approach. The aim of the study is to resolve the relaxation mechanism and forecast the lifetime of the exited state based on a microscopic description of the excitation and relaxation processes. The parameters characterizing the system are obtained from ab initio and Density Functional Theory-calculations with one parameter estimated from physical considerations and convergence studies. Two electronic states are considered and the relaxation is assumed to be due to the interaction of the excited adsorbate with electron hole pairs in the surface. Desorption probabilities and velocity distributions of the desorbing molecules are calculated and an exited state lifetime is predicted. Throughout this paper atomic units, i.e., ℏ = me = e = a0 = 1, have been used unless otherwise stated.

Theoretical study of the photodesorption mechanism of nitric oxide on a Ag(111) surface: A nonequilibrium Green’s function approach to hot-electron tunneling

The photoinduced desorption of NO molecules on a Ag surface was studied theoretically using a recently developed method based on the nonequilibrium Green's function approach combined with the density functional theory. Geometry optimizations for the stable NO dimer phase were carried out, and two structures of adsorbed dimers were identified. We calculated the reaction probabilities as a function of incident photon energy for each of the dimers and compared them with experimental action spectra. The two main features of the action spectra, (i) a long tail to the long wavelength (approximately 600 nm) and (ii) a rapid increase at approximately 350 nm, were well reproduced. By theoretical analysis, we found the importance of quantum interference for the interfacial charge transfer between the metal substrate and the adsorbate, as well as the contribution of secondary electrons. Our calculations suggest that the photoactive species is dimeric and that the resonant level is single for the photodesorption of NO.

Two-dimensional surrogate Hamiltonian investigation of laser-induced desorption of NO∕NiO(100)

The photodesorption of NO from NiO(100) is studied from first principles, with electronic relaxation treated by the use of the surrogate Hamiltonian approach. Two nuclear degrees of freedom of the adsorbate-substrate system are taken into account. To perform the quantum dynamical wave-packet calculations, a massively parallel implementation with a one-dimensional data decomposition had to be introduced. The calculated desorption probabilities and velocity distributions are in qualitative agreement with experimental data. The results are compared to those of stochastic wave-packet calculations where a sufficiently large number of quantum trajectories is propagated within a jumping wave-packet scenario.

Dissipative quantum dynamics with the surrogate Hamiltonian approach. A comparison between spin and harmonic baths

The dissipative quantum dynamics of an anharmonic oscillator coupled to a bath is studied with the purpose of elucidating the differences between the relaxation to a spin bath and to a harmonic bath. Converged results are obtained for the spin bath by the surrogate Hamiltonian approach. This method is based on constructing a system-bath Hamiltonian, with a finite but large number of spin bath modes, that mimics exactly a bath with an infinite number of modes for a finite time interval. Convergence with respect to the number of simultaneous excitations of bath modes can be checked. The results are compared to calculations that include a finite number of harmonic modes carried out by using the multiconfiguration time-dependent Hartree method of Nest and Meyer [J. Chem. Phys. 119, 24 (2003)]. In the weak coupling regime, at zero temperature and for small excitations of the primary system, both methods converge to the Markovian limit. When initially the primary system is significantly excited, the spin bath can saturate restricting the energy acceptance. An interaction term between bath modes that spreads the excitation eliminates the saturation. The loss of phase between two cat states has been analyzed and the results for the spin and harmonic baths are almost identical. For stronger couplings, the dynamics induced by the two types of baths deviate. The accumulation and degree of entanglement between the bath modes have been characterized. Only in the spin bath the dynamics generate entanglement between the bath modes.