Diffusion of point defects in crystalline silicon using the kinetic activation-relaxation technique method

We study point-defect diffusion in crystalline silicon using the kinetic activation-relaxation technique (k-ART), an off-lattice kinetic Monte Carlo method with on-the-fly catalog building capabilities based on the activation-relaxation technique (ART nouveau), coupled to the standard Stillinger-Weber potential. We focus more particularly on the evolution of crystalline cells with one to four vacancies and one to four interstitials in order to provide a detailed picture of both the atomistic diffusion mechanisms and overall kinetics. We show formation energies, activation barriers for the ground state of all eight systems, and migration barriers for those systems that diffuse. Additionally, we characterize diffusion paths and special configurations such as dumbbell complex, di-interstitial (IV-pair+2I) superdiffuser, tetrahedral vacancy complex, and more. This study points to an unsuspected dynamical richness even for this apparently simple system that can only be uncovered by exhaustive and systematic approaches such as the kinetic activation-relaxation technique.