Some aspects of forces and fields in atomic models of crack tips

This paper examines the stresses and displacement gradients in atomistic models of cracks based on an EAM potential devised for aluminum. Methods for computing these quantities are described. Results are presented for two models differing in terms of the orientations of the crack relative to the crystal, a [100] (010) orientation that behaves in a brittle fashion and a [111] (110) orientation that emits partial dislocations prior to extending. Both models display lattice trapping. The stresses in the brittle crack model are compared with the linear elastic prediction and found to be in remarkably good agreement to within distances of about one lattice parameter of the crack tip and at the free surface where contributions from sources other than strain energy (e.g., surface tension) influence the results. Similar results are observed for the ductile model until dislocation emission occurs. The largest stresses that develop just prior to crack extension or dislocation emission are used to estimate the ratio of theoretical tensile strength to shear strength in this material. Eshelby’s conservation integrals, F and M , are also computed. F is found to be essentially contour independent and in agreement with the linear elastic prediction in both models until dislocation emission occurs, at which point a large screening contribution arises from the emitted partials. The contour size dependence of M reveals some interesting features of the crack tip including a slight wobble of the crack tip inside its potential well with changing applied K and the existence of forces acting to move the crack faces apart as blunting occurs.