[論文レビュー] The projection of a nonlocal mechanical system onto the irreversible generalized Langevin equation, II: Numerical simulations

Marc Vogt, Rigoberto Hernandez|arXiv (Cornell University)|Jan 14, 2005

Spectroscopy and Quantum Chemical Studies被引用数 10

ひとこと要約

この論文は、非定常摩擦を確率的力学でモデル化する非可逆一般化ランジュバン方程式（iGLE）が、決定論的開放ハミルトニアン系への射影として物理的に実現可能であることを示している。数値シミュレーションにより、ハミルトニアン系における時間依存エネルギーおよび溶媒力相関がiGLEと正確に一致することが確認され、iGLEが現象論的構築物ではなく物理的根拠を持つモデルであることが検証された。非局所的散乱項がこの一致に不可欠であることが示された。

ABSTRACT

The irreversible generalized Langevin equation (iGLE) contains a nonstationary friction kernel that in certain limits reduces to the GLE with space-dependent friction. For more general forms of the friction kernel, the iGLE was previously shown to be the projection of a mechanical system with a time-dependent Hamiltonian. [R. Hernandez, J. Chem. Phys. 110, 7701 (1999)] In the present work, the corresponding open Hamiltonian system is further explored. Numerical simulations of this mechanical system illustrate that the time dependence of the observed total energy and the correlations of the solvent force are in precise agreement with the projected iGLE.

研究の動機と目的

• 非可逆一般化ランジュバン方程式（iGLE）を、時間依存ハミルトニアンを持つ決定論的開放ハミルトニアン系から射影することで、物理的に実現可能なモデルとして確立すること。
• iGLEの非定常摩擦カーネルおよびフラクチュエーション・ドレイン理論関係（FDR）が、時間依存ハミルトニアンを持つ機械的系から自然に導かれるかを検証すること。
• iGLEにおける非局所的散乱項が、基礎となるハミルトニアン系の力学を再現するために不可欠であることを示すこと。
• 全機械的系の数値積分によりiGLEの確率的挙動を検証し、速度自己相関関数および力相関関数において一致を示すこと。
• 現象論的モデルを超えて、非定常確率的力学の物理的基盤を拡張し、決定論的ハミルトニアンフレームワークと同等であることを示すこと。

提案手法

• iGLEは、摩擦カーネルが時間関数g(t)で変調される時間依存ハミルトニアンを持つ機械的系の射影として導出された。
• 系は、200個の調和振動子からなるバーストと結合したトレース粒子から構成され、全ハミルトニアンはiGLEの力学を正確に再現するように明示的に構築された。
• iGLEにおける非局所的散乱項は、バーストの時間発展する結合から導かれるが、正確なシミュレーションにはその明示的評価が必要である。
• エネルギーと長時間にわたる力学的挙動を保存するため、シンプレクティック積分器を用いて全ハミルトニアンの運動方程式を数値的に統合した。
• 各タイムステップで、FDRを満たすために、ξ(t) = v̇ + γ₀g(t)z₂(t) − F(t) の関係によりランダム力ξ(t)を計算した。
• 速度自己相関関数および力相関関数を、ハミルトニアン系と射影されたiGLEの間で直接比較した。

実験結果

リサーチクエスチョン

• RQ1 非可逆一般化ランジュバン方程式（iGLE）は、時間依存ハミルトニアンを持つ決定論的機械的系として物理的に実現可能か？
• RQ2 iGLEにおける非局所的散乱項は、基礎となるハミルトニアン系の正しい力学を再現するために必須の役割を果たすか？
• RQ3 ハミルトニアン系における速度自己相関関数および溶媒力相関関数は、iGLEが予測するものと定量的に同一か？
• RQ4 非定常条件下でも、ハミルトニアン系におけるランダム力に対してフラクチュエーション・ドレイン理論関係（FDR）が成立するか？
• RQ5 初期の摩擦強度およびg(t)の時間的変化が、ハミルトニアン系とiGLEとの一致に与える影響は何か？

主な発見

• 数値シミュレーションにより、すべてのケースおよび時間スケールにおいて、ハミルトニアン系におけるトレース粒子の速度自己相関関数とiGLEが予測する関数との間に正確な一致が確認された。
• ハミルトニアン系における溶媒力の自己相関関数は、すべての時間およびすべてのパrameter選択において、フラクチュエーション・ドレイン理論関係（FDR）を満たしている。
• 非局所的散乱項は不可欠である：その省略は、ハミルトニアン系とiGLEとの一致を破壊する。シミュレーションにより確認された。
• 自由粒子の場合、初期時刻（t = -4.0）における力学は、線形な速度自己相関関数によって示されるように、球状運動的挙動を示し、iGLEの予測と整合的である。
• 調和井戸の場合、閉じ込め効果により摩擦変化に対してやや感受性が低いが、非局所項を含めるとiGLEは依然として正確に力学を再現している。
• すべてのケースが長時間（t = 4.0）に同じ平衡状態に達しており、iGLEが系の長時間統計的挙動を正しく捉えていることが確認された。