Defect dynamics in two-dimensional black phosphorus under argon ion irradiation.

Fundamental understanding of the atomic-scale mechanisms underlying production, accumulation, and temporal evolution of defects in phosphorene during noble-gas ion irradiation is crucial to design efficient defect engineering routes to fabricate next-generation materials for energy technologies. Here, we employed classical molecular dynamics (CMD) simulations using a reactive force field to unravel the effect of defect dynamics on the structural changes in a monolayer of phosphorene induced by argon-ion irradiation, and its subsequent relaxation during post-radiation annealing treatment. Analysis of our CMD trajectories using unsupervised machine learning methods showed that radiation fluence strongly influences the types of defect that form, their dynamics, and their relaxation mechanisms during subsequent annealing. Low ion fluences yielded a largely crystalline sheet featuring isolated small voids (up to 2 nm), Stone-Wales defects, and mono-/di-vacancies; while large nanopores (∼10 nm) can form beyond a critical fluence of ∼10