[论文解读] Driven Langevin dynamics: heat, work and pseudo-work

David A. Sivak, John D. Chodera+1 more|arXiv (Cornell University)|Jul 14, 2011

Advanced Thermodynamics and Statistical Mechanics参考文献 9被引用 2

一句话总结

本文在有限时间步长的朗之万动力学中引入了“影子功”（shadow work）的概念，以校正由数值离散化引起的误差。通过将积分器视为一种受驱的非平衡过程，该方法利用非平衡功定理实现精确的自由能估算，通过影子功校正消除了与时间步长相关的偏差。

ABSTRACT

Common algorithms for computationally simulating Langevin dynamics must discretize the stochastic differential equations of motion. These resulting finite time step integrators necessarily have several practical issues in common: Microscopic reversibility is violated, the sampled stationary distribution differs from the desired equilibrium distribution, and the work accumulated in nonequilibrium simulations is not directly usable in estimators based on nonequilibrium work theorems. Here, we show that even with a time-independent Hamiltonian, finite time step Langevin integrators can be thought of as a driven, nonequilibrium physical process. Once an appropriate work-like quantity is defined -- here called the shadow work -- recently developed nonequilibrium fluctuation theorems can be used to measure or correct for the errors introduced by the use of finite time steps. In particular, we demonstrate that amending estimators based on nonequilibrium work theorems to include this shadow work removes the time step dependent error from estimates of free energies. We also quantify, for the first time, the magnitude of deviations between the sampled stationary distribution and the desired equilibrium distribution for equilibrium Langevin simulations of solvated systems of varying size. While these deviations can be large, they can be eliminated altogether by Metropolization or greatly diminished by small reductions in the time step. Through this connection with driven processes, further developments in nonequilibrium fluctuation theorems can provide additional analytical tools for dealing with errors in finite time step integrators.

研究动机与目标

• 解决有限时间步长朗之万积分器违反微观可逆性并采样错误平稳分布的根本问题。
• 识别并量化在溶剂化系统中采样的平稳分布与真实平衡分布之间的偏差。
• 开发一种类似功的量——影子功，以校正非平衡功定理中的时间步长误差。
• 证明在估算器中引入影子功可消除自由能计算中与时间步长相关的误差。
• 探讨通过Metropolization或减小时间步长，能否消除或最小化平衡模拟中的采样误差。

提出的方法

• 定义一种新的类似功的量，称为“影子功”，以量化有限时间步长离散化在朗之万动力学中引入的非平衡驱动效应。
• 将最近发展的非平衡涨落定理应用于影子功，以校正自由能估算器中的偏差。
• 通过不同尺寸的溶剂化系统的数值模拟，量化平稳分布偏离平衡状态的程度。
• 比较标准积分器与Metropolized版本的性能，以评估校正协议对误差减少的效果。
• 推导并分析由时间步长离散化诱导的有效非平衡动力学，将其表述为一种受驱过程。
• 证明影子功可实现与时间步长大小无关的一致自由能估算。

实验结果

研究问题

• RQ1 在哈密顿量不显含时间的情况下，朗之万动力学中的有限时间步长离散化如何引入非平衡效应？
• RQ2 在溶剂化系统中，采样的平稳分布与真实平衡分布之间的偏差程度有多大？
• RQ3 能否构建一个明确定义的类似功的量——影子功，以校正非平衡功定理中由时间步长引起的误差？
• RQ4 影子功在多大程度上可消除自由能估算器中与时间步长相关的偏差？
• RQ5 Metropolization和减小时间步长在缓解有限时间步长引入的采样误差方面有多有效？

主要发现

• 即使哈密顿量不显含时间，有限时间步长朗之万积分器仍会引入非平衡、受驱的过程，这是由于数值离散化所致。
• 在平衡模拟中，采样的平稳分布可能与真实平衡分布显著偏离，尤其在较小的溶剂化系统中更为明显。
• 对不同尺寸的溶剂化系统，量化了这些偏差的大小，表明在未校正时存在显著误差。
• Metropolization通过在积分器中强制满足细致平衡，完全消除了采样误差。
• 减小时间步长可大幅降低采样误差，但无法完全消除。
• 在非平衡功定理中引入影子功，成功消除了自由能估算中与时间步长相关的偏差。