[论文解读] Analysis of a Force-Based Quasicontinuum Approximation

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

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

一句话总结

本文对一维原子系统中具有次近邻Lennard-Jones相互作用的基于力的广义连续介质（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.

研究动机与目标

• 建立一维原子系统中基于力的广义连续介质（QCF）近似平衡解的存在性与唯一性。
• 分析求解QCF方程所用鬼力迭代方法的收敛行为。
• 量化外部载荷的范围，使QCF系统在接近断裂极限时仍保持适定性。
• 证明鬼力迭代为压缩映射，从而在适当条件下确保全局收敛。
• 为Lennard-Jones势模型提供鬼力迭代收敛速率的显式边界。

提出的方法

• 将平衡方程重新表述为原子间距而非原子位置的形式，以明确其数学结构。
• 通过添加非保守的鬼力以校正界面伪影，将基于能量的QCF方法修改为基于力的QCF近似。
• 利用不动点理论与矩阵范数分析，证明当外部力在指定范围内时，QCF系统具有唯一解。
• 将压缩映射理论应用于鬼力迭代，通过有界迭代矩阵的谱范数证明收敛性。
• 利用势函数导数的估计分析收敛速率，特别针对Lennard-Jones势。
• 通过选择收缩常数小于1的参数范围（例如 γ = 0.5），对Lennard-Jones模型数值验证收缩条件。

实验结果

研究问题

• RQ1 在何种外部载荷条件下，基于力的广义连续介质系统具有唯一解？
• RQ2 求解QCF方程的鬼力迭代方法是否具有全局收敛性？
• RQ3 鬼力迭代的收敛速率如何？其与系统参数的关系为何？
• RQ4 在接近断裂极限时，QCF方法在多大程度上仍能保持唯一解？
• RQ5 对于Lennard-Jones次近邻相互作用势，能否保证鬼力迭代的收敛性？

主要发现

• 当外部力的大小受约束时，基于力的广义连续介质方程具有唯一解，且为一维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) $ 内定义良好且收敛。