[論文レビュー] Nonlinear Elasticity from Atomistic Mechanics

Christoph Ortner, Florian Theil|arXiv (Cornell University)|Feb 17, 2012

Advanced Mathematical Modeling in Engineering参考文献 21被引用数 7

ひとこと要約

この論文は、静的および動的問題の両方における原子論的から連続体への結合を分析することで、非線形弾性におけるコーシー＝ボーン近似の厳密な収束限界を確立する。格子安定性および速やかに減衰する相互作用の仮定の下で、小規模なデータ条件下で、原子論的解がコーシー＝ボーン解に、変位勾配およびエネルギーの両方において $ε^2$ の誤差で収束することを証明する。これは、応力の局所的整合性と局所化技術を用いることで達成される。

ABSTRACT

We present sharp convergence results for the Cauchy--Born approximation of general classical atomistic interactions, for both static and dynamic problems, for small data.

研究の動機と目的

• 非線形弾性における一般の多体原子論的相互作用に対するコーシー＝ボーン近似を厳密に正当化すること。
• 静的問題から動的問題（時間発展方程式）への収束結果の拡張により、理論的理解のギャップを埋めること。
• 最小限の物理的仮定の下で、静的および動的弾性の両方における鋭い誤差推定を確立すること。
• 格子安定性が収束のための必要条件であることを、不安定性が無限大に発散する解を生じる反例を用いて示すこと。

提案手法

• Bravais格子 $\varepsilon\Lambda$ 上でスケーリングされた多体ポテンシャル $V$ を用いて原子論的エネルギー $E^\varepsilon$ を定式化し、$\Lambda = \mathbb{Z}^d$ とする。
• 連続体のコーシー＝ボーンエネルギー密度 $W(\mathbf{F}) = V(\{\mathbf{F}\rho\}_{\rho \in \Lambda \setminus \{0\}})$ を定義し、離散的および連続的運動学の間の橋渡しを図る。
• 離散的対象を連続的に表現するための局所化技術を用い、弱い定式化と応力の比較を可能にする。
• 原子論的応力とコーシー＝ボーンモデルにおけるパイオラ＝キルホフ応力の間で、局所的2次整合性を確立する。
• エネルギー推定とスペクトル解析を用いて、$L^2$ および $H^1$ ノルムにおける原子論的解とコーシー＝ボーン解の差を評価する。
• 近似固有関数を用いて、格子不安定性が指数関数的成長を引き起こすことを示し、安定性仮定の必要性を証明する。

実験結果

リサーチクエスチョン

• RQ1 静的非線形弾性問題において、コーシー＝ボーン近似が原子論的解に収束する条件は何か？
• RQ2 ニュートンの運動方程式に従う原子論的系の動的挙動を、コーシー＝ボーンモデルが正確に近似できるか？
• RQ3 小スケール格子間隔の極限 $\varepsilon \to 0$ において、原子論的解とコーシー＝ボーン解の間の収束速度は正確にどの程度か？
• RQ4 格子安定性は収束のための必要条件であるのか、それともより弱い仮定のもとで収束が成立するのか？
• RQ5 多体ポテンシャルの構造（例：二体、三体、四体）は、コーシー＝ボーン近似の整合性および収束にどのように影響するか？

主な発見

• 静的問題において、外部力および初期データが小さい条件下で、原子論的解とコーシー＝ボーン最小化子の勾配差は $C\varepsilon^2$ で有界である。
• 動的問題において、初期データの摂動が小さい条件下で、時間区間 $[0,T]$ において勾配の $L^2$ ノルムにおける解の差は $C\varepsilon^2$ で有界である。この $T$ は $\varepsilon$ に依存しない。
• 収束速度 $\varepsilon^2$ は鋭く、パイオラ＝キルホフ応力テンソルの局所的2次整合性に起因する。
• 格子安定性は収束のための十分条件であるだけでなく、必要条件でもある。安定性がないと、原子論的系は指数関数的成長を示す。
• 不安定な格子（例：$a_1 = -1, a_2 = 1/2$）を用いた反例では、原子論的エネルギーが下界を持たず、局所的最小解が無効になり、動的不安定性が生じることが示された。
• 本分析は、置換不変かつ安定性および減衰条件を満たす限り、一般の多体ポテンシャル（二体、三体、四体、二面角項、ねじれ項も含む）をカバーする。