Abstract: Using self-consistent KohnSham density-functional theory molecular dynamics simulations, we demonstrate the theoretical possibility to synthesize NiC60, the incarfullerene Ni@C60 and the heterofullerene C59Ni in an ion implantation setup. The corresponding formation mechanisms of all three complexes are elucidated as a function of the ion implantation energy and impact location, suggesting possible routes for selectively synthesizing these complexes.