M. Dobson and M. Luskin
SIAM Journal of Numerical Analysis, 47, 2455–2475 (2009).
We derive a model problem for quasicontinuum approximations that allows a simple, yet insightful, analysis of the optimal-order convergence rate in the continuum limit for both the energy-based quasicontinuum approximation and the quasi-nonlocal quasicontinuum approximation. For simplicity, the analysis is restricted to the case of second-neighbor interactions and is linearized about a uniformly stretched reference lattice. The optimal-order error estimates for the quasi-nonlocal quasicontinuum approximation are given for all strains up to the continuum limit strain for fracture. The analysis is based on an explicit treatment of the coupling error at the atomistic-to-continuum interface, combined with an analysis of the error due to the atomistic and continuum schemes using the stability of the quasicontinuum approximation.