[论文解读] Frozen Gaussian approximation with surface hopping for mixed quantum-classical dynamics: A mathematical justification of fewest switches surface hopping algorithms

Jianfeng Lu, Zhennan Zhou|arXiv (Cornell University)|Feb 20, 2016

Spectroscopy and Quantum Chemical Studies参考文献 28被引用 6

一句话总结

本文提出了在半经典 regime下混合量子-经典动力学的表面跳跃算法——冻结高斯近似与表面跳跃（FGA-SH）的数学上严格的推导。通过将解表示为路径积分的随机形式，该方法提供了多组绝热态下薛定谔方程的蒙特卡洛近似，在半经典极限下实现渐近正确性，其一致性误差与稳定性误差分别以 $\mathcal{O}(\varepsilon)$ 和 $\mathcal{O}(\varepsilon^2)$ 的速率衰减。

ABSTRACT

We develop a surface hopping algorithm based on frozen Gaussian approximation for semiclassical matrix Schrödinger equations, in the spirit of Tully's fewest switches surface hopping method. The algorithm is asymptotically derived from the Schrödinger equation with rigorous approximation error analysis. The resulting algorithm can be viewed as a path integral stochastic representation of the semiclassical matrix Schrödinger equations. Our results provide mathematical understanding to and shed new light on the important class of surface hopping methods in theoretical and computational chemistry.

研究动机与目标

• 为非绝热量子动力学中广泛使用的最少开关表面跳跃（FSSH）算法提供严格的数学基础。
• 解决尽管表面跳跃方法在计算化学中广泛应用，却长期缺乏理论依据的问题。
• 基于冻结高斯近似与路径积分随机表示，推导一种新算法 FGA-SH，确保在半经典极限下具有渐近正确性。
• 建立 FGA-SH 解与矩阵薛定谔方程真实解之间的定量误差估计。
• 通过将现有 FSSH 方法与解的明确定义的随机表示相联系，澄清其中存在的概念模糊性。

提出的方法

• 推导半经典矩阵薛定谔方程解的路径积分随机表示，将动力学视为相空间中的分段确定性轨迹。
• 将 FGA-SH 算法构建为一种用于计算路径积分的蒙特卡洛方法，其中轨迹在不同势能面上按经典哈密顿流演化。
• 在非绝热耦合点实施表面跳跃事件，跳跃概率由哈密顿量中的导数耦合项导出。
• 采用冻结高斯近似（Herman-Kluk传播器）作为波函数的试探解，通过相空间积分实现高效计算。
• 利用一种新颖的稳定性引理，建立一致性与稳定性估计，该引理将解的误差以近似薛定谔方程余项的形式进行有界控制。
• 证明 FGA-SH 近似中的总误差在一致性上受 $\mathcal{O}(\varepsilon)$ 限制，在稳定性上受 $\mathcal{O}(\varepsilon^2)$ 限制，常数依赖于时间和系统参数。

实验结果

研究问题

• RQ1 最少开关表面跳跃算法能否在半经典 regime 下从薛定谔方程严格推导而出？
• RQ2 包含非绝热跃迁的混合量子-经典动力学的路径积分表示背后的数学结构是什么？
• RQ3 如何将冻结高斯近似扩展至包含多态系统中表面跳跃的场景，并实现精确的误差控制？
• RQ4 FGA-SH 算法与量子-经典刘维尔方程之间存在何种关系？能否从后者对前者进行理论验证？
• RQ5 FGA-SH 近似在半经典参数 $\varepsilon$ 下的精确误差界是什么？

主要发现

• FGA-SH 算法被推导为半经典矩阵薛定谔方程解的路径积分随机表示，为表面跳跃提供了数学上严谨的基础。
• FGA-SH 近似的一致性误差受 $C\varepsilon \frac{(C_t)^n}{n!} + C\varepsilon^2 e^{C_t}$ 控制，当 $n \to \infty$ 时趋于零，确保了渐近正确性。
• 稳定性引理表明，近似解与真实解之间的误差最多以 $\frac{1}{\varepsilon} \int_0^t \mu^\varepsilon(s) ds$ 的速率增长，其中 $\mu^\varepsilon(s)$ 控制近似方程中的余项。
• FGA-SH 近似中的总误差在一致性上受 $\mathcal{O}(\varepsilon)$ 限制，在稳定性上受 $\mathcal{O}(\varepsilon^2)$ 限制，常数依赖于时间和系统参数。
• 该方法为最少开关表面跳跃算法提供了严格的理论依据，解决了其在非绝热动力学中长期存在的理论基础疑问。
• 该方法揭示，表面跳跃可被解释为对高维路径积分的直接蒙特卡洛求值，为理解该算法提供了全新的概念框架。