[論文レビュー] Pseudo generators of spatial transfer operators

Andreas Bittracher, Péter Koltai+1 more|arXiv (Cornell University)|Dec 4, 2014

Quantum chaos and dynamical systems参考文献 35被引用数 4

ひとこと要約

本稿では、半群でないが、運動量平均を用いず、4次までの明示的で well-defined なテイラー展開を許容する、分子動力学における空間的転送作用素の疑似生成子を導入する。この手法により、軌道ベースの手法と比較して計算コストを著しく削減しつつ、代わりに代用的動的解析の精度を保ちながら、メタ安定的ダイナミクスの効率的なコロケーションベース計算が可能になる。

ABSTRACT

Metastable behavior in dynamical systems may be a significant challenge for a simulation based analysis. In recent years, transfer operator based approaches to problems exhibiting metastability have matured. In order to make these approaches computationally feasible for larger systems, various reduction techniques have been proposed: For example, Schütte introduced a spatial transfer operator which acts on densities on configuration space, while Weber proposed to avoid trajectory simulation (like Froyland et al.) by considering a discrete generator. In this manuscript, we show that even though the family of spatial transfer operators is not a semigroup, it possesses a well defined generating structure. What is more, the pseudo generators up to order 4 in the Taylor expansion of this family have particularly simple, explicit expressions involving no momentum averaging. This makes collocation methods particularly easy to implement and computationally efficient, which in turn may open the door for further efficiency improvements in, e.g., the computational treatment of conformation dynamics. We experimentally verify the predicted properties of these pseudo generators by means of two academic examples.

研究の動機と目的

• 生物分子における希少なコンformational 変化遷移をシミュレートすることが計算的に非現実的であるという問題に取り組むこと。
• 大規模系において、高価な軌道シミュレーションを要する転送作用素手法の制限を克服すること。
• 半群でない空間的転送作用素に対して、生成子に基づくアプローチを構築し、効率的な数値近似を可能にすること。
• 運動量平均を避けることにより、計算に適した明示的表現を、4次までの疑似生成子に対して得ること。
• 高次元系におけるメタ安定的ダイナミクスの効率的かつ正確な近似を可能にするコロケーション法を実現すること。

提案手法

• 半群でないが、well-defined な生成構造を有する空間的転送作用素の族に対する疑似生成子を定義する。
• 位置空間における拡散係数とドリフト係数を用いて、空間的転送作用素の4次までの明示的テイラー展開係数を導出する。
• 小さな時間における空間的転送作用素を近似する「修復済み」作用素を構築し、元のダイナミクスと整合性を保つ。
• 疑似生成子に基づくメッシュフリーなコロケーション法を用い、軌道統合の必要性を排除する。
• 拡散行列の行列式の一様な下界（W₃）を証明することで、固有関数の滑らかさを保証し、遷移密度の正則性を保証する。
• ドリフト係数と拡散係数の周期的拡張を用いて、トーラス上での周期的系へと結果を拡張する。

実験結果

リサーチクエスチョン

• RQ1 半群でない空間的転送作用素に、標準的な無限小生成子を持たないにもかかわらず、生成子に基づくアプローチを適用できるか？
• RQ2 空間的転送作用素の疑似生成子に対して、4次までの明示的で計算に適した表現をどのように導出できるか？
• RQ3 これらの疑似生成子の構築において、運動量平均を避けるにはどうすればよいか？計算効率の向上に寄与する。
• RQ4 疑似生成子は、元の空間的転送作用素のメタ安定的構造と主要な時間スケールをどの程度正確に保持するか？
• RQ5 トーラス上での分子動力学のような周期的系において、空間的転送作用素の固有関数の滑らかさを厳密に確立できるか？

主な発見

• 半群でないにもかかわらず、空間的転送作用素の族は well-defined な生成構造を有し、疑似生成子の定義を可能にする。
• 4次までの疑似生成子は、運動量平均を必要とせず、明示的で閉形式の表現を持つ。実装が簡素化される。
• 本手法により、数値的軌道統合を必要としないコロケーションベースの近似が可能となり、計算コストが著しく削減される。
• ポテンシャルおよび拡散係数にややのRegularな条件下で、すべての t > 0 に対して空間的転送作用素の固有関数は滑らかである。
• トーラス上での周期的系では、遷移密度および固有関数が滑らかであり、導関数が指数的に減衰する。これにより、コロケーションスキームの適切な定義と収束性が保証される。
• 数値実験により、疑似生成子アプローチの予測された精度と効率性が、学術的テストケースで確認された。