Junginger Andrej, Craven Galen T, Bartsch Thomas, Revuelta F, Borondo F, Benito R M, Hernandez Rigoberto
Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, GA 30332-0400, USA.
Department of Mathematical Sciences, Loughborough University, Loughborough LE11 3TU, UK.
Phys Chem Chem Phys. 2016 Nov 9;18(44):30270-30281. doi: 10.1039/c6cp02519f.
Reaction rates across time-dependent barriers are difficult to define and difficult to obtain using standard transition state theory approaches because of the complexity of the geometry of the dividing surface separating reactants and products. Using perturbation theory (PT) or Lagrangian descriptors (LDs), we can obtain the transition state trajectory and the associated recrossing-free dividing surface. With the latter, we are able to determine the exact reactant population decay and the corresponding rates to benchmark the PT and LD approaches. Specifically, accurate rates are obtained from a local description regarding only direct barrier crossings and to those obtained from a stability analysis of the transition state trajectory. We find that these benchmarks agree with the PT and LD approaches for obtaining recrossing-free dividing surfaces. This result holds not only for the local dynamics in the vicinity of the barrier top, but also for the global dynamics of particles that are quenched at the reactant or product wells after their sojourn over the barrier region. The double-well structure of the potential allows for long-time dynamics related to collisions with the outside walls that lead to long-time returns in the low-friction regime. This additional global dynamics introduces slow-decay pathways that do not result from the local transition across the recrossing-free dividing surface associated with the transition state trajectory, but can be addressed if that structure is augmented by the population transfer of the long-time returns.
由于分隔反应物和产物的分隔面几何结构复杂,利用标准过渡态理论方法难以定义和获取跨越时间依赖势垒的反应速率。使用微扰理论(PT)或拉格朗日描述符(LD),我们可以获得过渡态轨迹以及相关的无再穿越分隔面。利用后者,我们能够确定准确的反应物布居数衰减以及相应的速率,以此来衡量PT和LD方法。具体而言,准确的速率可从仅关于直接穿越势垒的局部描述中获得,并与从过渡态轨迹稳定性分析中获得的速率进行比较。我们发现这些基准与用于获取无再穿越分隔面的PT和LD方法一致。这一结果不仅适用于势垒顶部附近的局部动力学,也适用于粒子在穿越势垒区域后在反应物或产物阱中猝灭的全局动力学。势的双阱结构允许与外壁碰撞相关的长时间动力学,这在低摩擦 regime 中会导致长时间返回。这种额外的全局动力学引入了缓慢衰减路径,这些路径并非由与过渡态轨迹相关的无再穿越分隔面的局部跃迁产生,但如果通过长时间返回的布居转移来增强该结构,就可以解决这个问题。