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.     

Transport properties of quantum-classical systems

Correlation function expressions for calculating transport coefficients for quantum-classical systems are derived. The results are obtained by starting with quantum transport coefficient expressions and replacing the quantum time evolution with quantum-classical Liouville evolution, while retaining the full quantum equilibrium structure through the spectral density function. The method provides a variety of routes for simulating transport coefficients of mixed quantum-classical systems, composed of a quantum subsystem and a classical bath, by selecting different but equivalent time evolution schemes of any operator or the spectral density. The structure of the spectral density is examined for a single harmonic oscillator where exact analytical results can be obtained. The utility of the formulation is illustrated by considering the rate constant of an activated quantum transfer process that can be described by a many-body bath reaction coordinate.     

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.     

Fully adaptive propagation of the quantum-classical Liouville equation

In mixed quantum-classical molecular dynamics few but important degrees of freedom of a dynamical system are modeled quantum-mechanically while the remaining ones are treated within the classical approximation. Rothe methods established in the theory of partial differential equations are used to control both temporal and spatial discretization errors on grounds of a global tolerance criterion. The TRAIL (trapezoidal rule for adaptive integration of Liouville dynamics) scheme [I. Horenko and M. Weiser, J. Comput. Chem. 24, 1921 (2003)] has been extended to account for nonadiabatic effects in molecular dynamics described by the quantum-classical Liouville equation. In the context of particle methods, the quality of the spatial approximation of the phase-space distributions is maximized while the numerical condition of the least-squares problem for the parameters of particles is minimized. The resulting dynamical scheme is based on a simultaneous propagation of moving particles (Gaussian and Dirac deltalike trajectories) in phase space employing a fully adaptive strategy to upgrade Dirac to Gaussian particles and, vice versa, downgrading Gaussians to Dirac-type trajectories. This allows for the combination of Monte-Carlo-based strategies for the sampling of densities and coherences in multidimensional problems with deterministic treatment of nonadiabatic effects. Numerical examples demonstrate the application of the method to spin-boson systems in different dimensionality. Nonadiabatic effects occurring at conical intersections are treated in the diabatic representation. By decreasing the global tolerance, the numerical solution obtained from the TRAIL scheme are shown to converge towards exact results.     

Mixed quantum-classical equilibrium: Surface hopping

We re-examine the analysis of the equilibrium limits of the fewest switches surface hopping algorithm for mixed quantum-classical dynamics. In contrast with previously reported results, we show that surface hopping does not, in general, exactly yield Boltzmann equilibrium, but that in practice the observed deviations are quite small. We also demonstrate that surface hopping does approach the exact equilibrium distribution in both the limits of small adiabatic splitting and/or strong nonadiabatic coupling. We verify these analytical results with numerical simulations for a simple two-level quantum system connected to a bath of classical particles.     

Quantum-classical Liouville dynamics of proton and deuteron transfer rates in a solvated hydrogen-bonded complex

Proton and deuteron transfer reactions in a hydrogen-bonded complex dissolved in a polar solution are studied using quantum-classical Liouville dynamics. Reactive-flux correlation functions that involve quantum-classical Liouville dynamics for species operators and quantum equilibrium sampling are used to calculate the rate constants. Adiabatic and nonadiabatic reaction rates are computed, compared, and analyzed. Large variations of the kinetic isotope effect (KIE) for this reaction have been observed in the literature, which depend on the nature of the approximate calculation used to estimate the proton and deuteron transfer rates. Our estimate of the KIE lies at the low end of the range of previously observed values, suggesting a rather small KIE for this reaction.     

多原子体系非绝热动力学的近似理论方法

非绝热动力学普遍存在于光物理和光化学过程中。描述非绝热跃迁需要处理电子-原子核间的相互耦合运动。由于计算量随体系尺度增大剧烈增长,准确的量子动力学计算目前只适用于描述小分子体系。为了研究多原子分子体系的非绝热过程,近年来发展了一些基于量子-经典动力学近似方法。本文将对典型的这类方法包括经典Ehrenfest方法、面跳跃方法、基于Wigner表示的混合量子-经典方法进行简要的介绍,并讨论如何将量子-经典动力学方法与电子结构从头算手段相结合,模拟非绝热过程。重点阐明各种方法的基本思想和优缺点,并对该领域的发展进行展望。     

Quantum-classical Liouville dynamics of nonadiabatic proton transfer

A proton transfer reaction in a linear hydrogen-bonded complex dissolved in a polar solvent is studied using mixed quantum-classical Liouville dynamics. In this system, the proton is treated quantum mechanically and the remainder of the degrees of freedom is treated classically. The rates and mechanisms of the reaction are investigated using both adiabatic and nonadiabatic molecular dynamics. We use a nonadiabatic dynamics algorithm which allows the system to evolve on single adiabatic surfaces and on coherently coupled pairs of adiabatic surfaces. Reactive-flux correlation function expressions are used to compute the rate coefficients and the role of the dynamics on the coherently coupled surfaces is elucidated.     

Army ants algorithm for rare event sampling of delocalized nonadiabatic transitions by trajectory surface hopping and the estimation of sampling errors by the bootstrap method

