Continuum limit semiclassical initial value representation for dissipative systems

In this paper, we consider a dissipative system in which the system is coupled linearly to a harmonic bath. In the continuum limit, the bath is defined via a spectral density and the classical system dynamics is given in terms of a generalized Langevin equation. Using the path integral formulation and factorized initial conditions, it is well known that one can integrate out the harmonic bath, leaving only a path integral over the system degrees of freedom. However, the semiclassical initial value representation treatment of dissipative systems has usually been limited to a discretized treatment of the bath in terms of a finite number of bath oscillators. In this paper, the continuum limit of the semiclassical initial value representation is derived for dissipative systems. As in the path integral, the action is modified with an added nonlocal term, which expresses the influence of the bath on the dynamics. The first order correction term to the semiclassical initial value approximation is also derived in the continuum limit.     

Optimizing hierarchical equations of motion for quantum dissipation and quantifying quantum bath effects on quantum transfer mechanisms

We present an optimized hierarchical equations of motion theory for quantum dissipation in multiple Brownian oscillators bath environment, followed by a mechanistic study on a model donor-bridge-acceptor system. We show that the optimal hierarchy construction, via the memory-frequency decomposition for any specified Brownian oscillators bath, is generally achievable through a universal pre-screening search. The algorithm goes by identifying the candidates for the best be just some selected Padé spectrum decomposition based schemes, together with a priori accuracy control criterions on the sole approximation, the white-noise residue ansatz, involved in the hierarchical construction. Beside the universal screening search, we also analytically identify the best for the case of Drude dissipation and that for the Brownian oscillators environment without strongly underdamped bath vibrations. For the mechanistic study, we quantify the quantum nature of bath influence and further address the issue of localization versus delocalization. Proposed are a reduced system entropy measure and a state-resolved constructive versus destructive interference measure. Their performances on quantifying the correlated system-environment coherence are exemplified in conjunction with the optimized hierarchical equations of motion evaluation of the model system dynamics, at some representing bath parameters and temperatures. Analysis also reveals the localization to delocalization transition as temperature decreases.     

Finite temperature application of the corrected propagator method to reactive dynamics in a condensed-phase environment

The recently proposed mixed quantum-classical method is extended to applications at finite temperatures. The method is designed to treat complex systems consisting of a low-dimensional quantum part (the primary system) coupled to a dissipative bath described classically. The method is based on a formalism showing how to systematically correct the approximate zeroth-order evolution rule. The corrections are defined in terms of the total quantum Hamiltonian and are taken to the classical limit by introducing the frozen Gaussian approximation for the bath degrees of freedom. The evolution of the primary system is governed by the corrected propagator yielding the exact quantum dynamics. The method has been tested on a standard model system describing proton transfer in a condensed-phase environment: a symmetric double-well potential bilinearly coupled to a bath of harmonic oscillators. Flux correlation functions and thermal rate constants have been calculated at two different temperatures for a range of coupling strengths. The results have been compared to the fully quantum simulations of Topaler and Makri [J. Chem. Phys. 101, 7500 (1994)] with the real path integral method.     

Optimal control simulation of the Deutsch-Jozsa algorithm in a two-dimensional double well coupled to an environment

The quantum Deutsch-Jozsa algorithm is implemented by using vibrational modes of a two-dimensional double well. The laser fields realizing the different gates (NOT, CNOT, and HADAMARD) on the two-qubit space are computed by the multitarget optimal control theory. The stability of the performance index is checked by coupling the system to an environment. Firstly, the two-dimensional subspace is coupled to a small number Nb of oscillators in order to simulate intramolecular vibrational energy redistribution. The complete (2+Nb)D problem is solved by the coupled harmonic adiabatic channel method which allows including coupled modes up to Nb=5. Secondly, the computational subspace is coupled to a continuous bath of oscillators in order to simulate a confined environment expected to be favorable to achieve molecular computing, for instance, molecules confined in matrices or in a fullerene. The spectral density of the bath is approximated by an Ohmic law with a cutoff for some hundreds of cm(-1). The time scale of the bath dynamics (of the order of 10 fs) is then smaller than the relaxation time and the controlled dynamics (2 ps) so that Markovian dissipative dynamics is used.     

