[论文解读] Bootstrapping time correlation functions of molecular dynamics

Nicolas Desbiens, Philippe Arnault+3 more|arXiv (Cornell University)|Jul 20, 2021

Spectroscopy and Quantum Chemical Studies参考文献 66被引用 5

一句话总结

本文提出了一种计算高效的自展法（BS）方法，用于量化分子动力学模拟中时间相关函数的不确定性，尤其针对二元离子混合物中速度和互扩散电流的自相关函数（ACFs）。通过块重采样方法对时间起点和粒子身份进行重采样，该方法在无需长时间模拟或重复运行的情况下，提供了可靠的误差条，从而仅通过一次模拟即可稳健评估Darken近似和有限尺寸效应。

ABSTRACT

Molecular dynamics is often considered as a numerical experiment. The error bars on the results are therefore mandatory, but sometimes difficult to determine and computationally demanding. As a low-cost approach, we describe the application of the bootstrap (BS) method to the quantification of uncertainties pertaining to the time correlation functions. We chose the autocorrelation functions of velocity and interdiffusion current for a binary ionic mixture as a test bed, and we assessed the merit of the Darken approximation relating both of them. The intrinsic errors related to phase space sampling is investigated comparing the BS method with the reference method of replica. We also study how the BS method can assist in addressing the finite size effects.

研究动机与目标

• 开发一种计算高效的分子动力学模拟时间相关函数不确定性量化方法。
• 通过不确定性量化评估Darken近似在互扩散系数与自扩散系数之间关系的可靠性。
• 利用自展法导出的误差条作为度量，研究分子动力学模拟中的有限尺寸效应。
• 为ab initio和经典MD模拟提供一种实用的后处理不确定性量化工具，尤其适用于无法进行重复模拟的情况。

提出的方法

• 通过重采样连续数据点块，将移动块自展法适配至时间序列数据，以保持自相关函数（ACFs）中的时间相关性。
• 对时间起点（块大小α = 1）和粒子身份进行自展重采样，以估计集体可观测量（如互扩散电流）的不确定性。
• 在重采样粒子时，对每个自展样本进行重归一化，以校正集体变量估计器中的偏差。
• 将自展法的误差估计与参考的重复运行方法进行比较，以验证其准确性和收敛性。
• 利用自展法生成的自相关函数分布计算标准误差，并评估如Darken关系等近似方法的统计显著性。
• 通过分析自展误差条随系统尺寸（N）的变化，监测向热力学极限收敛的过程。

实验结果

研究问题

• RQ1 自展法是否仅通过一次模拟即可为分子动力学中的时间相关函数提供可靠的不确定性估计？
• RQ2 自展法在量化自相关函数中固有的采样误差方面，与重复法相比表现如何？
• RQ3 自展误差条在多大程度上能够检测并量化分子动力学模拟中的有限尺寸效应？
• RQ4 当不确定性被准确量化时，互扩散的Darken近似是否可靠？
• RQ5 基于自展法的误差估计能否指导模拟设计的改进，如增大系统尺寸或延长模拟时间？

主要发现

• 自展法在时间相关函数上提供的误差条与参考的重复运行方法高度一致，验证了其在不确定性量化中的可靠性。
• 归一化互扩散电流自相关函数（JACF）在短时间延迟下的内在误差显著降低，与Zwanzig的理论预测一致。
• 在所有耦合强度（Γ = 0.4, 4, 40）下，Darken近似均位于JACF的1–2个标准误差范围内，证实其有效性，尽管在弱耦合（Γ = 0.4）时误差条更大。
• 有限尺寸效应在速度自相关函数（VACFs）中比在JACFs中更显著，且当系统尺寸从N = 250增加时，自展误差条保持稳定。
• 在集体可观测量中重采样粒子时，必须对每个自展样本进行重归一化，以避免偏差，这是实现准确不确定性估计的关键技术调整。
• 自展法使仅通过一次模拟即可直接监测向热力学极限的收敛过程，为昂贵的重复模拟或长时间模拟方法提供了实用的替代方案。