A complete quantum description of an ultrafast pump-probe charge transfer event in condensed phase

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.

Optimized sampling of mixed-state observables

Quantum dynamical simulations of statistical ensembles pose a significant computational challenge due to the fact that mixed states need to be represented. If the underlying dynamics is fully unitary, for example, in ultrafast coherent control at finite temperatures, then one approach to approximate time-dependent observables is to sample the density operator by solving the Schrödinger equation for a set of wave functions with randomized phases. We show that, on average, random-phase wave functions perform well for ensembles with high mixedness, whereas at higher purities a deterministic sampling of the energetically lowest-lying eigenstates becomes superior. We prove that minimization of the worst-case error for computing arbitrary observables is uniquely attained by eigenstate-based sampling. We show that this error can be used to form a qualitative estimate of the set of ensemble purities for which the sampling performance of the eigenstate-based approach is superior to random-phase wave functions. Furthermore, we present refinements to both schemes which remove redundant information from the sampling procedure to accelerate their convergence. Finally, we point out how the structure of low-rank observables can be exploited to further improve eigenstate-based sampling schemes.

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.

Quantum thermodynamics and open-systems modeling

A comprehensive approach to modeling open quantum systems consistent with thermodynamics is presented. The theory of open quantum systems is employed to define system bath partitions. The Markovian master equation defines an isothermal partition between the system and bath. Two methods to derive the quantum master equation are described: the weak coupling limit and the repeated collision model. The role of the eigenoperators of the free system dynamics is highlighted, in particular, for driven systems. The thermodynamical relations are pointed out. Models that lead to loss of coherence, i.e., dephasing are described. The implication of the laws of thermodynamics to simulating transport and spectroscopy is described. The indications for self-averaging in large quantum systems and thus its importance in modeling are described. Basic modeling by the surrogate Hamiltonian is described, as well as thermal boundary conditions using the repeated collision model and their use in the stochastic surrogate Hamiltonian. The problem of modeling with explicitly time dependent driving is analyzed. Finally, the use of the stochastic surrogate Hamiltonian for modeling ultrafast spectroscopy and quantum control is reviewed.

Exploiting Non-Markovianity for Quantum Control

Quantum technology, exploiting entanglement and the wave nature of matter, relies on the ability to accurately control quantum systems. Quantum control is often compromised by the interaction of the system with its environment since this causes loss of amplitude and phase. However, when the dynamics of the open quantum system is non-Markovian, amplitude and phase flow not only from the system into the environment but also back. Interaction with the environment is then not necessarily detrimental. We show that the back-flow of amplitude and phase can be exploited to carry out quantum control tasks that could not be realized if the system was isolated. The control is facilitated by a few strongly coupled, sufficiently isolated environmental modes. Our paradigmatic example considers a weakly anharmonic ladder with resonant amplitude control only, restricting realizable operations to SO(N). The coupling to the environment, when harnessed with optimization techniques, allows for full SU(N) controllability.

Dynamics of ground state absorption spectra in donor-acceptor pairs with ultrafast charge recombination

A theoretical approach to simulation of the transient spectra in molecular systems with ultrafast photoinduced nonradiative electronic transitions is developed. The evolution of the excited and ground state populations as well as the nonradiative transitions between them are calculated in the framework of the stochastic multichannel point-transition model involving the reorganization of the medium and the intramolecular high frequency vibrational modes. Simulations of transient spectra of donor-acceptor pairs excited in the charge-transfer band that are accompanied by ultrafast charge recombination into the ground state demonstrate a possibility of positive band appearance in the transient absorption spectrum caused by those systems in the ground state, which returned there from the excited state. The region of the parameters of the donor-acceptor systems where a positive ground state absorption signal can be observed is discussed. A qualitative comparison of the simulated transient spectra with the experimental data on betaine-30 is presented.

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.

Strong-pump strong-probe spectroscopy: effects of higher excited electronic states

The present paper is devoted to the simulation of (integral and dispersed) pump-probe signals in the nonperturbative regime for a series of material systems with multiple electronic states and excited-state absorption. We show that strong-pump strong-probe spectroscopy permits the probing of vibrational wavepackets in high-lying and/or short-lived excited electronic states with a time resolution which is not limited by the pulse durations. The field strength can be regarded as an additional experimentally controllable parameter, which can be tuned to maximize the spectroscopic information for a given material system.

Between superexchange and hopping: an intermediate charge-transfer mechanism in poly(A)-poly(T) DNA hairpins

We developed a model for hole migration along relatively short DNA hairpins with fewer that seven adenine (A):thymine (T) base pairs. The model was used to simulate hole migration along poly(A)-poly(T) sequences with a particular emphasis on the impact of partial hole localization on the different rate processes. The simulations, performed within the framework of the stochastic surrogate Hamiltonian approach, give values for the arrival rate in good agreement with experimental data. Theoretical results obtained for hairpins with fewer than three A:T base pairs suggest that hole transfer along short hairpins occurs via superexchange. This mechanism is characterized by the exponential distance dependence of the arrival rate on the donor/acceptor distance, k(a) ≃ e(-βR), with β = 0.9 Å(-1). For longer systems, up to six A:T pairs, the distance dependence follows a power law k(a) ≃ R(-η) with η = 2. Despite this seemingly clear signature of unbiased hopping, our simulations show the complete delocalization of the hole density along the entire hairpin. According to our analysis, the hole transfer along relatively long sequences may proceed through a mechanism which is distinct from both coherent single-step superexchange and incoherent multistep hopping. The criterion for the validity of this mechanism intermediate between superexchange and hopping is proposed. The impact of partial localization on the rate of hole transfer between neighboring A bases was also investigated.

