The study of rare events in molecular and atomic systems such as conformal changes and cluster rearrangements has been one of the most important research themes in chemical physics. Key challenges are associated with long waiting times rendering molecular simulations inefficient, high dimensionality impeding the use of PDE-based approaches, and the complexity or breadth of transition processes limiting the predictive power of asymptotic methods. Diffusion maps are promising algorithms to avoid or mitigate all these issues. We adapt the diffusion map with Mahalanobis kernel proposed by Singer and Coifman (2008) for the SDE describing molecular dynamics in collective variables in which the diffusion matrix is position-dependent and, unlike the case considered by Singer and Coifman, is not associated with a diffeomorphism. We offer an elementary proof showing that one can approximate the generator for this SDE discretized to a point cloud via the Mahalanobis diffusion map. We use it to calculate the committor functions in collective variables for two benchmark systems: alanine dipeptide, and Lennard-Jones-7 in 2D. For validating our committor results, we compare our committor functions to the finite-difference solution or by conducting a "committor analysis" as used by molecular dynamics practitioners. We contrast the outputs of the Mahalanobis diffusion map with those of the standard diffusion map with isotropic kernel and show that the former gives significantly more accurate estimates for the committors than the latter.
翻译:分子和原子系统中的稀有事件研究,如相近变化和集群重新排列,一直是化学物理中最重要的研究主题之一。关键挑战与漫长的等待时间相关,因为分子模拟效率低,高维阻碍使用基于PDE的方法,以及转型过程的复杂性或广度限制了无药可救方法的预测力。扩散图是避免或减轻所有这些问题的有希望的算法。我们用Singer和Coifman(2008年)提出的马哈拉诺比斯内核的分布图来调整SDE的分布图,描述在集体变量中的分子动态,在这些变量中,扩散矩阵取决于位置,而不像Singer和Coifman所考虑的案例中,分子模拟过程效率高,阻碍使用PDF-DE法的方法,以及转型过程的复杂或宽度,我们提供了一个基本证据,表明SDE的生成者可以通过Mahalanobism映射图接近点云。我们用它来计算两种基准系统中的集体变量中的承诺函数:alanine diptide,和Lennard-Jones-7D,其中的分子7D。为了确认我们之前的模型的流流流化结果,我们用前的流流流化分析,我们用前的流流流流流的流的流分析是用来进行较前的流分析,我们使用的流化的流化的流结果,我们用来进行较前的流的流的流的流的流的流的流的流的流的流结果的流的流的流的流的流。