Due to their algorithmic simplicity and high accuracy, force-based model coupling techniques are popular tools in computational physics. For example, the force-based quasicontinuum (QCF) approximation is the only known pointwise consistent quasicontinuum approximation for coupling a general atomistic model with a finite element continuum model. In this paper, we present a detailed stability and error analysis of this method. Our optimal order error estimates provide a theoretical justification for the high accuracy of the QCF approximation: they clearly demonstrate that the computational efficiency of continuum modeling can be utilized without a significant loss of accuracy if defects are captured in the atomistic region. The main challenge we need to overcome is the fact that the linearized QCF operator is typically not positive definite. Moreover, we prove that no uniform inf-sup stability condition holds for discrete versions of the W1,p-W1,q "duality pairing" with 1/p + 1/q = 1, if 1 ≤ p < ∞. However, we were able to establish an inf-sup stability condition for a discrete version of the W1,∞-W1,1 "duality pairing" which leads to optimal order error estimates in a discrete W1,∞-norm.
Bibliographical noteFunding Information:
This work was supported in part by DMS-0757355, DMS-0811039, the Department of Energy under Award Number DE-FG02-05ER25706, the Institute for Mathematics and Its Applications, the University of Minnesota Supercomputing Institute, the University of Minnesota Doctoral Dissertation Fellowship, and the EPSRC critical mass programme “New Frontier in the Mathematics of Solids”.