[论文解读] Continuum limit and preconditioned Langevin sampling of the path integral molecular dynamics

Jianfeng Lu, Yulong Lu+1 more|arXiv (Cornell University)|Nov 27, 2018

Quantum, superfluid, helium dynamics参考文献 44被引用 5

一句话总结

本文提出预处理朗之万动力学以实现路径积分分子动力学（PIMD）中的高效采样，利用环聚合物在 N→∞ 时的连续极限。它在二次势场中建立了与维度无关的指数收敛至平衡态，并提供了一个与非绝热多能级量子系统兼容的框架。

ABSTRACT

We investigate the continuum limit that the number of beads goes to infinity in the ring polymer representation of thermal averages. Studying the continuum limit of the trajectory sampling equation sheds light on possible preconditioning techniques for sampling ring polymer configurations with large number of beads. We propose two preconditioned Langevin sampling dynamics, which are shown to have improved stability and sampling accuracy. We present a careful mode analysis of the preconditioned dynamics and show their connections to the normal mode, the staging coordinate and the Matsubara mode representation for ring polymers. In the case where the potential is quadratic, we show that the continuum limit of the preconditioned mass modified Langevin dynamics converges to its equilibrium exponentially fast, which suggests that the finite-dimensional counterpart has a dimension-independent convergence rate. In addition, the preconditioning techniques can be naturally applied to the multi-level quantum systems in the nonadiabatic regime, which are compatible with various numerical approaches.

研究动机与目标

• 解决由于珠子数 N 增加时刚性弹簧势导致的路径积分分子动力学（PIMD）中采样效率低下的问题。
• 研究环聚合物动力学在 N→∞ 时的连续极限，以指导高珠子系统中的预处理策略。
• 开发预处理朗之万动力学，以提高PIMD模拟中的数值稳定性和采样精度。
• 在二次势场情况下，为预处理动力学建立理论收敛保证。
• 通过绝热基下的环聚合物表示，将该框架扩展至非绝热多能级量子系统。

提出的方法

• 将环聚合物的连续极限形式化为配置空间中的布朗回路，能量定义在环路空间 Lℝ^d 上。
• 通过 E(𝒬) = ∫₀^β [½|∂τ𝒬(τ)|² + V(𝒬(τ))] dτ 在环路空间上定义一个无穷维吉布斯测度 π。
• 提出两种预处理朗之万动力学：一种基于质量修改，另一种使用预处理矩阵 L^α。
• 利用模态分解分析动力学，将其与正则模态、分阶段坐标和马苏巴里模态联系起来。
• 推导预处理动力学的福克-普朗克方程，并验证不变测度与目标吉布斯测度一致。
• 通过绝热基下的环聚合物表示，将形式化扩展至两能级量子系统，给出有效哈密顿量和可观测量的权重函数。

实验结果

研究问题

• RQ1 当珠子数 N→∞ 时，环聚合物轨迹的连续极限行为如何？
• RQ2 在高珠子PIMD模拟中，哪些预处理策略可实现稳定与加速采样？
• RQ3 在二次势场情况下，预处理朗之万动力学能否实现与维度无关的收敛速率？
• RQ4 所提出的动力学与正则模态、分阶段坐标和马苏巴里模态等已有表示之间有何关联？
• RQ5 该预处理框架能否扩展至非绝热多能级量子系统？

主要发现

• 在二次势场情况下，预处理质量修改朗之万动力学以与珠子数无关的速率指数收敛至平衡态。
• 轨迹采样方程的连续极限揭示了结构洞见，可指导有限-N PIMD 模拟中的有效预处理。
• 所提出的预处理动力学被证明与多种数值方法兼容，并可自然推广至非绝热系统。
• 模态分析表明，预处理动力学与马苏巴里模态表示之间存在联系，提示其具有改进的谱性质。
• 对于两能级量子系统，使用绝热基的环聚合物表示可借助所提动力学实现热平均的一致采样。
• 预处理动力学的不变测度与目标吉布斯测度一致，证实了采样方案的正确性。