[论文解读] Canonical sampling through velocity-rescaling

Giovanni Bussi, Davide Donadio+1 more|arXiv (Cornell University)|Mar 28, 2008

Spectroscopy and Quantum Chemical Studies被引用 4

一句话总结

本文提出了一种新颖的分子动力学系综平衡器，通过从正态分布变量导出的随机因子对所有粒子速度进行随机重标度，以强制实现正则系综，确保对目标温度的正确采样。该方法保持了一个守恒的‘有效能量’量，可用于误差监控，并在Lennard-Jones和TIP4P水模型中表现出色，其动态性质在很大程度上与系综平衡器的弛豫时间无关。

ABSTRACT

We present a new molecular dynamics algorithm for sampling the canonical distribution. In this approach the velocities of all the particles are rescaled by a properly chosen random factor. The algorithm is formally justified and it is shown that, in spite of its stochastic nature, a quantity can still be defined that remains constant during the evolution. In numerical applications this quantity can be used to measure the accuracy of the sampling. We illustrate the properties of this new method on Lennard-Jones and TIP4P water models in the solid and liquid phases. Its performance is excellent and largely independent on the thermostat parameter also with regard to the dynamic properties.

研究动机与目标

• 开发一种简单、鲁棒的分子动力学平衡器，能够准确且高效地采样正则系综。
• 解决Berendsen和Langevin等流行随机平衡器缺乏守恒量的问题，从而阻碍误差监控与代码验证。
• 通过引入恰当选择的随机因子，推广速度重标度方法，以确保正确的统计采样，同时保持数值稳定性。
• 引入一种新的守恒量——‘有效能量’，以实现对有限时间步长误差和系综偏离的可靠评估。
• 在真实体系（包括Lennard-Jones和TIP4P水）中，通过固相和液相的测试，证明该方法的有效性与遍历性。

提出的方法

• 该算法通过从一个分布中抽取的随机因子α对所有粒子速度进行重标度，以确保动能与正则系综的目标值匹配。
• 重标度因子α源自一种类似于Ornstein-Uhlenbeck过程的随机过程，其弛豫时间τ控制随机性的强度。
• 通过证明一种新的守恒量——‘有效能量’——在动力学过程中保持不变，从而形式化地证明了该方法的合理性，类似于NVE模拟中的能量守恒。
• 有效能量被定义为动能与随机变量的函数，从而可直接监控系综偏离与时间步长误差。
• 该方法通过增加一个确保正确采样的随机分量，推广了Berendsen的平衡器，同时保留了其稳定性和简洁性。
• 该算法通过涉及高斯随机变量和伽马分布和的闭式表达式实现α²，从而实现高效的数值计算。

实验结果

研究问题

• RQ1 能否通过引入随机重标度，将速度重标度方法推广以正确采样正则系综？
• RQ2 引入随机重标度因子是否会产生一个守恒量，从而实现在分子动力学模拟中的误差监控？
• RQ3 与Nosé-Hoover和Berendsen等成熟方法相比，该新平衡器在遍历性与动态性质方面的性能如何？
• RQ4 系统的动态行为（如自扩散系数）在多大程度上受平衡器弛豫时间τ的影响？
• RQ5 该新引入的‘有效能量’量能否可靠检测有限时间步长误差，并确保NVT系综的正确采样？

主要发现

• 在遍历性假设下，该方法被分析证明可正确采样正则系综。
• 识别出一种新的守恒量——‘有效能量’，并证明其在模拟过程中保持恒定，从而实现稳健的误差监控与时间步长选择。
• 该方法在Lennard-Jones和TIP4P水模型的固相与液相中均表现出色，自扩散系数D等动态性质在很大程度上与平衡器弛豫时间τ无关。
• 在300 K的TIP4P水中，该新平衡器计算的自扩散系数D与文献值以及NVE模拟结果一致，仅在τ = 0.002 ps时出现微小偏差。
• 在数值测试中，该方法比Nosé-Hoover平衡器更具遍历性，表明其对相空间的采样更优。
• 该算法实现简单，并提供了一个统一框架，可将Berendsen、Nosé和当前提出的平衡器作为特例统一包含在内。