Atomistic simulation study on the crack growth stability of graphene under uniaxial tension and indentation

Abstract

Combining a series of atomistic simulations with fracture mechanics theory, we systematically investigate the crack growth stability of graphene under tension and indentation, with a pre-existing crack made by two methods: atom removal and (artificial) bonding removal. In the tension, the monotonically increasing energy release rate umentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G$$\end{document} is consistent with the unstable crack growth. In contrast, the non-monotonic G\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G$$\end{document} with a maximum for indentation explains the transition from unstable to stable crack growth when the crack length is comparable to the diameter of the contact zone. We also find that the crack growth stability within a stable crack growth regime can be significantly affected by the crack tip sharpness even down to a single atom scale. A crack made by atom removal starts to grow at a higher indentation force than the ultimately sharp crack made by bonding removal, which leads to a large force drop at the onset of the crack growth that can cause unstable crack growth under indentation with force control. In addition, we investigate the effect of the offset distance between the indenter and the crack to the indentation fracture force and find that the graphene with a smaller initial crack is more sensitive. The findings reported in this study can be applied to other related 2D materials because crack growth stability is determined primarily by the geometrical factors of the mechanical loading.