当前位置: X-MOL 学术Appl. Comput. Harmon. Anal. › 论文详情
Our official English website, www.x-mol.net, welcomes your feedback! (Note: you will need to create a separate account there.)
Data-driven efficient solvers for Langevin dynamics on manifold in high dimensions
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 MRp based on collected sample points {xi}i=1nM that probe the unknown manifold M. Through the diffusion map, we first learn the reaction coordinates {yi}i=1nN corresponding to {xi}i=1n, where N is a manifold diffeomorphic to M and isometrically embedded in R with p. The induced Langevin dynamics on N 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 N, we leverage the corresponding Fokker-Planck equation on the manifold N 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 N. Furthermore, we provide a weighted L2 convergence analysis of the finite volume scheme to the Fokker-Planck equation on N. The proposed finite volume scheme leads to a Markov chain on {yi}i=1n 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 N and the approximated Markov process enjoys detailed balance, ergodicity, and other good properties.



中文翻译:

高维流形上朗之万动力学的数据驱动高效求解器

我们研究了具有多种结构的物理系统的朗之万动力学Rp基于收集的样本点{X一世}一世=1n探索未知的流形. 通过扩散图,我们首先学习反应坐标{是的一世}一世=1nñ对应于{X一世}一世=1n, 在哪里ñ是流形微分同胚于并等距嵌入Rp. 诱导朗之万动力学ñ就反应坐标而言,它捕获了慢时间尺度的动力学,例如生化反应中的构象变化。为朗之万动力学构建一个有效且稳定的近似ñ,我们在流形上利用相应的 Fokker-Planck 方程ñ根据反应坐标y。我们为这个 Fokker-Planck 方程提出了一个可实现的、无条件稳定的、数据驱动的有限体积方案,它自动结合了ñ. 此外,我们提供加权大号2有限体积方案对 Fokker-Planck 方程的收敛性分析ñ. 所提出的有限体积方案导致马尔可夫链{是的一世}一世=1n在最近邻点之间具有近似的转移概率和跳跃率。在无条件稳定的显式时间离散化之后,数据驱动的有限体积方案给出了一个近似于朗之万动力学的马尔可夫过程ñ近似马尔可夫过程具有详细的平衡性、遍历性和其他良好性质。

更新日期:2022-09-26
down
wechat
bug