Stochastic Liouville equations for coherent multidimensional spectroscopy of excitons

Signatures of chemical exchange and spectral diffusion in 2D photon-echo line shapes of molecular aggregates are studied using model calculations for a dimer whose Hamiltonian parameters are stochastically modulated. Cross peaks induced by chemical exchange and by exciton transport have different dynamics and distinguish two models which have the same absorption spectrum (a two-state jump bath modulation model of a dimer and a four-state jump bath model of a single chromophore). Slow Gaussian-Markovian spectral diffusion of a symmetric dimer induces new peaks which are damped as the dipole moment is equilibrated. These effects require an explicit treatment of the bath and may not be described by lower-level theories such as the Redfield equations, which eliminate the bath.     

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.     

Analysis of kinetic isotope effects for nonadiabatic reactions

Factors influencing the rates of quantum mechanical particle transfer reactions in many-body systems are discussed. The investigations are carried out on a simple model for a proton transfer reaction that captures generic features seen in more realistic models of condensed phase systems. The model involves a bistable quantum oscillator coupled to a one-dimensional double-well reaction coordinate, which is in turn coupled to a bath of harmonic oscillators. Reactive-flux correlation functions that involve quantum-classical Liouville dynamics for chemical species operators and quantum equilibrium sampling are used to estimate the reaction rates. Approximate analytical expressions for the quantum equilibrium structure are derived. Reaction rates are shown to be influenced significantly by both the quantum equilibrium structure and nonadiabatic dynamics. Nonadiabatic dynamical effects are found to play the major role in determining the magnitude of the kinetic isotope effect for the model transfer reaction.     

Nonadiabatic chemical dynamics in an intense laser field: electronic wave packet coupled with classical nuclear motions

Dynamics of molecules in an intense laser field is studied in terms of the quantum electronic wave packet coupled with classical nuclear motions. The equations of motion are derived taking a proper account of molecular interactions with the vector potential of a classical electromagnetic field, along with the nonadiabatic interaction due to the breakdown of the Born-Oppenheimer approximation. With the aid of electronic structure calculations, the present method enables us to track, in an ab initio manner, the dynamics of polyatomic molecules in an intense field. Preliminary calculations are carried out for the vibrational state of LiF and a collision of Li+F under an intense laser pulse, which are limited to the domain of no ionization.     

Does ℏ play a role in multidimensional spectroscopy? Reduced hierarchy equations of motion approach to molecular vibrations

To investigate the role of quantum effects in vibrational spectroscopies, we have carried out numerically exact calculations of linear and nonlinear response functions for an anharmonic potential system nonlinearly coupled to a harmonic oscillator bath. Although one cannot carry out the quantum calculations of the response functions with full molecular dynamics (MD) simulations for a realistic system which consists of many molecules, it is possible to grasp the essence of the quantum effects on the vibrational spectra by employing a model Hamiltonian that describes an intra- or intermolecular vibrational motion in a condensed phase. The present model fully includes vibrational relaxation, while the stochastic model often used to simulate infrared spectra does not. We have employed the reduced quantum hierarchy equations of motion approach in the Wigner space representation to deal with nonperturbative, non-Markovian, and nonsecular system-bath interactions. Taking the classical limit of the hierarchy equations of motion, we have obtained the classical equations of motion that describe the classical dynamics under the same physical conditions as in the quantum case. By comparing the classical and quantum mechanically calculated linear and multidimensional spectra, we found that the profiles of spectra for a fast modulation case were similar, but different for a slow modulation case. In both the classical and quantum cases, we identified the resonant oscillation peak in the spectra, but the quantum peak shifted to the red compared with the classical one if the potential is anharmonic. The prominent quantum effect is the 1-2 transition peak, which appears only in the quantum mechanically calculated spectra as a result of anharmonicity in the potential or nonlinearity of the system-bath coupling. While the contribution of the 1-2 transition is negligible in the fast modulation case, it becomes important in the slow modulation case as long as the amplitude of the frequency fluctuation is small. Thus, we observed a distinct difference between the classical and quantum mechanically calculated multidimensional spectra in the slow modulation case where spectral diffusion plays a role. This fact indicates that one may not reproduce the experimentally obtained multidimensional spectrum for high-frequency vibrational modes based on classical molecular dynamics simulations if the modulation that arises from surrounding molecules is weak and slow. A practical way to overcome the difference between the classical and quantum simulations was discussed.     

