Be and Be-related impurities in diamond: density functional theory study

First-principles density functional simulations were employed to investigate the geometries, electrical properties, and hyperfine structures of various beryllium-doped diamond configurations, including interstitial (Bei), substitutional (Bes), and beryllium-nitrogen (Be-N) complexes. The incorporation of Be into the diamond lattice is more favorable as a substitutional dopant than as an interstitial dopant, although both processes are endothermic. Interstitial Be could potentially exhibit motional averaging from planar to axial symmetry with an activation energy of 0.1 eV. The most stable Bes configuration has Td symmetry with a spin state of S = 1. Co-doping with nitrogen reduces the formation energy of Bes-Nn (n = 1–4) complexes, which further decreases as the number of nitrogen atoms increases. This is attributed to the smaller covalent radius of nitrogen compared to carbon, resulting in reduced lattice distortion. Bes-N3 and Bes-N4 co-doping introduces shallow donors, while Bes exhibits n-type semiconductivity, but the deep donor level renders it impractical for room-temperature applications. These findings provide valuable insights into the behavior of beryllium as a dopant in diamond and highlight the potential of beryllium-nitrogen co-doping for achieving n-type diamond semiconductors.