[论文解读] Ergodicity of a Time-Reversibly Thermostated Harmonic Oscillator and the 2014 Ian Snook Prize

Wm.G. Hoover, C.G. Hoover|arXiv (Cornell University)|Jan 30, 2020

Advanced Physical and Chemical Molecular Interactions被引用 3

一句话总结

本文通过数值模拟与相空间分析，采用三种动力学模型——Marta-Klein-Tuckerman (MKT)、Hoover-Holian (HH) 和 Patra-Bhattacharya (PB)——研究了时间可逆系缚谐振子的遍历性。研究发现，MKT 与 HH 模型表现出与系综高斯分布一致的遍历行为，而 PB 模型则显示出非遍历特征，如环面结构和密度偏差，挑战了早期关于 MKT 非遍历性的说法，并为此类系统提供了 2014 年 Ian Snook 奖以证明遍历性的最终证据。

ABSTRACT

Shuichi Nos\'e opened up a new world of atomistic simulation in 1984. He formulated a Hamiltonian tailored to generate Gibbs' canonical distribution dynamically. This clever idea bridged the gap between microcanonical molecular dynamics and canonical statistical mechanics. Until then the canonical distribution was explored with Monte Carlo sampling. Nos\'e's dynamical Hamiltonian bridge requires the "ergodic" support of a space-filling structure in order to reproduce the entire distribution. For sufficiently small systems, such as the harmonic oscillator, Nos\'e's dynamical approach failed to agree with Gibbs' sampling and instead showed a complex structure, partitioned into a chaotic sea, islands, and chains of islands, that is familiar textbook fare from investigations of Hamiltonian chaos. In trying to enhance small-system ergodicity several more complicated "thermostated" equations of motion were developed. All were consistent with the canonical Gaussian distribution for the oscillator coordinate and momentum. The ergodicity of the various approaches has undergone several investigations, with somewhat inconclusive ( contradictory ) results. Here we illustrate several ways to test ergodicity and challenge the reader to find even more convincing algorithms or an entirely new approach to this problem.

研究动机与目标

• 测试不同动力学模型下时间可逆系缚谐振子的遍历性。
• 在 Patra 与 Bhattacharya 声称 MKT 非遍历的背景下，解决关于 MKT 与 HH 模型遍历性的相互矛盾的主张。
• 通过时间平均、投影、密度测量和周期轨道分析，提供数值与几何证据，以判断这些模型是否能采样完整的系综分布。
• 向研究社区发起挑战，要求开发更严格或新颖的算法以证明遍历性，并提供 2014 年 Ian Snook 奖作为激励。

提出的方法

• 使用四阶龙格-库塔法，时间步长为 0.001，对 MKT、HH 和 PB 模型进行一千万个时间步长的数值积分，以模拟轨迹。
• 对速度矩（p², p⁴, p⁶）进行时间平均，与系综高斯分布的期望值（1, 3, 15）进行比较。
• 对相空间在 (q,p)、(q,ζ)、(q,ξ)、(p,ζ)、(p,ξ) 和 (ζ,ξ) 平面上的投影，以可视化分布结构并检测空洞或对称性。
• 通过十亿个时间步长，在逐渐缩小的 (q,p) 圆形区域（面积从 0.01π 降至 0.01π/16）内测量密度，以评估原点附近的均匀性。
• 通过原点附近的初始条件识别周期轨道，特别在 PB 模型中发现嵌套环面结构。
• 使用随机或网格化初始条件进行系综测试，以检测双峰分布，表明非遍历性。

实验结果

研究问题

• RQ1 MKT、HH 和 PB 模型是否产生与系综高斯分布一致的时间平均速度矩？
• RQ2 这些模型的相空间投影是否表现出无空洞或不对称性的高斯结构，从而表明遍历性？
• RQ3 在四维相空间 (q,p,ζ,ξ) = (0,0,0,0) 附近，各模型的概率密度是否均匀？
• RQ4 PB 方程是否表现出周期轨道或环面结构，提示非遍历行为？
• RQ5 能否为 MKT 或 HH 模型构建出确定性的数学证明，以解决关于遍历性的相互矛盾主张？

主要发现

• MKT 与 HH 模型的速度矩（p², p⁴, p⁶）为（1.000, 3.000, 15.00），四位有效数字下与系综高斯分布完全一致。
• PB 模型的速度矩为（1.000, 3.800, 23.50），显著偏离高斯行为，表明非遍历性。
• MKT 与 HH 的相空间投影显示混沌海中存在原点附近的微小空洞，而 PB 显示 q 与 p 之间存在强烈正相关性，并呈现明显的环面结构。
• 在 (q,p) 圆形区域半径逐渐减小（降至原始面积的 1/16）的密度测量中，MKT 与 HH 未显示显著变化，表明在原点附近采样均匀。
• PB 模型的密度约为 MKT/HH 的一半，表明其相空间中约 50% 的测度为拟周期性，被排除在混沌采样之外。
• PB 模型在 (q,p,ζ) 空间中表现出清晰的环面结构，表明存在嵌套周期轨道及混沌海中的大空腔，确认了非遍历行为。