Surface-hopping dynamics and decoherence with quantum equilibrium structure

In open quantum systems, decoherence occurs through interaction of a quantum subsystem with its environment. The computation of expectation values requires a knowledge of the quantum dynamics of operators and sampling from initial states of the density matrix describing the subsystem and bath. We consider situations where the quantum evolution can be approximated by quantum-classical Liouville dynamics and examine the circumstances under which the evolution can be reduced to surface-hopping dynamics, where the evolution consists of trajectory segments exclusively evolving on single adiabatic surfaces, with probabilistic hops between these surfaces. The justification for the reduction depends on the validity of a Markovian approximation on a bath averaged memory kernel that accounts for quantum coherence in the system. We show that such a reduction is often possible when initial sampling is from either the quantum or classical bath initial distributions. If the average is taken only over the quantum dispersion that broadens the classical distribution, then such a reduction is not always possible.     

非平衡体系-热库纠缠定理与热输运

将最近建立的体系-热库纠缠定理(SBET)扩展到非平衡的情形. 其中, 任意体系与处于不同温度的多个高斯型热库环境相耦合. 现有的SBET将体系-热库的纠缠响应函数与体系的局域响应函数联系起来, 而扩展的理论则关注通过分子结的非平衡稳态量子输运流. 新理论是基于广义Langevin方程建立的, 它与量子情形下的非平衡热力学密切相关.     

Proton transfer reactions in model condensed-phase environments: Accurate quantum dynamics using the multilayer multiconfiguration time-dependent Hartree approach

The recently proposed multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) approach to evaluating reactive quantum dynamics is applied to two model condensed-phase proton transfer reactions. The models consist of a one-dimensional double-well "system" that is bilinearly coupled to a "bath" of harmonic oscillators parameterized to represent a condensed-phase environment. Numerically exact quantum-mechanical flux correlation functions and thermal rate constants are obtained for a broad range of temperatures and system-bath coupling strengths, thus demonstrating the efficacy of the ML-MCTDH approach. Particular attention is focused on the regime where low temperatures are combined with weak system-bath coupling. Under such conditions it is found that long propagation times are often required and that quantum coherence effects may prevent a rigorous determination of the rate constant.     

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.     

Condensed-phase relaxation of multilevel quantum systems. I. An exactly solvable model

An analytically solvable model of multilevel condensed-phase quantum dynamics relevant to vibrational relaxation and electron transfer is presented. Exact solutions are derived for the reduced system density matrix dynamics of a degenerate N-level quantum system characterized by nearest-neighbor hopping and off-diagonal coupling (which is linear in the bath coordinates) to a harmonic oscillator bath. We demonstrate that for N> 2 the long-time steady-state system site occupation probabilities are not the same for all sites; that is, they are distributed in a non-Boltzmann manner, which depends on the initial conditions and the number of levels in the system. Although the system-bath Hamiltonian considered here is restricted in form, the availability of an exact solution enables us to study the model in all regions of an extensive parameter space.     

Second-order quantized Hamilton dynamics coupled to classical heat bath

Starting with a quantum Langevin equation describing in the Heisenberg representation a quantum system coupled to a quantum bath, the Markov approximation and, further, the closure approximation are applied to derive a semiclassical Langevin equation for the second-order quantized Hamilton dynamics (QHD) coupled to a classical bath. The expectation values of the system operators are decomposed into products of the first and second moments of the position and momentum operators that incorporate zero-point energy and moderate tunneling effects. The random force and friction as well as the system-bath coupling are decomposed to the lowest classical level. The resulting Langevin equation describing QHD-2 coupled to classical bath is analyzed and applied to free particle, harmonic oscillator, and the Morse potential representing the OH stretch of the SPC-flexible water model.     

A local coherent-state approximation to system-bath quantum dynamics

