Comparison of different exact generalized Langevin equations with a non-linear potential of mean force and an observable-dependent mass and friction

技术摘要：具有非线性势能与观测量相关质量的精确广义朗之万方程比较研究

问题陈述
Mori-Zwanzig 投影形式理论为推导复杂系统中任意粗粒化反应坐标的广义朗之万方程（GLEs）提供了一个精确框架。虽然所得方程普遍由马尔可夫力、非马尔可夫摩擦项和正交（随机）力组成，但这些项的具体形式高度取决于投影算子的选择。当处理有效质量和摩擦力依赖于观测量本身，且平均力势（potential of mean-force）具有非线性的观测量时，会出现一个关键挑战。此外，关于记忆核与正交力之间的关系（即第二矩关系，SMR）通常被假设成立，但在一般情况下并不能得到保证。本文旨在系统地比较不同的精确 GLE 公式，以确定其关于有效势能中包含动能以及摩擦核对观测量依赖性的特性。

方法论
作者利用 Mori-Zwongiz 算子形式推导并比较了四种针对依赖于位置自由度的标量观测量 $A(t)$ 的不同精确 GLE。该推导利用了涉及一对投影算子（$P_1, Q_1$ 和 $P_2, Q_2$）的双重投影方案。

本研究重点关注观测量速度 $\dot{A}$ 满足条件维克定理（Wick's theorem，公式 21）的具体情况。研究采用了四种投影算子组合来生成四种不同的 GLE：

1. I-GLE：使用 $P_1 = P_{LV}$（线性速度）和 $P_2 = P_M$（Mori）。其特征是具有与观测量无关的摩擦核。
2. D-GLE：使用 $P_1 = P_2 = P_{LV}$。其特征是具有与观测量相关的摩擦核。
3. IKE-GLE：使用 $P_1 = P_{KE}$（动能）和 $P_2 = P_M$。其特征是具有与观测量无关的摩擦核，并包含了有效动能项在势能中。
4. DKE-GLE：使用 $P_1 = P_{KE}$ 和 $P_2 = P_{LV}$。其特征是具有与观测量相关的摩擦核，并包含了有效动能项在势能中。

作者定义了有效质量 $M_{eff}(A)$ 以及各种有效势能（$U_{eff}$ 和 $U^{KE}_{eff}$）。他们严格分析了记忆核与正交力之间满足第二矩关系（SMR）的条件。最后，这些理论框架被应用于 Villin 蛋白质折叠的分子动力学（MD）模拟数据，并以原生接触比例作为反应坐标。作者将 I-GLE、IKE-GLE 和恒定质量 GLE（CM-GLE）的性能与 MD 数据进行了对比。

主要贡献与结果

• 推导四种精确 GLE：本文系统地推导了四种精确 GLE。研究识别出一个主要的区别在于有效势能的定义：
  • I-GLE 和 D-GLE 使用仅依赖于位置的势能 $U_{eff}(A)$。
  • IKE-GLE 和 DKE-GLE 使用显式包含有效动能项（$\frac{1}{2}M_{eff}(A)\dot{A}^2$）的势能 $U^{KE}_{eff}(A, \dot{A})$。
• 包含动能的优势：作者证明，对于满足维克定理的观测量，IKE-GLE 和 DKE-GLE 公式仅通过有效势能 $U^{KE}_{eff}$ 的玻尔兹曼因子即可重现观测量 $A$ 及其速度 $\dot{A}$ 的精确联合分布。这一过程无需依赖非马尔可夫摩擦或正交力的贡献。这被视为模拟复杂数据的显著优势，因为它确保了能够自然地恢复正确的平衡分布。
• 第二矩关系（SMR）的有效性：研究澄清了 SMR（即记忆核与正交力的自相关函数成正比）在这些 GLE 中并不具有普适性。SMR 仅在限制性条件下有效，例如正交力的条件平均值为零，或者满足特定的偏微分方程。这一发现挑战了将 SMR 失效仅归因于非谐振势能的解释，表明其实际上是投影算子选择的结果。
• 观测量依赖型摩擦：本文确认，根据所选投影算子的不同，摩擦核可以是独立的（I-GLE, IKE-GLE）或依赖于观测量 $A$ 的（D-GLE, DKE-GLE）。
• 蛋白质折叠应用（Villin）：
  • 作者将这些 GLE 应用于 Villin 的折叠动力学。他们指出，原生接触比例这一反应坐标并不严格满足维克定理（公式 21），因为它涉及多个双线性组合的求和（$N_d > 1$）。
  • 尽管存在这种违反，数值结果显示 I-GLE 和 IKE-GLE 能够准确近似从 MD 数据中提取的联合分布 $P(A, \dot{A})$ 和有效质量 $M_{eff}(A)$。
  • 恒定质量 GLE（CM-GLE）无法重现联合分布，只能捕捉到 $A$ 的边缘分布。
  • 从 MD 数据中提取的三个 GLE（I-GLE, IKE-GLE, CM-GLE）的记忆核 $\Gamma(t)$ 被发现是非常相似的。
  • 对正交力的分析表明，虽然 CM-GLE 产生近似高斯分布的随机力，但 I-GLE 和 IKE-GLE 表现出显著的非高斯特性，反映了质量对观测量的强烈依赖性。

意义与主张
本文声称，在势能中包含有效动能（如 IKE-GLE 和 DKE-GLE 所示）对于分析复杂数据特别有利。当观测量满足维克定理时，这些公式保证了能够通过势能直接获得正确的位移-速度联合平衡分布，而无需依赖摩擦项或噪声项。

关于蛋白质折叠的应用，作者谦虚地指出，虽然 Villin 反应坐标并不严格满足理论前提（维克定理），但 I-GLE 和 IKE-GLE 公式仍可作为描述折叠动力学的“良好近似”。它们成功重现了关键动力学量（平均力势、构型相关质量和记忆核）以及 MD 模拟中观察到的联合分布 $P(A, \dot{A})$。这项工作表明，即使在严格的理论条件仅在近似层面满足的情况下，这些精确 GLE 公式仍是建模具有非线性势能和观测量依赖质量系统的鲁棒工具。