[论文解读] Canonical Temperature Control by Molecular Dynamics

William G. Hoover, Carol G. Hoover|arXiv (Cornell University)|Mar 3, 2024

Protein Structure and Dynamics被引用 1

一句话总结

本文采用一种构造性、基于概率的方法，对诺斯-胡贝尔分子动力学进行了教学性推导，将阻尼系数 ζ 视为扩展相空间中的动态变量。通过施加李想连续性方程，该方法确保了系综分布 f ∝ e^{-(q²+p²+ζ²)/2} 得到保持，提供了一种比诺斯原始的时间缩放形式更简单、更通用的替代方案。其主要贡献在于为分子动力学模拟中的系综温度控制提供了一个清晰、统一的框架。

ABSTRACT

"Pedagogical derivations for Nos\'e's dynamics can be developed in two different ways, (i) by starting with a temperature-dependent Hamiltonian in which the variable $s$ scales the time or the mass, or (ii) by requiring that the equations of motion generate the canonical distribution including a Gaussian distribution in the friction coefficient $\zeta$. Nos\'e's papers follow the former approach. Because the latter approach is not only constructive and simple, but also can be generalized to other forms of the equations of motion, we illustrate it here. We begin by considering the probability density $f(q,p,\zeta)$ in an extended phase space which includes $\zeta$ as well as all pairs of phase variables $q$ and $p$. This density $f(q,p,\zeta)$ satisfies the conservation of probability (Liouville's Continuity Equation)" $$(\partial f/\partial t) + \sum (\partial (\dot q f)/\partial q) + \sum (\partial (\dot p f)/\partial p) + \sum (\partial (\dot \zeta f)/\partial \zeta) = 0 \ . $$ The multi-authored ``review''\cite{b1} motivated our quoting the history of Nos\'e and Nos\'e-Hoover mechanics, aptly described on page 31 of Bill's 1986 {\it Molecular Dynamics} book, reproduced above\cite{b2}.

研究动机与目标

• 提供一种构造性、教学性的诺斯-胡贝尔动力学推导，避免时间缩放形式的复杂性。
• 证明通过将阻尼系数 ζ 视为扩展相空间中的动态变量，可以生成系综分布。
• 将此方法与诺斯原始的时间缩放方法进行对比，突出其简洁性和可推广性。
• 澄清文献中关于确定性系综器发展的历史误解。
• 支持将诺斯-胡贝尔力学作为非平衡和混沌系统模拟的基础工具。

提出的方法

• 该方法在包含坐标 q、动量 p 和阻尼系数 ζ 的扩展相空间中应用李想连续性方程：∂f/∂t + ∑(∂(¯q f)/∂q) + ∑(∂(¯p f)/∂p) + ∑(∂(¯ζ f)/∂ζ) = 0。
• 推导出诺斯-胡贝尔运动方程：¯q = p，¯p = -q - ζp，¯ζ = p² - 1，这些方程保持系综分布 f ∝ e^{-(q²+p²+ζ²)/2}。
• 该方法通过直接引入 ζ 作为具有自身运动方程的动态变量，避免了时间缩放。
• 通过验证在所推导的运动方程下连续性方程成立，来确认稳态分布的正确性。
• 通过一个包含四个固定散射体和光滑势 φ(r < 1) = (1 - r²)² 的二维单元模型说明该方法。
• 使用四阶龙格-库塔法进行数值积分，模拟轨迹，验证能量守恒和相空间分布的一致性。

实验结果

研究问题

• RQ1 如何在不依赖时间缩放形式的前提下，构造性地推导出系综分布？
• RQ2 阻尼系数 ζ 作为动态变量在生成系综系综中起什么作用？
• RQ3 为何诺斯-胡贝尔方法比原始诺斯时间缩放方法更简单且更具通用性？
• RQ4 连续性方程如何确保在扩展相空间中保持系综分布？
• RQ5 该方法对模拟非平衡和混沌系统具有何种意义？

主要发现

• 诺斯-胡贝尔运动方程：¯q = p，¯p = -q - ζp，¯ζ = p² - 1，精确保持系综分布 f ∝ e^{-(q²+p²+ζ²)/2}。
• 在这些动力学下，连续性方程成立，证实了扩展相空间中概率守恒。
• 与诺斯的时间缩放方法相比，该方法提供了更简单、更直接的系综动力学推导。
• 数值模拟显示 200,000 个时间步长下能量守恒精度达到十一位小数，验证了方程的有效性。
• 该方法通过澄清基于 ζ 的形式比时间缩放更简单、更具通用性，解决了文献中的历史混淆问题。
• 该方法使混沌和非平衡系统（包括六维相空间中具有奇异吸引子的系统）的稳健模拟成为可能。