Cubic, wurtzite, and 4H-BN band structures calculated using GW methods and maximally localized Wannier functions interpolation

Quasiparticle band structures of cubic boron nitride (c-BN), wurtzite boron nitride (w-BN), and 4H-BN are calculated using GW methods. Maximally localized Wannier functions are used to interpolate quasiparticle energies at arbitrary points along high symmetry lines in the first Brillouin zone based on quasiparticle energies and wavefunctions on a uniform grid. Quasiparticle self-consistent GW (QPscGW) predictions are compared with results obtained using the one-shot G0W0 approximation and density functional theory. All of the polymorphs of BN considered in this paper have indirect band-gaps. The band-gap energies predicated using the one-shot G0W0 are 6.18 eV for c-BN, 6.86 eV for w-BN, and 7.48 eV for 4H-BN. The QPscGW calculations, which are based on the pseudopotential approximation, predict larger band gaps for these BN polymorphs.