Roux Benoît
Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago, Illinois 60637, United States.
Department of Chemistry, The University of Chicago, 5735 S. Ellis Avenue, Chicago, Illinois 60637, United States.
J Phys Chem A. 2021 Sep 2;125(34):7558-7571. doi: 10.1021/acs.jpca.1c04110. Epub 2021 Aug 18.
The kinetics of a dynamical system comprising two metastable states is formulated in terms of a finite-time propagator in phase space (position and velocity) adapted to the underdamped Langevin equation. Dimensionality reduction to a subspace of collective variables yields familiar expressions for the propagator, committor, and steady-state flux. A quadratic expression for the steady-state flux between the two metastable states can serve as a robust variational principle to determine an optimal approximate committor expressed in terms of a set of collective variables. The theoretical formulation is exploited to clarify the foundation of the string method with swarms-of-trajectories, which relies on the mean drift of short trajectories to determine the optimal transition pathway. It is argued that the conditions for Markovity within a subspace of collective variables may not be satisfied with an arbitrary short time-step and that proper kinetic behaviors appear only when considering the effective propagator for longer lag times. The effective propagator with finite lag time is amenable to an eigenvalue-eigenvector spectral analysis, as elaborated previously in the context of position-based Markov models. The time-correlation functions calculated by swarms-of-trajectories along the string pathway constitutes a natural extension of these developments. The present formulation provides a powerful theoretical framework to characterize the optimal pathway between two metastable states of a system.
一个由两个亚稳态组成的动力学系统的动力学,是根据相空间(位置和速度)中适应欠阻尼朗之万方程的有限时间传播子来表述的。将维度降低到集体变量子空间,可得到传播子、反应坐标和稳态通量的常见表达式。两个亚稳态之间稳态通量的二次表达式,可作为一种稳健的变分原理,用于确定用一组集体变量表示的最优近似反应坐标。利用该理论公式来阐明具有轨迹群的弦方法的基础,该方法依赖于短轨迹的平均漂移来确定最优过渡路径。有人认为,在集体变量子空间内,任意短时间步长可能不满足马尔可夫性条件,只有在考虑较长滞后时间的有效传播子时才会出现适当的动力学行为。具有有限滞后时间的有效传播子适用于特征值 - 特征向量谱分析,如之前在基于位置的马尔可夫模型背景下所阐述的那样。沿着弦路径由轨迹群计算的时间相关函数构成了这些进展的自然扩展。本公式提供了一个强大的理论框架,用于表征系统两个亚稳态之间的最优路径。