[论文解读] An Electronically Non-Adiabatic Generalization of Ring Polymer Molecular Dynamics

Timothy J. H. Hele|arXiv (Cornell University)|Aug 19, 2013

Advanced Physical and Chemical Molecular Interactions参考文献 7被引用 17

一句话总结

本文提出了一种非绝热环聚合物分子动力学（RPMD）速率理论的推广，能够准确模拟电子非绝热体系中的量子效应和穿垒动力学。该方法通过在修正的势能面上引入电子态耦合，将经典RPMD扩展至非绝热体系，即使在深度隧穿区域，其计算结果与精确量子速率的偏差也在两倍以内，且计算成本随体系规模线性增长。

ABSTRACT

In this thesis I generalize Ring Polymer Molecular Dynamics (RPMD) rate theory to electronically non-adiabatic systems, followed by application to two one-dimensional curve crossing models and a multidimensional spin-boson model.

研究动机与目标

• 将环聚合物分子动力学（RPMD）速率理论推广至核运动发生在多个电子势能面上的电子非绝热体系。
• 开发一种计算高效的计算方法，以捕捉非绝热过程中诸如隧穿和反复穿越等量子效应。
• 确保该方法在单电子态极限下退化为传统的绝热RPMD，并在高温极限下保持物理一致性。
• 在一维曲线交叉模型和具有精确量子基准的多维自旋-玻色子模型上验证该方法的有效性。
• 证明该方法的计算成本随体系规模线性增长，同时在准确性上优于其他近似方法。

提出的方法

• 通过将量子体系表示为通过谐弹簧连接的类经典粒子环链，推广经典同构理论，扩展相空间以包含电子自由度。
• 推导出一种新的环聚合物有效势能面，通过哈密顿量中修改的耦合项引入电子非绝热耦合。
• 应用Bennett-Chandler分解方法，将RPMD速率分解为在分割面上的质心概率密度与产物区域中正动量质心通量的乘积。
• 在非对称曲线交叉模型中采用变分法确定最优分割面。
• 在扩展相空间上使用经典分子动力学实现该方法，反应速率通过环聚合物质心的时间相关函数计算得出。
• 通过将发散的量子过渡态理论（QTST）速率分解为两个物理解释明确且有限的贡献，解决了自旋-玻色子模型中QTST速率发散的挑战。

实验结果

研究问题

• RQ1 是否可以将RPMD速率理论推广至电子非绝热体系，同时保持其准确性和计算效率？
• RQ2 在对称与非对称耦合的一维曲线交叉模型中，非绝热RPMD速率与精确量子速率相比如何？
• RQ3 该方法能否准确描述具有德拜谱密度的自旋-玻色子模型，特别是在反应坐标质量趋近于零的极限下？
• RQ4 当电子态解耦时，非绝热RPMD方法是否能正确恢复绝热RPMD极限？
• RQ5 在标准QTST因质量因子发散或速率未定义而失效的体系中，RPMD速率是否可被可靠计算？

主要发现

• 在一维对称曲线交叉模型中，即使在100 K时，经典速率失效数量级，RPMD速率仍与精确量子速率保持在两倍以内。
• 在非对称曲线交叉模型中，RPMD对精确量子速率的近似效果甚至优于对称情况，偏差程度与绝热极限下的偏差相当。
• 对于具有德拜谱密度的自旋-玻色子模型，尽管由于质量因子发散导致QTST未定义，但通过收敛分解，RPMD速率仍可计算。
• 非绝热RPMD速率与精确量子速率的偏差小于几乎所有其他测试的近似方法，包括流行的均场和半经典方法。
• 当电子耦合关闭时，该方法正确退化为传统的绝热RPMD，验证了与已有理论的一致性。
• 计算成本随体系规模线性增长，使其适用于大规模电子非绝热动力学模拟。