Molecular dynamics study on atomic elastic stiffness at mode I crack tip in Si: Precursor instability in their eigenvalue before crack propagation

Abstract

We performed molecular dynamics simulations on the [001](010), [001](110) and [1¯1¯2](111) mode I through cracks in Si, and discussed the unstable crack propagation based on the determinant and eigenvalue of the atomic elastic stiffness (AES), Bijα, under the finite temperature of 300K. Two significant findings are obtained; the one is the existence of extremely high detBijα atoms, or “highly stable” zone, at the crack tip although the tip surface actually has negative detBijα atoms. The other finding is the precursor instability of the 1st eigenvalue of [Bijα]{Δεj}=ηα{Δεj}. Although there is a little difference in the responses of the [001](010), [001](110) and [1¯1¯2](111) cracks, extremely negative eigenvalue of ηα(1) is observed before the unstable crack propagation. The corresponding eigenvector shows large positive Δε2 (crack opening or mode I) in the [001](010) crack, large negative Δε1 (intrusion in the x-direction or mode II) in the [001](110) crack and large Δε5 (shear in the zx or mode III) in the [1¯1¯2](111) crack. However, the fracture process in the later stage is rather complicated so that the final fracture mode is not correspond to these instability mode.