Automatic adaptivity in the fully nonlocal quasicontinuum method for coarse-grained atomistic simulations

Abstract

The quasicontinuum (QC) method is a concurrent scale-bridging technique that extends atomistic accuracy to significantly larger length scales by reducing the full atomic ensemble to a small set of representative atoms and using interpolation to recover the motion of all lattice sites where full atomistic resolution is not necessary. While traditional QC methods thereby create interfaces between fully resolved and coarse-grained regions, the recently introduced fully nonlocal QC framework does not fundamentally differentiate between atomistic and coarsened domains. Adding adaptive refinement enables us to tie atomistic resolution to evolving regions of interest such as moving defects. However, model adaptivity is challenging because large particle motion is described based on a reference mesh (even in the atomistic regions). Unlike in the context of, for example, finite element meshes, adaptivity here requires that (i) all vertices lie on a discrete point set (the atomic lattice), (ii) model refinement is performed locally and provides sufficient mesh quality, and (iii) Verlet neighborhood updates in the atomistic domain are performed against a Lagrangian mesh. With the suite of adaptivity tools outlined here, the nonlocal QC method is shown to bridge across scales from atomistics to the continuum in a truly seamless fashion, as illustrated for nanoindentation and void growth.