[論文レビュー] (In-)Stability and Stabilisation of QNL-Type Atomistic-to-Continuum Coupling Methods

Christoph Ortner, Alexander V. Shapeev+1 more|arXiv (Cornell University)|Aug 18, 2013

Numerical methods for differential equations参考文献 19被引用数 6

ひとこと要約

本稿は、1次元および2次元におけるQNL型原子論的から連続体への（a/c）カップリング法の包括的な安定性理論を確立する。1次元では普遍的に安定なa/cカップリングスキームを導入し、不安定なスキームの安定化機構を提示する一方で、2次元では分岐点付近での一貫性を損なう、新たな深刻な不安定性を特定し、既存のゴースト力フリー法に根本的な限界を露呈する。

ABSTRACT

We study the stability of ghost force-free energy-based atomistic-to-continuum coupling methods. In 1D we essentially complete the theory by introducing a universally stable a/c coupling as well as a stabilisation mechanism for unstable coupling schemes. We then present a comprehensive study of a two-dimensional scalar planar interface setting, as a step towards a general 2D/3D vectorial analysis. Our results point out various new challenges. For example, we find that none of the ghost force-free methods known to us are universally stable (i.e., stable under general interaction and general loads). We then explore to what extent our 1D stabilisation mechanism can be extended.

研究の動機と目的

• a/cカップリング法における長年の未解決問題である安定性、特に保存的QNL型スキームに関する問題を解消すること。
• 一般の相互作用および外力下でも安定を保つ、1次元における普遍的安定a/cカップリングスキームを構築すること。
• 1次元とは異なり、2次元におけるQNL型a/cカップリングに内在する、新たな本質的不安定性を同定・分析すること。
• 1次元の安定化機構が2次元にどの程度拡張可能かを調査し、特に分岐点付近における一貫性損失を定量化すること。
• 1次元を超えるa/cカップリング法における安定性と一貫性のトレードオフの理論的基盤を提供すること。

提案手法

• 有限差分作用素と反射対称性を持つ多体サイトエネルギー関数を用いて1次元QNL定式化を構築する。
• 一般条件の下で安定定数が正となるように構成された、普遍的に安定な新しいa/cカップリングスキームを導入する。
• コンcentration compactnessの議論を適用して、QNLエネルギー汎関数のヘッセ行列の極限を分析し、安定定数が欠損部とQNL安定定数の最小値に収束することを証明する。
• 変位場を局所化された成分と減衰する成分に分解することで、大系におけるヘッセ行列の下界を推定する。
• ポテンシャルの2階微分のリプシッツ連続性を用いて、エネルギーヘッセ行列における摂動を制御する。
• 変形およびポテンシャルに関する厳密な仮定の下で、2次元スカラー平面界面問題を分析し、1次元には存在しない新たな不安定性の原因を同定する。

実験結果

リサーチクエスチョン

• RQ1 一般の相互作用ポテンシャルおよび荷重に依存しない、1次元における普遍的安定なQNL型a/cカップリングスキームを構築可能か？
• RQ2 2次元におけるQNL型a/cカップリングの不安定性の性質と起源は何か？1次元の不安定性とはどのように異なるか？
• RQ3 1次元の安定化機構は2次元にどの程度まで拡張可能か？その結果、一貫性にどのような影響が生じるか？
• RQ4 2次元における不安定性はQNL型スキームに特有のものか、より広い一般化幾何再構成法のクラスにも影響を及ぼすか？
• RQ5 2次元における安定化は、特に分岐点付近において、その方法の一貫性にどのような影響を与えるか？

主な発見

• 一般の相互作用ポテンシャルおよび外部荷重下でも安定を保証する、新たな普遍的安定1次元a/cカップリングスキームが構築された。
• 2次元スカラー平面界面モデルは、1次元に存在しない新たな不安定性要因を明らかにしたが、これは本質的であり、単純な修正では除去できない。
• この2次元不安定性は普遍的であり、QNL型スキームに限らず、広範な一般化幾何再構成法のクラスに影響を及ぼす。
• 2次元におけるQNL型スキームの安定化は、系が分岐点に近づくと、一貫性を著しく損なう。
• 1次元におけるヘッセ行列の安定定数の極限が、欠損部とQNL安定定数の最小値に収束することを示し、安定性解析の厳密な基礎を提供した。
• 分析により、既存のゴースト力フリー法が2次元では普遍的に安定ではないことが確認され、現在のa/cカップリング理論における重要なギャップが浮き彫りになった。