[论文解读] Adding quantum effects to the semi-classical molecular dynamics simulations

Siyang Yang|arXiv (Cornell University)|Oct 18, 2011

Spectroscopy and Quantum Chemical Studies参考文献 2被引用 1

一句话总结

本文提出了一套用于半经典的分子动力学的新型运动方程（EOM），通过引入全局量子力，将诸如隧穿等量子效应纳入其中。该方法保持了量子力学能量守恒，并使轨迹能够隧穿高于其能量的势垒，从而在非绝热体系中提升了精度，同时增强了量子与经典自由度之间的自洽性。

ABSTRACT

Simulating the molecular dynamics (MD) using classical or semi-classical trajectories provides important details for the understanding of many chemical reactions, protein folding, drug design, and solvation effects. MD simulations using trajectories have achieved great successes in the computer simulations of various systems, but it is difficult to incorporate quantum effects in a robust way. Therefore, improving quantum wavepacket dynamics and incorporating nonadiabatic transitions and quantum effects into classical and semi-classical molecular dynamics is critical as well as challenging. In this paper, we present a MD scheme in which a new set of equations of motion (EOM) are proposed to effectively propagate nuclear trajectories while conserving quantum mechanical energy which is critical for describing quantum effects like tunneling. The new quantum EOM is tested on a one-state one-dimensional and a two-state two-dimensional model nonadiabatic systems. The global quantum force experienced by each trajectory promotes energy redistribution among the bundle of trajectories, and thus helps the individual trajectory tunnel through the potential barrier higher than the energy of the trajectory itself. Construction of the new quantum force and EOM also provides a better way to treat the issue of back-reaction in mixed quantum-classical (MQC) methods, i.e. self-consistency between quantum degrees of freedom (DOF) and classical DOF.

研究动机与目标

• 为在经典和半经典的分子动力学模拟中稳健地引入量子效应（尤其是隧穿）提供解决方案。
• 改进混合量子-经典（MQC）框架中非绝热跃迁和能量重分布的处理方法。
• 通过确保量子与经典自由度之间的自洽性，解决反作用问题。
• 开发一种系统化的方法，用于传播核轨迹，同时保持量子力学能量守恒。
• 在模型体系上测试新的运动方程，以验证其在半经典环境中描述量子现象的能力。

提出的方法

• 提出了一套新的运动方程（EOM），其中引入了源自量子波包动力学的全局量子力。
• 引入一种作用于每个经典轨迹的量子力，通过在轨迹束中重新分配能量，使隧穿能够穿过经典禁阻的势垒。
• 通过在整个传播过程中强制保持量子力学能量期望值，确保能量守恒。
• 构建量子力，使其源自量子态的动力学，并与经典核运动一致耦合。
• 将该方法应用于单态一维和双态二维非绝热模型体系，以测试其性能。
• 利用该框架处理反作用效应，确保量子自由度以自洽方式影响经典轨迹。

实验结果

研究问题

• RQ1 如何能有效将诸如隧穿等量子效应整合到半经典的分子动力学模拟中？
• RQ2 能否构建一种全局量子力，使单个轨迹能够隧穿高于其能量的势垒？
• RQ3 所提出的EOM如何改善混合量子-经典动力学中量子与经典自由度之间的自洽性？
• RQ4 通过量子力实现的能量重分布对模型体系中非绝热动力学精度有何影响？
• RQ5 新EOM能否在传播经典轨迹的同时保持量子力学能量守恒？

主要发现

• 新运动方程成功传播了核轨迹，同时保持了量子力学能量守恒，实现了对量子效应的精确描述。
• 全局量子力使单个轨迹能够隧穿高于其能量的势垒，这是量子隧穿的关键特征。
• 通过量子力在轨迹束中实现的能量重分布，增强了系统探索经典禁阻区域的能力。
• 该方法通过确保量子与经典自由度之间的一致耦合，改进了混合量子-经典动力学中反作用效应的处理。
• 该方法在单维和双维非绝热模型体系上得到验证，证明其在捕捉量子动力学方面的可行性与有效性。
• 该框架提供了一种系统且稳健的方法，可在不依赖完整量子动力学的前提下，将量子效应整合到半经典的MD中。