[論文レビュー] Position-dependent memory kernel in generalized Langevin equations: theory and numerical estimation

Hadrien Vroylandt, Pierre Monmarché|arXiv (Cornell University)|Jan 7, 2022

Theoretical and Computational Physics参考文献 50被引用数 45

ひとこと要約

本稿は、平衡ハミルトニアン系における非線形で粗粒度化された力学の一般化ランジュバン方程式（GLE）を、位置に依存する記憶核を用いて厳密に導出する。位置に依存するフラクチュアーション・ディスシペーショング関係を確立し、すべての原子シミュレーションから核を推定するボルテラに基づく数値的手法を導入し、位置に依存する核が、標準的仮定を越えて正確なモデリングに不可欠であることを示している。

ABSTRACT

Generalized Langevin equations with non-linear forces and position-dependent linear friction memory kernels, such as commonly used to describe the effective dynamics of coarse-grained variables in molecular dynamics, are rigorously derived within the Mori-Zwanzig formalism. A fluctuation-dissipation theorem relating the properties of the noise to the memory kernel is shown. The derivation also yields Volterra-type equations for the kernel, which can be used for a numerical parametrization of the model from all-atom simulations.

研究の動機と目的

• モーリ・ツァンツィグ形式内において、非線形力および位置に依存する線形摩擦核を有する一般化ランジュバン方程式（GLE）を厳密に導出すること。
• ノイズの共分散を記憶核に関連付ける位置に依存するフラクチュエーション・ディスシペーショング関係（FDT）を確立すること。
• すべての原子シミュレーションデータから記憶核を推定するための、ボルテラ型積分微分方程式に基づく数値スキームを開発すること。
• 系に位置に依存する摩擦が存在する場合、時間のみに依存する核を仮定する標準的なGLEが、正確な力学を捉えられないことを示すこと。
• 1次元系における数値実験を用いて、スプライン、最小、線形GLEの各形式を比較し、モデルを検証すること。

提案手法

• 集合変数の集合への射影を用いて、モーリ・ツァンツィグ形式を用いて位置に依存する核を有するGLEを導出する。
• 射影形式から導かれるボルテラ型積分微分方程式を満たす記憶核を定式化する。
• すべての原子シミュレーションの時系列データから核を推定するための数値的逆ボルテラ法を実装する。
• フラクチュエーション・ディスシペーショング関係を用いて、ノイズの共分散を記憶核に関連付け、平衡統計力学と整合性を保つ。
• 平均力ポテンシャルを有する1次元モデル系にこの手法を適用し、スプライン、最小、線形GLEの各形式の結果を比較する。
• ノイズの自己相関関数と核の定数部を比較することで、ノイズモデルを検証し、ノイズヒストグラムの非ガウス的尾部を評価する。

実験結果

リサーチクエスチョン

• RQ1 非線形力および位置に依存する摩擦核を有する一般化ランジュバン方程式は、モーリ・ツァンツィグ形式から厳密に導出可能か？
• RQ2 このような系において、ノイズの共分散を記憶核に関連付ける位置に依存するフラクチュエーション・ディスシペーショング関係が成立するか？
• RQ3 ボルテラ型方程式を用いて、すべての原子シミュレーションデータから記憶核を正確に推定できるか？
• RQ4 スプライン、最小、線形GLEといった異なるGLE形式は、真の力学およびノイズ統計を再構築する上でどのように比較されるか？
• RQ5 ノイズ分散の位置依存性が、FDTの有効性および全体的なモデル精度にどの程度影響を与えるか？

主な発見

• 非線形で粗粒度化された力学を正確に記述するためには、GLEにおける記憶核が位置に依存する必要がある。これは、標準的な時間のみに依存するという仮定に反する。
• 平衡統計力学と整合的な形で、ノイズ共分散を記憶核に関連付ける位置に依存するフラクチュエーション・ディスシペーショング関係が厳密に導出された。
• ボルテラ方程式に基づく数値的推定手法により、すべての原子軌道から記憶核を効果的に回復できた。特にスプラインGLEが参照データと最も良好な一致を示した。
• 最小GLEおよびスプラインGLEにおけるノイズは非ガウス的尾部を示しており、単純なガウス的ノイズモデルでは正確な力学を再現できないことが示された。
• 最小GLEの記憶核には定数係数が存在しないため、標準的なFDTと整合しない。一方、スプラインおよび線形GLEでは、ノイズ自己相関関数と核係数の間に良好な一致が得られた。
• ノイズの分散は、ノイズデータをビン別にした結果、強く位置に依存することが確認され、位置に依存するノイズモデリングの必要性が強調された。