  This paper presents a consistent approach to prescribe traction boundary
conditions in atomistic models. Due to the typical multiple-neighbor
interactions, finding an appropriate boundary condition that models a desired
traction is a non-trivial task. We first present a one-dimensional example,
which demonstrates how such boundary conditions can be formulated. We further
analyze the stability, and derive its continuum limit. We also show how the
boundary conditions can be extended to higher dimensions with an application to
a dislocation dipole problem under shear stress.