Optical N-wave-mixing spectroscopy with strong and temporally well-separated pulses: the doorway-window representation

We have extended the doorway-window representation of optical pump-probe spectroscopy with weak pulses toward N-wave-mixing spectroscopy with temporally well-separated pulses of arbitrary strength. The expressions for the signals in the strong-pulse doorway-window representation are derived in the framework of the nonperturbative theory of N-wave-mixing spectroscopy. The strong-pulse doorway-window representation is complementary to the equation-of-motion phase-matching approach. The latter fully accounts for pulse-overlap effects in signals induced by weak pulses but is computationally more expensive. The performance of the doorway-window approximation for temporally well-separated strong pulses is illustrated for an electronic two-level system with an underdamped Condon-active vibrational mode.

A nonperturbative calculation of nonlinear spectroscopic signals in liquid solution

Nonlinear spectroscopic signals in liquid solution were calculated without treating the field-matter interaction in a perturbative manner. The calculation is based on the assumption that the intermolecular degrees of freedom can be treated classically, while the time evolution of the electronic state is treated quantum mechanically. The calculated overall electronic polarization is then resolved into its directional components via the method of Seidner et al. [J. Chem. Phys. 103, 3998 (1995)]. It is shown that the time dependence of the directional components is independent of laser intensity in the impulsive pulse regime, which allows for flexibility in choosing the procedure for calculating optical response functions. The utility and robustness of the nonperturbative procedure is demonstrated in the case of a two-state chromophore solvated in a monoatomic liquid, by calculating nonlinear time-domain signals in the strong-field, weak-field, impulsive, and nonimpulsive regimes.

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.

Minimizing broadband excitation under dissipative conditions

Optimal control theory is employed for the task of minimizing the excited-state population of a dye molecule in solution. The spectrum of the excitation pulse is contained completely in the absorption band of the molecule. Only phase control is studied which is equivalent to optimizing the transmission of the pulse through the medium. The molecular model explicitly includes two electronic states and a single vibrational mode. The other degrees of freedom are classified as bath modes. The surrogate Hamiltonian method is employed to incorporate these bath degrees of freedom. Their influence can be classified as electronic dephasing and vibrational relaxation. In accordance with experimental results, minimal excitation is associated with a negatively chirped pulses. Optimal pulses with more complex transient structure are found to be superior to linearly chirped pulses. The difference is enhanced when the fluence is increased. The improvement degrades when dissipative effects become more dominant.

Dissipative dynamics of a system passing through a conical intersection: Ultrafast pump-probe observables

The dynamics of a system incorporating a conical intersection, in the presence of a dissipative environment, is studied with the purpose of identifying observable ultrafast spectroscopic signatures. A model system consisting of two vibronically coupled electronic states with two nuclear degrees of freedom is constructed. Dissipation is treated by two different methods, Lindblad semigroup formalism and the surrogate Hamiltonian approach. Pump-probe experimental expectation values such as transient emission and transient absorption are calculated and compared to the adiabatic and diabatic population transfer. The ultrafast population transfer reflecting the conical intersection is not mirrored in transient absorption measurements such as the recovery of the bleach. Emission from the excited state can be suppressed on the ultrafast time scale, but the existence of a conical intersection is only one of the possible mechanisms that can provide ultrafast damping of emission.

Transient Phenomena in Time- and Frequency-Gated Spontaneous Emission

The effect of overlapping pump and gate pulses on time- and frequency-gated spontaneous emission spectra is explored for a model of material dynamics that accounts for strong nonadiabatic and electron-vibrational coupling effects, vibrational relaxation, and optical dephasing, thus representing characteristic features of photoinduced excited-state dynamics in large molecules in the gas phase or in condensed phases. The behaviors of the sequential, coherent, and doorway-window contributions to the spontaneous emission spectrum are studied separately. The interrelation between the sequential and coherent contributions is demonstrated to be sensitive to the carrier frequencies of the pump and gate pulses and also to the optical dephasing rate, opening the possibility of an experimental determination of the latter. The coherent contribution is shown to dominate the spectrum at specific emission frequencies.

Time- and frequency-resolved fluorescence spectra of nonadiabatic dissipative systems: What photons can tell us

The monitoring of the excited-state dynamics by time- and frequency-resolved spontaneous emission spectroscopy has been studied in detail for a model exhibiting an excited-state curve crossing. The model represents characteristic aspects of the photoinduced ultrafast dynamics in large molecules in the gas or condensed phases and accounts for strong nonadiabatic and electron-vibrational coupling effects, as well as for vibrational relaxation and optical dephasing. A comprehensive overview of the dependence of spontaneous emission spectra on the characteristics of the excitation and detection processes (such as carrier frequencies, pump/gate pulse durations, as well as optical dephasing) is presented. A systematic comparison of ideal spectra, which provide simultaneously perfect time and frequency resolution and thus contain maximal information on the system dynamics, with actually measurable time- and frequency-gated spectra has been carried out. The calculations of real time- and frequency-gated spectra demonstrate that complementary information on the excited-state dynamics can be extracted when the duration of the gate pulse is varied.

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.

Femtosecond photodesorption of small molecules from surfaces: a theoretical investigation from first principles

A microscopic model for the excitation and relaxation processes in photochemistry at surfaces is developed. Our study is based on ab initio calculations and the surrogate Hamiltonian method treating surface electron-hole pairs as a bath of two-level systems. Desorption probabilities and velocities in the experimentally observed range are obtained. The excited state lifetime is calculated, and a dependence of observables on pulse length is predicted.