[論文レビュー] Exact rate calculations by trajectory parallelization and twisting

Eric Vanden‐Eijnden, Maddalena Venturoli|ArXiv.org|Apr 23, 2009

Quantum chaos and dynamical systems被引用数 9

ひとこと要約

本稿では、遷移経路理論を用いて、平衡および非平衡系において正確な反応速度を計算するための軌道並列化およびねじり手法を提示する。ボロノイタイルレーションとフラックスマッチング境界条件を活用することで、反応的軌道の並列シミュレーションが可能となり、前方および後方の速度定数を正確に計算できる。2次元ミューラー・ポテンシャルのテストケースにおいて、ブルートフォース・シミュレーションと非常に良好な一致が得られた。

ABSTRACT

A sampling procedure to compute exactly the rate of activated processes arising in systems at equilibrium or nonequilibrium steady state is presented. The procedure is a generalization of the method in [A. Warmflash, P. Bhimalapuram, and A. R. Dinner, J. Chem. Phys., 127, 154112 (2007); A. Dickson, A. Warmflash, and A. R. Dinner, J. Chem. Phys., 130, 074104 (2009)] in which one performs simulations restricted into cells by using a reinjection rule at the boundaries of the cells which is consistent with the exact probability fluxes through these boundaries. Our generalization uses results from transition path theory which indicate how to twist the dynamics to calculate reaction rates.

研究の動機と目的

• 任意の平衡または非平衡系において、統計的定常状態にある反応速度を正確に計算する一般化された手法を開発すること。
• フラックスマッチング境界条件を組み込むことで、従来の軌道並列化手法を非平衡ダイナミクスに拡張すること。
• ボロノイ細胞を通じて反応的軌道の並列計算を可能にしつつ、正確な統計的性質を保持すること。
• マイルストーン、TIS、FFSを一般化し、非単調な界面および任意の細胞配置を許容すること。
• ブルートフォース・シミュレーションと正確な比較が可能な2次元ミューラー・ポテンシャル系での手法の妥当性を検証すること。

提案手法

• 生成点の集合を用いて状態空間をボロノイ細胞に分割し、状態空間の各点が最も近い細胞に属することを保証する。
• 各細胞内で独立に反応的軌道をシミュレートし、出口点を記録して、フラックスマッチング再注入則に従い隣接細胞に再初期化する。
• 遷移経路理論を用いてダイナミクスをねじることで、反応体と生成物状態を結ぶ軌道へのサンプリングをバイアスし、正確な速度計算を可能にする。
• 非バイアスおよびねじったダイナミクスの両方で、各細胞あたり10⁸ステップをΔt = 10⁻⁴の2次精度ラングジュアン積分法でシミュレートする。
• 軌道集合から定常確率παおよびπαᴬを計算し、Aが源でBがシンクである状況での非平衡分布παᴬを表す。
• 時間反転対称性や単調な界面順序の仮定を避け、複雑で非単調な反応経路に対しても適用可能である。

実験結果

リサーチクエスチョン

• RQ1 並列化されたセルベースのシミュレーションフレームワークを用いて、非平衡系における正確な反応速度を計算できるか？
• RQ2 単調な界面順序の仮定をせず、任意のボロノイ細胞に一般化された軌道並列化は可能か？
• RQ3 遷移経路理論によるダイナミクスのねじりが、速度計算における正確な統計的性質をどの程度保持するか？
• RQ4 セルが最適でない配置にある場合や反応座標に一致しない場合でも、この手法は正確性を維持できるか？
• RQ5 計算された速度は、長時間の非バイアス直接シミュレーションで得られたものとどのように比較できるか？

主な発見

• ねじったダイナミクスを用いた反応フラックスνRは5.9×10⁻³であり、ブルートフォース値の5.8×10⁻³と非常に近い。
• 提案手法による前方速度定数kA,Bは7.4×10⁻³であり、直接シミュレーションの7.1×10⁻³と比較的一致している。
• 本手法による後方速度定数kB,Aは3.1×10⁻³であり、直接シミュレーションの3.2×10⁻³と良好に一致している。
• 本手法で計算された平衡確率παは、長時間の非バイアス軌道からの結果と良好に一致しており、図5で確認された。
• ねじったダイナミクス下での非平衡確率παᴬは、期待どおりBに向かって単調に減少し、特にπ₃₆ᴬは構成上0である。
• 非最適なセル配置に対しても本手法は頑健であり、36個のランダムに配置されたボロノイ細胞であっても正確性に影響を及ぼさなかった。