[論文レビュー] Coarse-graining molecular dynamics models using an extended Galerkin projection

Xiantao Li|arXiv (Cornell University)|Oct 16, 2012

Advanced Mathematical Modeling in Engineering参考文献 44被引用数 4

ひとこと要約

本稿では、結晶性固体の粗粒度化分子動力学モデルに対して、欠陥付近で粗粒度スケール部分空間を拡張することで、原子的/連続体界面におけるフォノンの反射を最小化する拡張されたガレルキン射影法を提案する。従来のガレルキン法とは異なり、一般化ランジュビン方程式における記憶関数の事前計算を回避しつつ、変分構造と力学的平衡の安定性を保持する。

ABSTRACT

We present a new framework for coarse-graining molecular dynamics models for crystalline solids. The reduction method is based on a Galerkin projection to a subspace, whose dimension is much smaller than that of the full atomistic model. The subspace is expanded by adding more coarse-grain variables near the interface between lattice defects and the surrounding regions. This effectively minimizes reflection of phonons at the interface. In this approach, there is no need to pre-compute the memory function in the generalized Langevin equations, a typical model of interface conditions. Moreover, the variational formulation preserves the stability of mechanical equilibria.

研究の動機と目的

• 粗粒度化分子動力学における原子的/連続体界面における人工的フォノン反射を解消すること。
• 一般化ランジュビン方程式における記憶核の事前計算に伴う計算負荷を克服すること。
• 粗粒度化モデルにおける元のハミルトニアン力学の変分構造と力学的平衡の安定性を保持すること。
• 結晶性固体における非平衡ダイナミクス（例：不履歴の運動や衝撃波）の正確なシミュレーションを可能にすること。
• ノイズの組み込みにより有限温度系へ適応可能なフレームワークを提供すること。

提案手法

• 基底関数で張られる次元が全原子系に比べてはるかに小さい縮小された部分空間にガレルキン射影を適用する。
• フォノン反射を抑制するために、格子欠陥および界面付近で追加の自由度を粗粒度化部分空間に追加する。
• 変分定式化を用いて力学的平衡の安定性を維持し、元のハミルトニアン力学と整合性を保つ。
• 中点台形則と有限要素離散化を用いて粗粒度方程式を近似し、コーシー＝ボーンの構成関係を有する弾性動力学モデルを導出する。
• インターフェースで重複ノードを繰り返し導入することで部分空間を拡張し、一般化ランジュビン方程式に依存せずに遷移領域を効果的にモデル化する。
• 局所的粗粒度化の改善を図るため、標準ガレルキン射影を欠陥付近で局所的リファインメントを組み込んだ拡張版に置き換える。

実験結果

リサーチクエスチョン

• RQ1 粗粒度化分子動力学モデルにおいて、原子的/連続体界面におけるフォノン反射をどのように最小化できるか？
• RQ2 一般化ランジュビン方程式における記憶関数の事前計算を回避しつつも、精度を維持できるか？
• RQ3 元の分子動力学モデルの変分構造を、粗粒度化形式でもどのように保持できるか？
• RQ4 欠陥付近での部分空間の拡張が、粗粒度化シミュレーションの精度と安定性に与える影響は何か？
• RQ5 有限温度系へ拡張可能か？その場合、揺動散逸定理を満たすか？

主な発見

• 拡張ガレルキン法は、不履歴ダイポールのシミュレーションで示されるように、従来のガレルキン法と比較して界面におけるフォノン反射を顕著に低減する。
• 従来のガレルキン法は、特に原子的領域で界面反射に起因する大きな数値誤差を生じる。
• Krylov部分空間におけるL=8の拡張ガレルキン法は、L=4よりも著しく低い誤差を達成しており、部分空間の拡張に伴い精度が向上していることが示された。
• 一般化ランジュビン方程式における記憶核の事前計算を回避でき、実装の簡素化と計算コストの低減が達成された。
• 拡張ガレルキン射影を用いて導出した粗粒度化モデルは、完全な分子動力学と一貫しており、特に不履歴コアのダイナミクスと波動伝播の捉え方において優れた結果を示した。
• 変分定式化により力学的平衡の安定性が保証され、これは多くの既存の結合手法では保証されない。