Applied and Computational Harmonic Analysis ( IF 2.5 ) Pub Date : 2022-09-26 , DOI: 10.1016/j.acha.2022.09.003 Yuan Gao, Jian-Guo Liu, Nan Wu
We study the Langevin dynamics of a physical system with manifold structure based on collected sample points that probe the unknown manifold . Through the diffusion map, we first learn the reaction coordinates corresponding to , where is a manifold diffeomorphic to and isometrically embedded in with . The induced Langevin dynamics on in terms of the reaction coordinates captures the slow time scale dynamics such as conformational changes in biochemical reactions. To construct an efficient and stable approximation for the Langevin dynamics on , we leverage the corresponding Fokker-Planck equation on the manifold in terms of the reaction coordinates y. We propose an implementable, unconditionally stable, data-driven finite volume scheme for this Fokker-Planck equation, which automatically incorporates the manifold structure of . Furthermore, we provide a weighted convergence analysis of the finite volume scheme to the Fokker-Planck equation on . The proposed finite volume scheme leads to a Markov chain on with an approximated transition probability and jump rate between the nearest neighbor points. After an unconditionally stable explicit time discretization, the data-driven finite volume scheme gives an approximated Markov process for the Langevin dynamics on and the approximated Markov process enjoys detailed balance, ergodicity, and other good properties.
中文翻译:
高维流形上朗之万动力学的数据驱动高效求解器
我们研究了具有多种结构的物理系统的朗之万动力学基于收集的样本点探索未知的流形. 通过扩散图,我们首先学习反应坐标对应于, 在哪里是流形微分同胚于并等距嵌入和. 诱导朗之万动力学就反应坐标而言,它捕获了慢时间尺度的动力学,例如生化反应中的构象变化。为朗之万动力学构建一个有效且稳定的近似,我们在流形上利用相应的 Fokker-Planck 方程根据反应坐标y。我们为这个 Fokker-Planck 方程提出了一个可实现的、无条件稳定的、数据驱动的有限体积方案,它自动结合了. 此外,我们提供加权有限体积方案对 Fokker-Planck 方程的收敛性分析. 所提出的有限体积方案导致马尔可夫链在最近邻点之间具有近似的转移概率和跳跃率。在无条件稳定的显式时间离散化之后,数据驱动的有限体积方案给出了一个近似于朗之万动力学的马尔可夫过程近似马尔可夫过程具有详细的平衡性、遍历性和其他良好性质。