Pomyalov Anna, Tannor David J
Department of Chemical Physics, The Weizmann Institute of Science, Rehovot, 76100 Israel.
J Chem Phys. 2005 Nov 22;123(20):204111. doi: 10.1063/1.2121649.
The calculation of chemical reaction rates in the condensed phase is a central preoccupation of theoretical chemistry. At low temperatures, quantum-mechanical effects can be significant and even dominant; yet quantum calculations of rate constants are extremely challenging, requiring theories and methods capable of describing quantum evolution in the presence of dissipation. In this paper we present a new approach based on the use of a non-Markovian quantum master equation (NM-QME). As opposed to other approximate quantum methods, the quantum dynamics of the system coordinate is treated exactly; hence there is no loss of accuracy at low temperatures. However, because of the perturbative nature of the NM-QME it breaks down for dimensionless frictions larger than about 0.1. We show that by augmenting the system coordinate with a collective mode of the bath, the regime of validity of the non-Markovian master equation can be extended significantly, up to dimensionless frictions of 0.5 over the entire temperature range. In the energy representation, the scaling goes as the number of levels in the relevant energy range to the third power. This scaling is not prohibitive even for chemical systems with many levels; hence we believe that the current method will find a useful place alongside the existing techniques for calculating quantum condensed-phase rate constants.
凝聚相中化学反应速率的计算是理论化学的核心关注点。在低温下,量子力学效应可能很显著甚至占主导;然而,速率常数的量子计算极具挑战性,需要能够描述存在耗散时量子演化的理论和方法。在本文中,我们提出一种基于使用非马尔可夫量子主方程(NM-QME)的新方法。与其他近似量子方法不同,系统坐标的量子动力学得到精确处理;因此在低温下不会损失精度。然而,由于NM-QME的微扰性质,对于大于约0.1的无量纲摩擦它会失效。我们表明,通过用浴的集体模式扩充系统坐标,非马尔可夫主方程的有效范围可以显著扩展,在整个温度范围内可达无量纲摩擦为0.5。在能量表象中,标度与相关能量范围内的能级数量的三次方成正比。即使对于具有许多能级的化学系统,这种标度也并非过高;因此我们相信,当前方法将在计算量子凝聚相速率常数的现有技术中找到有用的位置。