Statistical investigations on nitrogen-vacancy center creation

Experiments show that shallow nitrogen implantations (<10keV) result in a negatively charged nitrogen-vacancy center (NV-) yield of 0.01-0.1%. The most succesful technique for introducing NV- centers in the carbon matrix is ion implantation followed by annealing at 1100K. We investigated the influence of channeling effects during shallow implantation and statistical diffusion using molecular dynamics (MD) and Monte Carlo (MC) approaches. Energy barriers for the diffusion process were calculated using the density functional theory (DFT). Our simulations show a significant difference in the NV yield compared to the experiment. Statistically, 25% of the implanted nitrogens form a NV center after annealing.