The most widely used algorithm for Monte Carlo sampling of electronic transitions in trajectory surface hopping (TSH) calculations is the so-called anteater algorithm, which is inefficient for sampling low-probability nonadiabatic events. We present a new sampling scheme (called the army ants algorithm) for carrying out TSH calculations that is applicable to systems with any strength of coupling. The army ants algorithm is a form of rare event sampling whose efficiency is controlled by an input parameter. By choosing a suitable value of the input parameter the army ants algorithm can be reduced to the anteater algorithm (which is efficient for strongly coupled cases), and by optimizing the parameter the army ants algorithm may be efficiently applied to systems with low-probability events. To demonstrate the efficiency of the army ants algorithm, we performed atom-diatom scattering calculations on a model system involving weakly coupled electronic states. Fully converged quantum mechanical calculations were performed, and the probabilities for nonadiabatic reaction and nonreactive deexcitation (quenching) were found to be on the order of 10(-8). For such low-probability events the anteater sampling scheme requires a large number of trajectories ( approximately 10(10)) to obtain good statistics and converged semiclassical results. In contrast by using the new army ants algorithm converged results were obtained by running 10(5) trajectories. Furthermore, the results were found to be in excellent agreement with the quantum mechanical results. Sampling errors were estimated using the bootstrap method, which is validated for use with the army ants algorithm.     

Mixed quantum-classical description of spectroscopy of dissipative systems

Mixed quantum-classical statistical mechanics is employed to calculate dipole moment correlation function and linear absorption spectra. A quantum two-level subsystem interacting with quantum vibrations (primary oscillators) which in turn are coupled to a classical bath composed of infinite set of harmonic oscillators is used as a dissipative system. Starting with mixed quantum-classical Liouville equation for the evaluation of the mixed quantum-classical dipole moment correlation function and using coherent states and the inverse of Baker-Campbell-Hausdorf formula to evaluate the trace over the primary oscillators, whereby, a closed analytical expression for the electronic dipole moment correlation function is obtained. Illustrations of several absorption spectra at different temperatures are provided. An approximate optical four-point correlation is obtained in the high temperature limit. A strategy for deriving an exact optical four-point correlation is suggested.     

混合量子-经典方法计算电荷转移速率及其在实际体系中的应用(英文)

混合量子-经典方法在复杂分子体系动力学过程的模拟方面有重要应用.我们采用Ehrenfest方法、surfacehopping方法和混合量子经典Liouville方程计算了在非绝热极限下的电荷转移速率.然后将这三种方法应用于有机半导体材料电荷转移速率的计算.研究结果发现,Ehrenfest方法和surface hopping方法可能严重偏离正确的结果.偏离的原因是这两种方法没有正确处理相干项的运动,而且这种偏离在涉及到高频模式时显得更加严重.     

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.     

On the properties of a primitive semiclassical surface hopping propagator for nonadiabatic quantum dynamics

A previously developed nonadiabatic semiclassical surface hopping propagator [M. F. Herman J. Chem. Phys. 103, 8081 (1995)] is further studied. The propagator has been shown to satisfy the time-dependent Schrodinger equation (TDSE) through order h, and the O(h2) terms are treated as small errors, consistent with standard semiclassical analysis. Energy is conserved at each hopping point and the change in momentum accompanying each hop is parallel to the direction of the nonadiabatic coupling vector resulting in both transmission and reflection types of hops. Quantum mechanical analysis and numerical calculations presented in this paper show that the h2 terms involving the interstate coupling functions have significant effects on the quantum transition probabilities. Motivated by these data, the h2 terms are analyzed for the nonadiabatic semiclassical propagator. It is shown that the propagator can satisfy the TDSE for multidimensional systems by including another type of nonclassical trajectories that reflect on the same surfaces. This h2 analysis gives three conditions for these three types of trajectories so that their coefficients are uniquely determined. Besides the nonadiabatic semiclassical propagator, a numerically useful quantum propagator in the adiabatic representation is developed to describe nonadiabatic transitions.     

WavePacket: A Matlab package for numerical quantum dynamics. III. Quantum-classical simulations and surface hopping trajectories

WavePacket is an open-source program package for numerical simulations in quantum dynamics. Building on the previous Part I (Schmidt and Lorenz, Comput. Phys. Commun. 2017, 213, 223] and Part II (Schmidt and Hartmann, Comput. Phys. Commun. 2018, 228, 229] which dealt with quantum dynamics of closed and open systems, respectively, the present Part III adds fully classical and mixed quantum-classical propagation techniques to WavePacket. There classical phase-space densities are sampled by trajectories which follow (diabatic or adiabatic) potential energy surfaces. In the vicinity of (genuine or avoided) intersections of those surfaces, trajectories may switch between them. To model these transitions, two classes of stochastic algorithms have been implemented: (1) Tully's fewest switches surface hopping and (2) Landau–Zener-based single switch surface hopping. The latter one offers the advantage of being based on adiabatic energy gaps only, thus not requiring nonadiabatic coupling information any more. The present work describes the MATLAB version of WavePacket 6.1.0, which is essentially an object-oriented rewrite of previous versions, allowing to perform fully classical, quantum-classical and quantum-mechanical simulations on an equal footing, that is, for the same physical system described by the same WavePacket input. The software package is hosted and further developed at the Sourceforge platform, where also extensive Wiki-documentation as well as numerous worked-out demonstration examples with animated graphics are available. © 2019 Wiley Periodicals, Inc.