▎ 摘　　要

The fracture behavior of trilayer graphene was investigated using molecular dynamics simulations, with the emphasis on the effect of vacancies and interlayer bonds. While both defects make the fracture happen earlier and reduce the ultimate strength, the exact details of how they change the fracture behavior are different. Interlayer bonds can make the three graphene layers fracture in a more simultaneous way by transferring and distributing load more uniformly between these layers when loaded in the zigzag direction. When loaded in the armchair direction, however, the interlayer bonds can shift the fracture plane from a bifurcating mode to a zigzag direction-dominated mode. In contrast, with increasing tensile strain, vacancies can change the fracturing manner from a catastrophic mode to a successive one, leading to a stepwise stress strain curve and the formation of carbon atom chains and polygonal rings during the fracture process. The mixture of vacancies and interlayer bonds breaks the trilayer graphene in a way coupling the effect of both defects, depending on the chirality and defect density. The results presented in this work can provide a guideline for the design and application of multilayer graphene.