[论文解读] Stable and accurate schemes for Langevin dynamics with general kinetic energies

Gabriel Stoltz, Zofia Trstanova|arXiv (Cornell University)|Sep 9, 2016

Markov Chains and Monte Carlo Methods参考文献 3被引用 5

一句话总结

本文提出了一种稳定且弱二阶精度的数值格式，用于具有非标准动能形式的朗之万动力学，例如自适应约束或非全局Lipschitz形式的动能。该方法通过结合分裂积分器与Metropolis校正，确保了精度与稳定性。该方法可实现对Boltzmann-Gibbs分布的有效采样，并减少复杂系统中的亚稳态行为。

ABSTRACT

We study Langevin dynamics with a kinetic energy different from the standard, quadratic one in order to accelerate the sampling of the Boltzmann-Gibbs distribution. We consider two cases: kinetic energies which are local perturbations of the standard kinetic energy around the origin, where they vanish (this corresponds to the so-called adaptively restrained Langevin dynamics); and more general non-globally Lipschitz energies. We develop numerical schemes which are stable and of weak order two, by considering splitting strategies where the discretizations of the fluctuation/dissipation are corrected by a Metropolis procedure. We use the newly developped schemes for two applications: optimizing the shape of the kinetic energy for the adaptively restrained Langevin dynamics, and reducing the metastability of some toy models with non-globally Lipschitz kinetic energies.

研究动机与目标

• 为在朗之万动力学中使用非标准动能形式时，高效采样Boltzmann-Gibbs分布提供解决方案。
• 开发在非全局Lipschitz或标准二次形式的局部扰动下仍保持稳定与准确的数值格式。
• 通过数值实验，优化自适应约束朗之万动力学中动能形式的形状，以提升采样性能。
• 通过改进的数值积分方法，减少具有非全局Lipschitz动能形式的模型中的亚稳态行为。

提出的方法

• 设计分裂积分器，将动力学分解为确定性与随机分量，并对涨落/耗散项施加校正。
• 在每个积分步骤后应用Metropolis-Hastings校正，以保证细致平衡并确保弱二阶精度。
• 通过精心平衡分裂与校正步骤，构建即使在动能具有奇点或非Lipschitz行为时也保持稳定的格式。
• 利用所提出的格式，通过数值实验优化自适应约束朗之万动力学中动能函数的形式。
• 将该格式应用于具有非全局Lipschitz动能形式的玩具模型，以证明其可减少亚稳态并实现更快收敛。

实验结果

研究问题

• RQ1 当动能非全局Lipschitz时，如何使朗之万动力学的数值格式保持稳定且具有弱二阶精度？
• RQ2 在非标准动能形式下，使用Metropolis校正确对分裂积分器的稳定性与精度有何影响？
• RQ3 是否可通过优化自适应约束朗之万动力学中动能形状来加速Boltzmann-Gibbs分布的采样？
• RQ4 所提出的格式在具有非全局Lipschitz动能形式的模型中，能在多大程度上减少亚稳态行为？

主要发现

• 所提出的Metropolis校正分裂格式实现了弱二阶精度，并在非全局Lipschitz动能形式下仍保持稳定。
• 该格式在具有病态动能函数的玩具模型中成功减少了亚稳态，实现了更快的平衡收敛。
• 通过数值实验表明，自适应约束朗之万动力学中动能形状的优化可显著提升采样效率。
• 该方法保持了目标Boltzmann-Gibbs分布的不变测度，确保了长期统计行为的正确性。