[論文レビュー] Energy-based ghost force removing techniques for the quasicontinuum method

Ping Lin, Alexander V. Shapeev|ArXiv.org|Sep 29, 2009

Force Microscopy Techniques and Applications参考文献 24被引用数 12

ひとこと要約

本稿では、準連続体（QC）シミュレーションにおける局所的・非局所的界面におけるゴースト力の除去を、有限要素法に類似したノードフレームワークに原子的自由度を射影することによって行うエネルギーに基づく手法、すなわち準連続体射影（QCP）法を提案する。QCP法は、QNL や GCR と同等またはそれ以上の精度を達成し、誤差は同等または低く、計算コストも同じである。また、標準的な有限要素理論に依存し、自由パラメータが存在しないため、実装や解析がより簡単である。

ABSTRACT

This paper studies numerical methods for accurate treatment of the interface between the local and the nonlocal region in a QC approximation of atomistic materials. Only the energy-based methods are considered. Particularly, a quasicontinuum projection (QCP) method based on the idea of finite elements is shown to be accurate and efficient for this problem. We analyse the QCP method and study its relation to the existing methods, such as the quasinonlocal quasicontinuum method and the geometrically consistent reconstruction-based method. The analysis and the results of numerical tests confirm that the projection-based QC method successfully removes the ghost force with the same computational cost as the other methods. In all computed examples the error of QCP is either the same or lower as the error of the other methods. The performance of these methods in treating interfaces of elements in the local region is also examined.

研究の動機と目的

• 準連続体（QC）法における局所的・非局所的界面に生じるゴースト力の持続的問題に取り組み、アトムスティックから連続体へのマルチスケールシミュレーションにおける精度を損なう要因を解消すること。
• ゴースト力の除去を実現するエネルギーに基づく手法を開発し、余分な力やエネルギー保存則の違反を引き起こさないようにすること。
• 特に、Cauchy-Born 外挿法に起因する界面誤差が顕著に現れる2次元（2D）ケースにおいて、QC シミュレーションの精度を向上させること。
• GCR や他の既存手法と比較して、問題固有のパラメータを必要とせず、標準的な有限要素概念に依存するため、実装や解析がより簡単である手法を提供すること。
• QCP 法が、QNL や GCR といった最先端の手法に匹敵する収束速度と誤差水準を達成できることを示すこと。

提案手法

• QCP 法は、有限要素法に類似したフレームワークを用いて、原子的自由度をノード点に射影することで、一貫性のあるエネルギーに基づく定式化を可能にする。
• 非ノード原子の位置をノード原子の位置の線形結合として表現することで、系の全ポテンシャルエネルギーを定式化し、エネルギーの一貫性を確保する。
• 射影されたエネルギーに変分原理を適用して力および剛性行列を導出し、界面におけるゴースト力の発生を回避する。
• 局所領域における Cauchy-Born 外挿法を回避することで、特に遅く減衰する（1/r のように）摂動が生じる2次元領域で、端面および界面誤差を低減する。
• QCPm 変種では、効率性を高めるために局所領域で Cauchy-Born 外挿法を許容しつつも、界面の精度を維持する。
• この手法は有限要素理論に数学的に根ざしており、近似誤差理論を用いた収束解析が可能である。

実験結果

リサーチクエスチョン

• RQ1 エネルギーに基づく手法は、局所的・非局所的界面におけるゴースト力を、余分な力の発生なしに効果的に除去できるか？
• RQ2 1次元および2次元のテストケースにおいて、QCP 法は QNL や GCR 法と比較して、精度と計算コストの面でどのように異なるか？
• RQ3 局所領域における Cauchy-Born 外挿法の回避が、特に2次元において顕著に界面誤差を低減するのか？
• RQ4 自由パラメータがなく、標準的な有限要素概念に依存するため、GCR よりも実装や解析が簡単であるか？
• RQ5 QCP 法は、金属合金のような非一様格子を持つ複雑な材料へ一般化可能か？

主な発見

• QCP 法は、QCE や QNL、GCR と同等の計算コストで、局所的・非局所的界面におけるゴースト力を効果的に除去した。
• すべての数値例において、QCP 法の誤差は QCE や QNL と同等またはそれ以下であり、2次元ケースでは GCR と同等またはそれ以下であった。
• QCP 法は、自由度（DoF）の増加に伴い、誤差が安定して減少する傾向を示した。これに対して、QCE や QNL はゼロへの収束に失敗した。
• 2次元シミュレーションでは、Cauchy-Born 外挿法を用いる手法に比べ、QCP 法が端面および界面誤差が低減したため、優れた性能を示した。特に、摂動がゆっくり減衰する（1/r のように）領域で顕著であった。
• QCP 法は、特定の GCR 法のインスタンスと数学的に同等であるが、事前計算された再構成係数や界面距離の計算を必要としないため、実装や解析がより簡単である。
• QCPm 変種は、局所領域で Cauchy-Born 外挿法を用いるが、GCR とほぼ同一の結果を達成する一方で、実装性と一般化可能性が向上している。