[論文レビュー] Driven Langevin dynamics: heat, work and pseudo-work

David A. Sivak, John D. Chodera+1 more|arXiv (Cornell University)|Jul 14, 2011

Advanced Thermodynamics and Statistical Mechanics参考文献 9被引用数 2

ひとこと要約

この論文は、数値離散化に起因する誤差を補正するため、有限時間刻みのランジュバンダイナミクスにおける「シャドウ仕事」の概念を導入する。積分子を駆動される非平衡過程とみなすことにより、非平衡仕事定理を用いて正確な自由エネルギー推定が可能となり、シャドウ仕事による補正によって時間刻み依存のバイアスが解消される。

ABSTRACT

Common algorithms for computationally simulating Langevin dynamics must discretize the stochastic differential equations of motion. These resulting finite time step integrators necessarily have several practical issues in common: Microscopic reversibility is violated, the sampled stationary distribution differs from the desired equilibrium distribution, and the work accumulated in nonequilibrium simulations is not directly usable in estimators based on nonequilibrium work theorems. Here, we show that even with a time-independent Hamiltonian, finite time step Langevin integrators can be thought of as a driven, nonequilibrium physical process. Once an appropriate work-like quantity is defined -- here called the shadow work -- recently developed nonequilibrium fluctuation theorems can be used to measure or correct for the errors introduced by the use of finite time steps. In particular, we demonstrate that amending estimators based on nonequilibrium work theorems to include this shadow work removes the time step dependent error from estimates of free energies. We also quantify, for the first time, the magnitude of deviations between the sampled stationary distribution and the desired equilibrium distribution for equilibrium Langevin simulations of solvated systems of varying size. While these deviations can be large, they can be eliminated altogether by Metropolization or greatly diminished by small reductions in the time step. Through this connection with driven processes, further developments in nonequilibrium fluctuation theorems can provide additional analytical tools for dealing with errors in finite time step integrators.

研究の動機と目的

• 有限時間刻みのランジュバン積分子が微視的可逆性を破り、誤った定常分布をサンプリングするという根本的問題に対処すること。
• 溶媒化系におけるサンプリングされた定常分布と真の平衡分布との間のずれを特定および定量すること。
• 時間刻みに起因する誤差を補正可能な「シャドウ仕事」と呼ばれる、仕事に類似した量を考案すること。
• 推定器にシャドウ仕事を取り入れることで、自由エネルギー計算における時間刻み依存の誤差が除去されることを実証すること。
• メトロポリス化または時間刻みの短縮が、平衡シミュレーションにおけるサンプリング誤差をどのように排除または最小限に抑えるかを検討すること。

提案手法

• 有限時間刻みの離散化によって生じる非平衡駆動効果を定量化する新たな仕事に類似した量、すなわち「シャドウ仕事」を定義する。
• 最近開発された非平衡フラクチュエーション定理をシャドウ仕事に適用し、自由エネルギー推定器におけるバイアスを補正する。
• サイズの異なる溶媒化系の数値シミュレーションを用いて、定常分布の平衡からのずれの大きさを定量する。
• 標準的積分子とメトロポリス化されたバージョンの性能を比較し、補正プロトコルによる誤差低減を評価する。
• 時間刻みの離散化が誘発する有効な非平衡ダイナミクスを導出し、それを駆動過程として定式化する。
• シャドウ仕事が時間刻みの大きさに依存しない一貫性のある自由エネルギー推定を可能にすることを示す。

実験結果

リサーチクエスチョン

• RQ1 時間に依存しないハミルトニアンを持つにもかかわらず、ランジュアンダイナミクスにおける有限時間刻みの離散化がどのように非平衡効果を生じるのか？
• RQ2 溶媒化系において、サンプリングされた定常分布と真の平衡分布との間のずれの大きさはどの程度か？
• RQ3 時間刻みに起因する誤差を補正可能な明確に定義された仕事に類似した量、すなわちシャドウ仕事は構築可能か？
• RQ4 シャドウ仕事は、自由エネルギー推定器における時間刻み依存バイアスをどの程度まで除去できるか？
• RQ5 メトロポリス化と時間刻みの短縮は、有限時間刻みに起因するサンプリング誤差をどの程度効果的に低減できるか？

主な発見

• 有限時間刻みのランジュアン積分子は、数値的離散化のため、時間に依存しないハミルトニアンであっても非平衡的で駆動される過程を引き起こす。
• 平衡シミュレーションにおけるサンプリングされた定常分布は、特に小さな溶媒化系では真の平衡分布から著しくずれることがある。
• サイズの異なる溶媒化系に対して、これらのずれの大きさを定量した結果、補正がなければ顕著な誤差が生じることが示された。
• メトロポリス化は、積分子に詳細平衡を強制することで、サンプリング誤差を完全に排除する。
• 時間刻みを短縮することで、サンプリング誤差は著しく小さくなるが、完全には解消されない。
• 非平衡仕事定理にシャドウ仕事を取り入れることで、自由エネルギー推定における時間刻み依存バイアスが成功裏に除去された。