Grunwald Robbie, Kapral Raymond
Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada.
J Chem Phys. 2007 Mar 21;126(11):114109. doi: 10.1063/1.2567164.
The conditions under which quantum-classical Liouville dynamics may be reduced to a master equation are investigated. Systems that can be partitioned into a quantum-classical subsystem interacting with a classical bath are considered. Starting with an exact non-Markovian equation for the diagonal elements of the density matrix, an evolution equation for the subsystem density matrix is derived. One contribution to this equation contains the bath average of a memory kernel that accounts for all coherences in the system. It is shown to be a rapidly decaying function, motivating a Markovian approximation on this term in the evolution equation. The resulting subsystem density matrix equation is still non-Markovian due to the fact that bath degrees of freedom have been projected out of the dynamics. Provided the computation of nonequilibrium average values or correlation functions is considered, the non-Markovian character of this equation can be removed by lifting the equation into the full phase space of the system. This leads to a trajectory description of the dynamics where each fictitious trajectory accounts for decoherence due to the bath degrees of freedom. The results are illustrated by computations of the rate constant of a model nonadiabatic chemical reaction.
研究了量子 - 经典刘维尔动力学可简化为主方程的条件。考虑了可划分为与经典热库相互作用的量子 - 经典子系统的体系。从密度矩阵对角元的精确非马尔可夫方程出发,推导了子系统密度矩阵的演化方程。该方程的一项包含一个记忆核的热库平均值,该记忆核考虑了系统中的所有相干性。它被证明是一个快速衰减的函数,这促使在演化方程中对该项进行马尔可夫近似。由于热库自由度已从动力学中投影出去,所得的子系统密度矩阵方程仍然是非马尔可夫的。如果考虑非平衡平均值或关联函数的计算,通过将该方程提升到系统的全相空间,可以消除该方程的非马尔可夫特性。这导致了动力学的轨迹描述,其中每个虚拟轨迹考虑了由于热库自由度引起的退相干。通过计算一个模型非绝热化学反应的速率常数来说明这些结果。