Force-based atomistic/continuum blending for multilattices

We formulate the blended force-based quasicontinuum method for multilattices and develop rigorous error estimates in terms of the approximation parameters: choice of atomistic region, blending region, and continuum finite element mesh. Balancing the approximation parameters yields a convergent atomistic/continuum multiscale method for multilattices with point defects, including a rigorous convergence rate in terms of the computational cost. The analysis is illustrated with numerical results for a Stone–Wales defect in graphene.