[论文解读] The Spectrum of the Force-Based Quasicontinuum Operator for a Homogeneous Periodic Chain

Matthew Dobson, Christoph Ortner+1 more|Warwick Research Archive Portal (University of Warwick)|Apr 20, 2010

Numerical methods in engineering参考文献 16被引用 7

一句话总结

本文严格证明了在第二近邻相互作用情况下，均匀周期性链的线性化力基准连续（QCF）算子具有正谱，且其谱与准非局部（QNL）算子完全相同，并建立了其特征向量矩阵条件数在系统尺寸和原子区域尺寸下的一致有界性。这些结果验证了先前的数值猜想，并为GMRES求解器的收敛性提供了严格估计，同时为QCF方法提供了新的稳定性界。

ABSTRACT

We show under general conditions that the linearized force-based quasicontinuum (QCF) operator has a positive spectrum, which is identical to the spectrum of the quasinonlocal quasicontinuum (QNL) operator in the case of second-neighbour interactions. Moreover, we establish a bound on the condition number of a matrix of eigenvectors that is uniform in the number of atoms and the size of the atomistic region. These results establish the validity of and improve upon recent conjectures ([arXiv:0907.3861, Conjecture 2] and [arXiv:0910.2013, Conjecture 8]) which were based on numerical experiments. As immediate consequences of our results we obtain rigorous estimates for convergence rates of (preconditioned) GMRES algorithms, as well as a new stability estimate for the QCF method.

研究动机与目标

• 严格验证关于线性化力基准连续（QCF）算子谱和特征向量条件数的数值猜想。
• 在离散Sobolev型范数下，建立QCF方法在系统尺寸和原子区域尺寸上均独立的统一稳定性界。
• 为应用于QCF平衡方程的迭代求解器（如GMRES）的收敛性提供理论基础。
• 将谱分析从最近邻相互作用扩展到一般有限范围相互作用的周期链。
• 在$υ^{2,2}$范数下，推导出QCF方法的新稳定性估计，其独立于原子数$N$和原子区域$τ$。

提出的方法

• 分析基于周期链中$N$个原子在均匀形变状态下的QCF算子线性化。
• 通过谱理论和在$×^N$中的矩阵表示，将QCF算子的谱与基于能量的准非局部（QNL）算子的谱进行比较。
• 显式构造一个特征向量矩阵，并通过谱分解和算子范数估计，统一地在$N$和原子区域尺寸上界定了其条件数。
• 对于有限范围相互作用，通过势能二阶导数的界，证明QCF算子的谱位于原子和连续算子谱之间。
• 利用统一的条件数界和预处理算子的谱性质，推导出预处理GMRES的收敛速率。
• 通过在均值为零的离散Sobolev空间中分析QCF算子的逆，获得$υ^{2,2}$范数下的新稳定性估计。

实验结果

研究问题

• RQ1 在第二近邻相互作用情况下，线性化QCF算子的谱是否与QNL算子的谱一致？
• RQ2 QCF算子特征向量矩阵的条件数是否在原子数和原子区域尺寸上一致有界？
• RQ3 QCF算子的谱性质能否从最近邻相互作用推广到一般有限范围相互作用？
• RQ4 求解QCF平衡方程的预处理GMRES方法的收敛速率如何？
• RQ5 QCF算子在离散$υ^{2,2}$范数下是否稳定，且其稳定性是否在系统尺寸和原子区域尺寸上一致？

主要发现

• 在第二近邻相互作用情况下，线性化QCF算子的谱为正，且与QNL算子的谱完全相同。
• QCF算子特征向量矩阵的条件数在$N$和原子区域尺寸上一致有界，优于先前的猜想。
• 对于有限范围相互作用，QCF算子的谱位于原子和连续算子谱之间，且界不依赖于$N$和$τ$。
• 预处理GMRES方法求解QCF系统时，残差满足界$\|L_{1}^{-1}r^{(m)}\| \leq C N^3 q^m \|L_{1}^{-1}r^{(0)}\|$和$\|L_{1}^{-1/2}r^{(m)}\| \leq C N^2 q^m \|L_{1}^{-1/2}r^{(0)}\|$，其中$C>0$，$q\in(0,1)$，且均不依赖于$N$和$\mathcal{A}$。
• QCF算子在$\mathcal{U}^{2,2}$范数下稳定：$\|(L^{\rm qcf}_0)^{-1}\|_{L(\mathcal{U}^{0,2},\mathcal{U}^{2,2})}$有界，其界仅依赖于势参数，而不依赖于$N$或$\mathcal{A}$。