Bounds on the range of density-functional-theory point-defect levels in semiconductors and insulators

• Published in: Computational materials science 2014-09, Vol.92 (C), p.431-438
• Main Authors: Modine, N.A., Wright, A.F., Lee, S.R.
• Format: Article
• Language: English
• Subjects: Charge; Condensed matter: electronic structure, electrical, magnetic, and optical properties; Defects; Density functional theory; Electron states; Energy use; Exact sciences and technology; Functionals; Impurity and defect levels; Insulators; Mathematical analysis; Physics; Semiconductors; Thunderstorms
• ISSN: 0927-0256; 1879-0801
• DOI: 10.1016/j.commatsci.2014.05.032

[Display omitted] •Derive supercell-size dependent bounds on DFT calculated defect levels.•Identify when defect levels involve delocalized band-edge-like states.•Show that finite supercells help mitigate a overly small Kohn–Sham band gap. Defects in semiconductors and insulators are characterized by their levels, which are defined as the values of the Fermi level at which the charge state of the defect changes. Kohn–Sham density functional theory calculations for charged defects have been widely and successfully used to predict defect levels. Due to their lower computational cost and demonstrated ability to predict levels spanning the measured band gap, semilocal exchange–correlation functionals are widely used in these calculations. However, there is a potential pitfall in using semilocal functionals: although they often predict accurate energies for adding or removing electrons from states that are localized near the defect, the famous band gap error results in overly small energies for adding or removing electrons from extended band edge states. As a result, electrons (or holes) that should occupy localized states may become partially or fully delocalized. In order to help detect and analyze such cases, we introduce bounds on the defect levels that can be obtained using a given functional, supercell, and Brillouin zone sampling. Since these bounds correspond to the charge transition levels of the corresponding defect-free supercell, comparison with the bounds reveals when a calculated level is behaving in a bulk-like rather than defect-like manner. We find that the bounds depend significantly on supercell size due to band-filling effects that arise from the finite charge density created when one electron is added to or removed from a finite-sized supercell, and this size dependence helps explain the success of defect level calculations using semilocal functionals.