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

技術要約：非線形ポテンシャルおよび観測量依存質量を伴う厳密な一般化ランジュバン方程式の比較

問題提起
モリー・ツィグ（Mori-Zwanzig）射影形式は、複雑な系における任意の粗視化反応座標に対して、厳密な一般化ランジュバン方程式（GLE）を導出するための厳密な枠組みを提供する。得られる方程式は、マルコフ的な力、非マルコフ的な摩擦項、および直交（ランダム）力によって構成されることが普遍的であるが、その各項の具体的な形式は、射影演算子の選択に大きく依存する。有効質量と摩擦が観測量自体に依存し、かつ平均力ポテンシャルが非線形である観測量を扱う場合には、重大な課題が生じる。さらに、メモリカーネルと直交力の関係（第二モーメント関係、またはSMR）は、しばしば成立すると仮定されているが、一般的なケースでは保証されていない。本論文は、異なる厳密なGLE定式化を体系的に比較し、有効ポテンシャルへの運動エネルギーの包含や、摩擦カーネルの観測量依存性に関する特性を決定する必要性に対処するものである。

手法
著者らは、位置の自由度に依存するスカラー観測量 $A(t)$ に対して、4つの異なる厳密なGLEを導出し、比較するためにモリー・ツィグ演算子形式を採用している。導出には、一対の射影演算子（$P_1, Q_1$ および $P_2, Q_2$）を用いた二重射影スキームを用いている。

本研究は、観測量の速度 $\dot{A}$ が条件付きウィックの定理（式21）を満たす特定のケースに焦点を当てている。4つの異なるGLEを生成するために、以下の4つの射影演算子の組み合わせが利用されている：

1. I-GLE: $P_1 = P_{LV}$（線形速度）および $P_2 = P_M$（モリー）を使用。観測量に依存しない摩擦カーネルを特徴とする。
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データと比較する。

主要な貢献および結果

• 4つの厳密なGLEの導出: 本論文は、4つの厳密な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）にも依存（D-GLE, DKE-GLE）にもなり得ることを確認している。
• タンパク質折り畳み（Villin）への適用:
  • 著者らは、これらのGLEをVillinの折り畳みダイナミクスに適用している。彼らは、ネイティブコンタクトの割合という反応座標は、複数の双線形結合（$N_d > 1$）を含むため、厳密にはウィックの定理（式21）を満たさないことを指摘している。
  • この違反があるにもかかわらず、数値結果は、I-GLEおよびIKE-GLEが、MDデータから抽出された同時分布 $P(A, \dot{A})$ および有効質量 $M_{eff}(A)$ を正確に近似できることを示している。
  • 定質量GLE（CM-GLE）は、同時分布を再現できず、 $A$ の周辺分布のみを捉えるにとどまる。
  • 3つのGLE（I-GLE, IKE-GLE, CM-GLE）すべてを用いてMDデータから抽出されたメモリカーネル $\Gamma(t)$ は、極めて類似していることが判明した。
  • 直交力の解析により、CM-GLEが近似的にガウス型のランダム力分布を与える一方で、I-GLEおよびIKE-GLEは、質量の強い観測量依存性を反映した、かなりの非ガウス的特性を示すことが明らかになった。

意義および主張
本論文は、有効的な運動エネルギーをポテンシャルに含めること（IKE-GLEおよびDKE-GLEに見られるように）が、複雑なデータの解析において特に有利であると主張している。観測量がウィックの定理を満たす場合、これらの定式化は、摩擦項やノイズ項に依存することなく、位置と速度の正しい平衡同時分布を保証する。

タンパク質の折り畳みへの適用に関して、著者らは、Villinの反応座標は理論的前提（ウィックの定理）を厳密には満たさないものの、I-GLEおよびIKE-GLEの定式化は、折り畳みダイナミクスを記述するための「優れた近似」として機能すると控えめに主張している。彼らは、平均力ポテンシャル、構成依存質量、およびメモリカーネルといった主要な速度論的量、ならびに $P(A, \dot{A})$ の同時分布を、MDシミュレーションから観察された通りに再現することに成功した。本研究は、厳密な理論的条件が近似的にしか満たされない場合であっても、これらの厳密なGLE定式化が、非線形ポテンシャルや観測量依存質量を持つ系をモデル化するための堅牢なツールであることを示唆している。