A novel quantum method to deal with typical system-bath dynamical problems is introduced. Subsystem discrete variable representation and bath coherent-state sets are used to write down a multiconfigurational expansion of the wave function of the whole system. With the help of the Dirac-Frenkel variational principle, simple equations of motion--a kind of Schrodinger-Langevin equation for the subsystem coupled to (pseudo) classical equations for the bath--are derived. True dissipative dynamics at all times is obtained by coupling the bath to a secondary, classical Ohmic bath, which is modeled by adding a friction coefficient in the derived pseudoclassical bath equations. The resulting equations are then solved for a number of model problems, ranging from tunneling to vibrational relaxation dynamics. Comparison of the results with those of exact, multiconfiguration time-dependent Hartree calculations in systems with up to 80 bath oscillators shows that the proposed method can be very accurate and might be of help in studying realistic problems with very large baths. To this end, its linear scaling behavior with respect to the number of bath degrees of freedom is shown in practice with model calculations using tens of thousands of bath oscillators.     

Multidimensional vibrational spectroscopy for tunneling processes in a dissipative environment

Simulating tunneling processes as well as their observation are challenging problems for many areas. In this study, we consider a double-well potential system coupled to a heat bath with a linear-linear (LL) and square-linear (SL) system-bath interactions. The LL interaction leads to longitudinal (T1) and transversal (T2) homogeneous relaxations, whereas the SL interaction leads to the inhomogeneous dephasing (T2*) relaxation in the white noise limit with a rotating wave approximation. We discuss the dynamics of the double-well system under infrared (IR) laser excitations from a Gaussian-Markovian quantum Fokker-Planck equation approach, which was developed by generalizing Kubo's stochastic Liouville equation. Analytical expression of the Green function is obtained for a case of two-state-jump modulation by performing the Fourier-Laplace transformation. We then calculate a two-dimensional infrared signal, which is defined by the four-body correlation function of optical dipole, for various noise correlation time, system-bath coupling parameters, and temperatures. It is shown that the bath-induced vibrational excitation and relaxation dynamics between the tunneling splitting levels can be detected as the isolated off-diagonal peaks in the third-order two-dimensional infrared (2D-IR) spectroscopy for a specific phase matching condition. Furthermore, this spectroscopy also allows us to directly evaluate the rate constants for tunneling reactions, which relates to the coherence between the splitting levels; it can be regarded as a novel technique for measuring chemical reaction rates. We depict the change of reaction rates as a function of system-bath coupling strength and a temperature through the 2D-IR signal.     

Theory of multichromophoric coherent resonance energy transfer: a polaronic quantum master equation approach

The approach of second order time local quantum master equation in the polaron picture, which has been employed for a theory of coherent resonance energy transfer, is extended for general multichromophore systems. Explicit expressions for all the kernel and inhomogeneous terms are derived, which can be calculated by any standard numerical procedure. The theory is then applied to a model of donor-bridge-acceptor system moderately coupled to bosonic bath. The results are compared with those based on the theory of Fo?rster's resonance energy transfer. It is shown that coherently coupled multichromophores can speed up the transfer of energy substantially and in a way insensitive to the disorder.     

Time-dependent density functional theory of open quantum systems in the linear-response regime

Time-dependent density functional theory (TDDFT) has recently been extended to describe many-body open quantum systems evolving under nonunitary dynamics according to a quantum master equation. In the master equation approach, electronic excitation spectra are broadened and shifted due to relaxation and dephasing of the electronic degrees of freedom by the surrounding environment. In this paper, we develop a formulation of TDDFT linear-response theory (LR-TDDFT) for many-body electronic systems evolving under a master equation, yielding broadened excitation spectra. This is done by mapping an interacting open quantum system onto a noninteracting open Kohn-Sham system yielding the correct nonequilibrium density evolution. A pseudoeigenvalue equation analogous to the Casida equations of the usual LR-TDDFT is derived for the Redfield master equation, yielding complex energies and Lamb shifts. As a simple demonstration, we calculate the spectrum of a C(2 +) atom including natural linewidths, by treating the electromagnetic field vacuum as a photon bath. The performance of an adiabatic exchange-correlation kernel is analyzed and a first-order frequency-dependent correction to the bare Kohn-Sham linewidth based on the Go?rling-Levy perturbation theory is calculated.     

Hybrid quantum/classical molecular dynamics for a proton transfer reaction coupled to a dissipative bath

