Trstanova Z, Leimkuhler B, Lelièvre T
School of Mathematics, University of Edinburgh, James Clerk Maxwell Building, Peter Guthrie Tait Road, Edinburgh EH9 3FD, UK.
CERMICS (ENPC), INRIA, Marne-la-Vallée 77455, France.
Proc Math Phys Eng Sci. 2020 Jan;476(2233):20190036. doi: 10.1098/rspa.2019.0036. Epub 2020 Jan 15.
Diffusion maps approximate the generator of Langevin dynamics from simulation data. They afford a means of identifying the slowly evolving principal modes of high-dimensional molecular systems. When combined with a biasing mechanism, diffusion maps can accelerate the sampling of the stationary Boltzmann-Gibbs distribution. In this work, we contrast the local and global perspectives on diffusion maps, based on whether or not the data distribution has been fully explored. In the global setting, we use diffusion maps to identify metastable sets and to approximate the corresponding committor functions of transitions between them. We also discuss the use of diffusion maps the metastable sets, formalizing the locality via the concept of the quasi-stationary distribution and justifying the convergence of diffusion maps within a local equilibrium. This perspective allows us to propose an enhanced sampling algorithm. We demonstrate the practical relevance of these approaches both for simple models and for molecular dynamics problems (alanine dipeptide and deca-alanine).
扩散映射从模拟数据中近似朗之万动力学的生成器。它们提供了一种识别高维分子系统缓慢演化的主模式的方法。当与一种偏置机制相结合时,扩散映射可以加速平稳玻尔兹曼-吉布斯分布的采样。在这项工作中,我们基于是否充分探索了数据分布,对比了扩散映射的局部和全局视角。在全局设置中,我们使用扩散映射来识别亚稳态集,并近似它们之间跃迁的相应反应坐标函数。我们还讨论了扩散映射在亚稳态集方面的应用,通过准平稳分布的概念形式化局部性,并证明扩散映射在局部平衡内的收敛性。这种视角使我们能够提出一种增强采样算法。我们展示了这些方法对于简单模型和分子动力学问题(丙氨酸二肽和十肽丙氨酸)的实际相关性。