Molecular dynamics investigation of dynamic crack stability

A series of molecular-dynamics simulations has been performed in order to evaluate the effects of several physical factors on dynamic crack stability. These factors are the crystalline structure and the interatomic interaction modeled by various empirical potentials. For brittle crack propagation at low temperature we find that steady-state crack velocities are limited to a band of accessible values. Increasing the overload beyond K-lc, the crack can propagate with a steady-state velocity, which quickly reaches the terminal velocity of about 0.4 of the Rayleigh wave speed. The magnitude of the terminal velocity can be related to the nonlinearity of the interatomic interaction. Further increasing the overload does not change the steady-state velocity dramatically, but significantly increases the amplitude of acoustic emission from the crack tip. Loading the crack even further leads to instabilities which take the form of cleavage steps, dislocation emission, or branching. The instability is closely related to the buildup of a localized coherent, phononlike field generated by the bond-breaking events. The form of the instability depends critically on crystal structure and on the crystallographic orientation of the crack system but can also be correlated with the relative ease of dislocation generation (and motion).