A hybrid quantum/classical molecular dynamics approach is applied to a proton transfer reaction represented by a symmetric double well system coupled to a dissipative bath. In this approach, the proton is treated quantum mechanically and all bath modes are treated classically. The transition state theory rate constant is obtained from the potential of mean force, which is generated along a collective reaction coordinate with umbrella sampling techniques. The transmission coefficient, which accounts for dynamical recrossings of the dividing surface, is calculated with a reactive flux approach combined with the molecular dynamics with quantum transitions surface hopping method. The hybrid quantum/classical results agree well with numerically exact results in the spatial-diffusion-controlled regime, which is most relevant for proton transfer in proteins. This hybrid quantum/classical approach has already been shown to be computationally practical for studying proton transfer in large biological systems. These results have important implications for future applications to hydrogen transfer reactions in solution and proteins.     

A derivation of the mixed quantum-classical Liouville equation from the influence functional formalism

We show that the mixed quantum-classical Liouville equation is equivalent to linearizing the forward-backward action in the influence functional. Derivations are provided in terms of either the diabatic or adiabatic basis sets. An application of the mixed quantum-classical Liouville equation for calculating the memory kernel of the generalized quantum master equation is also presented. The accuracy and computational feasibility of such an approach is demonstrated in the case of a two-level system nonlinearly coupled to an anharmonic bath.     

Reduced density matrix hybrid approach: an efficient and accurate method for adiabatic and non-adiabatic quantum dynamics

We present a new approach to calculate real-time quantum dynamics in complex systems. The formalism is based on the partitioning of a system's environment into "core" and "reservoir" modes with the former to be treated quantum mechanically and the latter classically. The presented method only requires the calculation of the system's reduced density matrix averaged over the quantum core degrees of freedom which is then coupled to a classically evolved reservoir to treat the remaining modes. We demonstrate our approach by applying it to the spin-boson problem using the noninteracting blip approximation to treat the system and core, and Ehrenfest dynamics to treat the reservoir. The resulting hybrid methodology is accurate for both fast and slow baths, since it naturally reduces to its composite methods in their respective regimes of validity. In addition, our combined method is shown to yield good results in intermediate regimes where neither approximation alone is accurate and to perform equally well for both strong and weak system-bath coupling. Our approach therefore provides an accurate and efficient methodology for calculating quantum dynamics in complex systems.     

Exact Generator and its High Order Expansions in Time-Convolutionless Generalized Master Equation: Applications to Spin-Boson Model and Exictation Energy Transfer?

The time-convolutionless (TCL) quantum master equation provides a powerful tool to simulate reduced dynamics of a quantum system coupled to a bath. The key quantity in the TCL master equation is the so-called kernel or generator, which describes effects of the bath degrees of freedom. Since the exact TCL generators are usually hard to calculate analytically, most applications of the TCL generalized master equation have relied on approximate generators using second and fourth order perturbative expansions. By using the hierarchical equation of motion (HEOM) and extended HEOM methods, we present a new approach to calculating the exact TCL generator and its high order perturbative expansions. The new approach is applied to the spin-boson model with different sets of parameters, to investigate the convergence of the high order expansions of the TCL generator. We also discuss circumstances where the exact TCL generator becomes singular for the spin-boson model, and a model of excitation energy transfer in the Fenna-Matthews-Olson complex.     

Stochastic surrogate Hamiltonian

The surrogate Hamiltonian is a general scheme to simulate the many body quantum dynamics composed of a primary system coupled to a bath. The method has been based on a representative bath Hamiltonian composed of two-level systems that is able to mimic the true system-bath dynamics up to a prespecified time. The original surrogate Hamiltonian method is limited to short time dynamics since the size of the Hilbert space required to obtain convergence grows exponentially with time. By randomly swapping bath modes with a secondary thermal reservoir, the method can simulate quantum dynamics of the primary system from short times to thermal equilibrium. By averaging a small number of realizations converged values of the system observables are obtained avoiding the exponential increase in resources. The method is demonstrated for the equilibration of a molecular oscillator with a thermal bath.     

Exact quantum master equation via the calculus on path integrals

An exact quantum master equation formalism is constructed for the efficient evaluation of quantum non-Markovian dissipation beyond the weak system-bath interaction regime in the presence of time-dependent external field. A novel truncation scheme is further proposed and compared with other approaches to close the resulting hierarchically coupled equations of motion. The interplay between system-bath interaction strength, non-Markovian property, and required level of hierarchy is also demonstrated with the aid of simple spin-boson systems.