[论文解读] Ab initio instanton rate theory made efficient using Gaussian process regression

Gabriel Laude, Danilo Calderini+2 more|arXiv (Cornell University)|May 7, 2018

Advanced Chemical Physics Studies被引用 4

一句话总结

本文提出了一种高斯过程回归（GPR）加速的从头算瞬子方法，显著减少了准确预测量子反应速率所需的高水平电子结构计算次数。通过使用GPR在瞬子路径附近局部拟合势能面，该方法在仅需标准方法约1/10次从头算计算的情况下，即可将结果收敛至基准值的1%以内，从而实现了对H + C2H6等反应在UCCSD(T)-F12b水平下的高效、高精度速率计算。

ABSTRACT

Ab initio instanton rate theory is a computational method for rigorously including tunnelling effects into calculations of chemical reaction rates based on a potential-energy surface computed on the fly from electronic-structure theory. This approach is necessary to extend conventional transition-state theory into the deep-tunnelling regime, but is also more computationally expensive as it requires many more ab initio calculations. We propose an approach which uses Gaussian process regression to fit the potential-energy surface locally around the dominant tunnelling pathway. The method can be converged to give the same result as from an on-the-fly ab initio instanton calculation but requires far fewer electronic-structure calculations. This makes it a practical approach for obtaining accurate rate constants based on high-level electronic-structure methods. We show fast convergence to reproduce benchmark H + CH4 results and evaluate new low-temperature rates of H + C2H6 in full dimensionality at a UCCSD(T)-F12b/cc-pVTZ-F12 level.

研究动机与目标

• 降低从头算瞬子理论的计算成本，该理论因需要沿完整瞬子路径获取能量、梯度和Hessian矩阵而需进行大量高水平电子结构计算。
• 克服即时计算（on-the-fly）从头算瞬子计算的低效性，其计算成本高于过渡态理论，原因在于所需电子结构计算次数过多。
• 开发一种在显著减少昂贵计算次数的同时，仍能保持高水平电子结构方法高精度的方法。
• 通过使高阶理论水平下的计算变得可行，实现从头算瞬子理论在复杂多原子反应中的实际应用。

提出的方法

• 使用高斯过程回归（GPR）构建代理势能面（PES），仅基于沿瞬子路径采样的从头算数据点进行训练。
• GPR模型基于沿环聚合物瞬子路径的关键点上计算的少量、自适应的高水平电子结构计算数据（能量、梯度、Hessian矩阵）进行训练。
• 该方法根据预测不确定性自适应选择新的训练点，以确保收敛至真实瞬子路径和速率常数。
• GPR代理模型使得在远少于标准即时计算瞬子方法所需的高水平电子结构计算次数下，实现精确的瞬子速率计算。
• 通过将收敛速率与基准全从头算瞬子计算结果进行比较，验证了该方法，结果显示结果在1%以内一致。
• 该方法被应用于H + CH4和H + C2H6反应，在仅约6个Hessian矩阵计算下即实现了UCCSD(T)-F12b/cc-pVTZ-F12水平下的全维瞬子速率计算。

实验结果

研究问题

• RQ1 能否使用高斯过程回归在仅使用极少高水平电子结构计算的情况下，准确近似瞬子路径上的势能面？
• RQ2 为使瞬子速率收敛至全即时计算结果的1%以内，所需的最少从头算计算次数是多少？
• RQ3 与传统从头算瞬子理论相比，GPR加速的瞬子方法在H + CH4等基准反应中的精度和效率如何？
• RQ4 该方法能否扩展至H + C2H6等复杂多原子反应，在高阶理论水平下实现计算成本可控的实际应用？
• RQ5 能否将GPR框架与低水平电子结构数据结合，进一步减少高水平计算次数而不损失精度？

主要发现

• 对于H + CH4反应，GPR加速的瞬子方法仅使用标准方法所需电子结构计算次数的一小部分，即可将结果收敛至基准全从头算瞬子速率的1%以内。
• 对于H + C2H6反应，该方法在仅约6个Hessian矩阵计算下，即实现了UCCSD(T)-F12b/cc-pVTZ-F12水平下的全维瞬子速率计算，其效率接近经典过渡态理论（TST）。
• 即使瞬子路径与最小能量路径显著偏离，该方法仍保持高精度，表明其在复杂反应机理中的鲁棒性。
• 使用GPR可实现对环聚合物粒子数和训练集大小的系统性收敛，确保结果的可靠性。
• 该方法使得在瞬子理论中应用高水平电子结构方法成为可能，而此前这些方法对许多体系而言计算成本过高。
• 该方法为在复杂反应网络中系统性地引入隧道效应打开了大门，仅需对关键步骤进行高效、高精度的瞬子计算。