Molecular Dynamics Simulations of Vacancy Generation and Migration near a Monocrystalline Silicon Surface during Energetic Cluster Ion Implantation

The process of ion implantation often involves vacancy generation and migration. The vacancy generation and migration near a monocrystalline silicon surface during three kinds of energetic Si35 cluster ion implantations were investigated by molecular dynamics simulations in the present work. The patterns of vacancy generation and migration, as well as the implantation-induced amorphous structure, were analyzed according to radial distribution function, Wigner–Seitz cell, and identify diamond structure analytical methods. A lot of vacancies rapidly generate and migrate in primary directions and form an amorphous structure in the first two picoseconds. The cluster with higher incident kinetic energy can induce the generation and migration of more vacancies and a deeper amorphous structure. Moreover, boundaries have a loading–unloading effect, where interstitial atoms load into the boundary, which then acts as a source, emitting interstitial atoms to the target and inducing the generation of vacancies again. These results provide more insight into doping silicon via ion implantation.