[論文レビュー] Energy landscapes and rare events

E Weinan, Weiqing Ren+1 more|ArXiv.org|Dec 1, 2002

stochastic dynamics and bifurcation参考文献 12被引用数 4

ひとこと要約

本稿では、複雑で粗いポテンシャルを有する高次元エネルギー・ランドスケープにおいて、最小エネルギー経路（MEP）および遷移レートを計算するための有限温度ストリング法を導入する。反応座標をストリング（一次元曲線）として表現し、ストリングに沿った局所的平衡分布を組み合わせることで、元の滑らかなポテンシャルの明示的知識がなくても、遷移機構とレートを効率的に同定できる。これにより、材料やバイオ分子系におけるレアイベントの正確な計算が可能になる。

ABSTRACT

Many problems in physics, material sciences, chemistry and biology can be abstractly formulated as a system that navigates over a complex energy landscape of high or infinite dimensions. Well-known examples include phase transitions of condensed matter, conformational changes of biopolymers, and chemical reactions. The energy landscape typically exhibits multiscale features, giving rise to the multiscale nature of the dynamics. This is one of the main challenges that we face in computational science. In this report, we will review the recent work done by scientists from several disciplines on probing such energy landscapes. Of particular interest is the analysis and computation of transition pathways and transition rates between metastable states. We will then present the string method that has proven to be very effective for some truly complex systems in material science and chemistry.

研究の動機と目的

• 多スケール特徴と粗いポテンシャルを有する高次元エネルギー・ランドスケープにおける遷移経路とレートの計算という課題に対処すること。
• 事前に定義された反応座標と滑らかなランドスケープに依存する伝統的な遷移状態理論の限界を克服すること。
• 元の滑らかなポテンシャルの事前知識が不要な、自己整合的に反応座標を同定し、遷移メカニズムを計算する計算フレームワークの構築。
• 時間スケールの大きな差異を示す系、例えばバイオポリマーのコンformational 変化や材料の相転移において、レアイベントの正確なモデリングを可能にすること。
• 熱揺らぎを伴う確率的ダイナミクスにおける遷移レート計算のための、大偏差理論の数値的に効率的な代替手法を提供すること。

提案手法

• 有限温度ストリング法は、状態空間内の一次元曲線（ストリング）を反応座標として用い、ストリング上の点がポテンシャル力の法線成分に従って進化する。
• この手法は、ストリングに垂直な運動を強制する射影作用素、熱揺らぎを表すガウスノイズ、および均等なパrameter化を維持するためのラグランジュ乗数を含む確率微分方程式（5.3）を用いる。
• ストリングの統計的定常状態は、ストリング上の位置が与えられたときの系の位置の条件付き期待値に一致し、これはストリングに直交する超平面におけるボルツマン分布から導かれる。
• この手法は、ストリングを進化させ、各点における局所的平衡分布（μα）を自由エネルギー・ランドスケープを用いて推定することで、有効な最小エネルギー経路（MEP）を計算する。
• スムーズなポテンシャルの明示的計算を避けるために、微小スケールの特徴が熱平均によって平滑化されるという事実を活用し、完全な粗いポテンシャル V(x) を直接扱う。
• ストリングは定期的に再パrameter化され、均等な間隔を維持する。この手法は、摂動を加えたミューラー・ポテンシャルに対して検証され、ノイズが存在する中でも真のMEPを正確に回復することが示された。

実験結果

リサーチクエスチョン

• RQ1 微小スケールの熱揺らぎを含む粗いポテンシャルを有する高次元エネルギー・ランドスケープにおいて、最小エネルギー経路（MEP）をどのように計算できるか。
• RQ2 反応座標が事前に未知である場合に、メタ安定状態間の遷移レートを効率的かつ安定的に計算する数値的手法は何か。
• RQ3 ポテンシャルエネルギー・ランドスケープが熱エネルギー（O(kBT)）と同程度のノイズによって摂動される系において、ストリング法をどのように適応可能にするか。
• RQ4 事前に定義された反応座標や鞍点解析を必要とせずに、ストリング法が自己整合的に反応座標を同定できるか。
• RQ5 エントロピー障壁を有する系において、計算されたストリングと真の自由エネルギー・ランドスケープとの関係は何か。

主な発見

• 有限温度ストリング法は、摂動を加えたミューラー・ポテンシャルに対して、最小エネルギー経路（MEP）を正確に計算できた。計算されたMEP（図4の黒線）は、滑らかな成分の真のMEP（赤線）とよく一致した。
• この手法はストリング周辺の局所的揺らぎを捉えており、図4の緑線で示されるように、ストリング上の各点におけるボルツマン分布の期待される熱幅と整合的である。
• ストリングダイナミクスの統計的定常状態（式5.4）は、ストリング上の位置が与えられたときの位置の条件付き期待値に正確に一致し、統計力学との整合性が確認された。
• ストリングに沿った定常分布は、ストリングに直交する超平面におけるボルツマン加重マージナルに等しい（式5.5）。これにより、自由エネルギーおよび遷移レートの直接計算が可能になる。
• この手法はポテンシャルのノイズに対して頑健であり、|δV| ≤ O(kBT) の範囲で摂動が加えられた場合でも、ポテンシャルの明示的分解が不要なままに正確なMEP回復が達成された。
• 大偏差理論の複雑な境界値問題を解く必要がなく、精度を保ちつつも計算が効率的であるため、この手法はその代替手段として有効である。