Instability of rock-salt cubic NbN in density functional calculations

All-nitride semiconductor/superconductor heterojunctions utilizing cubic niobium nitride (NbN) are a promising approach to superconducting quantum circuits for next-generation quantum-information systems. However, there are fundamental open questions on the atomic structure of NbN. In calculations with several levels of density functional theory (DFT), we find that the cubic rocksalt structure NbN (Fm-3m, 225) is energetically very unstable against the ground state hexagonal NbN in a tungsten carbide (P-6m2, 187) type lattice. To better understand the appearance of a cubic phase in numerous experiments, we perform a DFT study on possible NbN structures, determining the energy ordering between different polymorphs from databases and structure prediction. We perform supercell calculations of disordered NbN, finding that the rocksalt structure is dynamically unstable and relaxes to a lower energy monoclinic phase (C2/m, 12), which retains an approximate average cubic symmetry. However, the associated energy gain is not substantial enough for a plausible explanation for cubic NbN. We further investigate the role of external factors such as in-plane strain during epitaxial growth as well as the presence of off-stoichiometry and impurity doping on the energy ordering.