Nonadiabatic quantum-classical reaction rates with quantum equilibrium structure

Time correlation function expressions for quantum reaction-rate coefficients are computed in a quantum-classical limit. This form for the correlation function retains the full quantum equilibrium structure of the system in the spectral density function but approximates the time evolution of the operator by quantum-classical Liouville dynamics. Approximate analytical expressions for the spectral density function, which incorporate quantum effects in the many-body environment and reaction coordinate, are derived. The results of numerical simulations of the reaction rate are presented for a reaction model in which a two-level system is coupled to a bistable oscillator which is, in turn, coupled to a bath of harmonic oscillators. The nonadiabatic quantum-classical dynamics is simulated in terms of an ensemble of surface-hopping trajectories and the effects of the quantum equilibrium structure on the reaction rate are discussed.     

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.     

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.     

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.     

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.     

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.     

Correlated driving-and-dissipation equation for non-Condon spectroscopy with the Herzberg–Teller vibronic coupling

Correlated driving-and-dissipation equation (CODDE) is an optimized complete second-order quantum dissipation approach, which is originally concerned with the reduced system dynamics only. However, one can actually extract the hybridized bath dynamics from CODDE with the aid of dissipaton-equation-of-motion theory, a statistical quasi-particle quantum dissipation formalism. Treated as a one–dissipaton theory, CODDE is successfully extended to deal with the Herzberg–Teller vibronic couplings in dipole–field interactions. Demonstrations will be carried out on the non-Condon spectroscopies of a model dimer system.     

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.     

Hybrid quantum/classical studies of photodissociation and recombination of I2 (A) in rare gas matrices: A linear chain model

A hybrid quantum/classical model is developed for the photodissociation and recombination dynamics of an I₂ molecule in low-temperature rare-gas (Rg) matrices. The simplified model consists of an I₂ molecule embedded in a linear chain of Rg atoms. The aggregate is partitioned into a quantum system and a classical bath, which are self-consistently coupled. Two partitioning schemes are used. The first treats the I-I coordinate quantum mechanically and the Rg coordinates classically. The second and more reliable scheme includes in the quantum system both the I-I mode and the symmetric motion of the two nearest Rg atoms. Both models show substantial energy transfer from the dissociating I₂ to the solvent, followed by coherent vibrational motion of the recombined I₂. It is found that the one-dimensional quantum/classical scheme is consistent with its higher dimensional counterpart, although the latter shows much faster dephasing. © 1996 John Wiley & Sons, Inc.     

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.     

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.     

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.     

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.