[論文レビュー] Analysis of a Force-Based Quasicontinuum Approximation

Matthew Dobson, Mitchell Luskin|arXiv (Cornell University)|Nov 17, 2006

Force Microscopy Techniques and Applications参考文献 17被引用数 5

ひとこと要約

本稿は、次近接原子間相互作用を伴う1次元原子格子系における力に基づく準連続体（QCF）法について、厳密な数学的解析を提供する。外部荷重に明示的な境界を設ける条件下で解の存在と一意性を確立し、ゴースト力反復法が収縮写像として全空間に収束することを証明し、実用的パラメータ領域における収束速度を定量的に評価している。

ABSTRACT

We analyze a force-based quasicontinuum approximation to a one-dimensional system of atoms that interact by a classical atomistic potential. This force-based quasicontinuum approximation is derived as the modification of an energy-based quasicontinuum approximation by the addition of nonconservative forces to correct nonphysical ``ghost'' forces that occur in the atomistic to continuum interface. We prove that the force-based quasicontinuum equations have a unique solution under suitable restrictions on the loads. For Lennard-Jones next-nearest-neighbor interactions, we show that unique solutions exist for loads in a symmetric region extending nearly to the tensile limit. We give an analysis of the convergence of the ghost force iteration method to solve the equilibrium equations for the force-based quasicontinuum approximation. We show that the ghost force iteration is a contraction and give an analysis for its convergence rate.

研究の動機と目的

• 1次元原子格子系における力に基づく準連続体（QCF）近似の平衡解の存在と一意性を確立すること。
• QCF方程式を解くために用いられるゴースト力反復法の収束挙動を分析すること。
• QCF系が適切に定式化される外部荷重の範囲を特定し、特に引張強度の限界付近での挙動を特定すること。
• ゴースト力反復法が収縮写像であることを示し、適切な条件下で全空間への収束を保証すること。
• Lennard-Jonesポテンシャルモデルにおいて、ゴースト力反復の収束速度に対して明示的な境界を導出すること。

提案手法

• 原子位置ではなく原子間距離を変数として用いることで、方程式の数学的構造を明確化するための再定式化。
• 界面に生じるアーティファクトを補正する非保存的ゴースト力を追加することで、エネルギーに基づくQCF法の変種として力に基づくQCF近似を導出する。
• 固定点理論と行列ノルム解析を用いて、外部力が所定の境界内にある場合にQCF系に一意解が存在することを証明する。
• ゴースト力反復に対して収縮写像理論を適用し、反復行列のスペクトルノルムを評価することで収束を証明する。
• ポテンシャル関数の微分の推定値を用いて収束速度を分析し、特にLennard-Jonesポテンシャルに対して評価する。
• Lennard-Jonesモデルに対して数値的に収縮条件を検証し、収縮定数が1未満となるパラメータ範囲（例：γ = 0.5）を選定する。

実験結果

リサーチクエスチョン

• RQ1 外部荷重にどのような条件を課すと、力に基づく準連続体系に一意解が存在するか？
• RQ2 QCF方程式を解くためのゴースト力反復法は、全空間に収束するか？
• RQ3 ゴースト力反復の収束速度はどのようであり、システムのパラメータにどのように依存するか？
• RQ4 QCF法は、引張強度の限界に近づいた状態でも、一意解を保つことができるか？
• RQ5 Lennard-Jonesの次近接相互作用ポテンシャルに対して、ゴースト力反復の収束を保証できるか？

主な発見

• 外部力の大きさが境界内にある場合、力に基づく準連続体方程式は一意解を持つ。1次元Lennard-Jonesモデルに対して明示的な境界が導出された。
• Lennard-Jonesの次近接相互作用を伴う系では、外部荷重が引張強度の限界に近い範囲まで拡張されても、一意解が存在する。
• ゴースト力反復法は収縮写像であり、外部荷重の境界内であれば、一意解への全空間への収束が保証される。
• ゴースト力反復の収束速度は $ \frac{8|\theta'(2r_L)|}{\theta'(r_U) - 5|\theta'(2r_L)|} $ で抑えられ、システムパラメータが適切な条件下にあればこの定数は1未満である。
• 収縮定数 $ \gamma = 0.5 $ の場合、反復は線形収束を示し、$ ||\mathbf{r}^{n+1} - \mathbf{r}||_{\infty} \leq \frac{1}{2} ||\mathbf{r}^n - \mathbf{r}||_{\infty} $ を満たす。これは高速な収束を示している。
• 解析により、ゴースト力反復法が対称的配置において、中程度の荷重下で間隔 $ (0.9706, 1.0771) $ の原子間距離の範囲で適切に定義され、収束することが確認された。