Local stress calculation in simulations of multicomponent systems

The virial and Hardy methods provide accurate local stresses for single component materials such as monatomic metals. In contrast to the elemental material case, both methods provide poor estimates of the local stress for multicomponent materials. Using binary materials such as CaO, SiC and AlN and homogeneous strain, we demonstrate that there are several sources for the slow convergence of the virial and Hardy local stresses to the bulk values. Different approaches such as enforced stoichiometry, atomic localization functions and the atomic voronoi volume are used to improve the convergence and increase the spatial resolution of the local stress. The virial method with enforced stoichiometry and atomic voronoi volumes is the most accurate, giving exact stress values by the first atomic shell. In the general case, not assuming stoichiometry, the virial method with localization functions converge to 93% of the bulk value by the third atomic shell. This work may be particularly useful for the real-time description of stresses in simulations of shock waves and deformation dynamics.