Collisional line shapes for low frequency vibrations of adsorbates on a metal surface

The dynamics of atoms or molecules adsorbed on a metal surface, and excited by collisions with an atomic beam, are treated within a theory that includes energy dissipation into lattice vibrations by means of a frequency and temperature dependent friction function. The theory provides dynamic structure factors for energy transfer derived from collisional time correlation functions. It describes the relaxation of a vibrationally excited atom or molecule within a model of a damped quantum harmonic oscillator bilinearly coupled to a bath of lattice oscillators. The collisional time correlation function is generalized to include friction effects and is applied to the vibrational relaxation of the frustrated translation mode of Na adsorbed on a Cu(001) surface, CO on Cu(001), and CO on Pt(111), following excitation by collisions with He atoms. Results for the frequency shift and width of line shapes versus surface temperature are in very good agreement with experimental measurements of inelastic He atom scattering. Our interpretation of the experimental results provides insight on the relative role of phonon versus electron-hole relaxation.     

The femtosecond dynamics of electron transfer in a modified photosynthesis reaction center: Quantum, classical, and kinetic analyses

The dynamics of electron transfer in a modified photosynthesis reaction center in which electron transfer from the bridge to the acceptor is blocked is considered. A microscopic model of the process is suggested. Within this model, the diabatic electronic states of the donor and bridge are described by one-dimensional displaced harmonic oscillators. The dynamics of the population of electronic states is calculated by the quantum method of wave packets and classical and kinetic modeling. The suggested model is used to study the qualitative dependence of the dynamics of electron transfer on the nonadiabatic interaction potential. The parameters of the model are determined by comparing the experimental and calculated absorption spectra of the product of electron transfer. It is shown that kinetic models can be used to approximately describe the dynamics of electron transfer in reaction centers. The boundaries of the applicability of the kinetic method are considered.     

Computer simulation of quantum dynamics in a classical spin environment

In this paper, a formalism for studying the dynamics of quantum systems coupled to classical spin environments is reviewed. The theory is based on generalized antisymmetric brackets and naturally predicts open-path off-diagonal geometric phases in the evolution of the density matrix. It is shown that such geometric phases must also be considered in the quantum–classical Liouville equation for a classical bath with canonical phase space coordinates; this occurs whenever the adiabatics basis is complex (as in the case of a magnetic field coupled to the quantum subsystem). When the quantum subsystem is weakly coupled to the spin environment, non-adiabatic transitions can be neglected and one can construct an effective non-Markovian computer simulation scheme for open quantum system dynamics in classical spin environments. In order to tackle this case, integration algorithms based on the symmetric Trotter factorization of the classical-like spin propagator are derived. Such algorithms are applied to a model comprising a quantum two-level system coupled to a single classical spin in an external magnetic field. Starting from an excited state, the population difference and the coherences of this two-state model are simulated in time while the dynamics of the classical spin is monitored in detail. It is the author’s opinion that the numerical evidence provided in this paper is a first step toward developing the simulation of quantum dynamics in classical spin environments into an effective tool. In turn, the ability to simulate such a dynamics can have a positive impact on various fields, among which, for example, nanoscience.     

Dynamics of linear exchange reactions. Quasi-classical model of high-energy vibrational inversion

The dynamics of exothermic exchange reactions in natural reaction coordinates (NRC) is investigated analytically. Such a study is possible because (1) there is a small parameter in several chemical systems, i.e., the product of the curvature of the reaction path and the classical amplitude of transversal vibrations, and (2) the motion along the reaction coordinate may be considered quasi-classically. The quasi-classical equations of motion are derived, which have an exact solution. The transition probabilities between vibrational levels are investigated in detail. The quasi-classical analytic distributions are in satisfactory agreement with the numerical results of quantum calculations made for model potential surfaces. The quasi-classical model is shown to fit for real reaction systems with mass ratio L + HH (L is light atom, and H is heavy atom). The main factors which determine the vibrational distribution of products for these reactions are the curvature of the reaction path and the number of open channels at the point of maximum curvature.