Effects of band edge positions on defect structure in lead halide perovskites: A case study on the Br vacancy in CsPbBr3

In practical density-functional-theory-based simulations of defects in lead halide perovskites (LHPs), it is often assumed that spin-orbit coupling (SOC) and self-interaction corrections (SIC) mostly affect the electronic structure, whereas their effects on the structural properties are minor. Therefore, SOC and SIC are usually excluded during structural relaxation or molecular dynamics simulations. With such an assumption, previous studies predicted several unusual features of Br vacancies (${V}_{\text{Br}}$) in $\mathrm{Cs}\mathrm{Pb}{\mathrm{Br}}_{3}$, including the formation of a Pb dimer, negative-$U$ character, and highly dynamic defect level. In this work, it is shown that SOC and SIC can play important roles in determining the equilibrium geometry and potential energy surface of ${V}_{\text{Br}}$ in $\mathrm{Cs}\mathrm{Pb}{\mathrm{Br}}_{3}$. Including SOC and SIC for structural relaxation results in a normal shallow level feature of ${V}_{\text{Br}}$, instead of the negative-$U$ character, and the Pb dimer becomes unstable in this case. Moreover, the highly dynamic defect level of ${V}_{\text{Br}}$ is absent when SOC is included for molecular dynamic trajectory generation. The contrasting results compared to previous studies are understood by the sensitive dependence of the occupation of the ${V}_{\text{Br}}$-associated states on the band edge positions of $\mathrm{Cs}\mathrm{Pb}{\mathrm{Br}}_{3}$, which is greatly affected by the strong SOC and SIC. These results thus highlight the importance of the correct prediction of the band edge positions for accurate modeling of defects in LHPs.