Linear Stationary Iterative Methods for the Force-based Quasicontinuum Approximation

Force-based multiphysics coupling methods have become popular since they provide a simple and efficient coupling mechanism, avoiding the difficulties in formulating and implementing a consistent coupling energy. They are also the only known pointwise consistent methods for coupling a general atomistic model to a finite element continuum model. However, the development of efficient and reliable iterative solution methods for the force-based approximation presents a challenge due to the non-symmetric and indefinite structure of the linearized force-based quasicontinuum approximation, as well as to its unusual stability properties. In this paper, we present rigorous numerical analysis and computational experiments to systematically study the stability and convergence rate for a variety of linear stationary iterative methods.