▎ 摘　　要

The evolution of multiple vacancies (V(n)s) in graphene under electron irradiation (EI) was explored systematically by long time non-equilibrium molecular dynamics simulations, with n varying from 4 to 40. The simulations showed that the V(n)s form haeckelites in the case with small n, while forming holes as n increases. The scale of the haeckelites, characterized by the number of pentagon-heptagon pairs, grows linearly with n. Such a linear relationship can be interpreted as a consequence of compensating the missing area, caused by the V(n)s, in order to maintain the area of the perfect sp(2) network by self-healing. Beyond that, the scale of the haeckelite vs. the density of missing atoms is predicted to be S-h similar to 6D(n), where S-h and D-n are the percentage of non-hexagonal rings and missing atoms, respectively. This study provides an intuitive picture of the formation of amorphous graphene under EI and the quantitative understanding of the mechanism.