“…However, while the low computational cost of semilocal functionals has very much contributed to the success of DFT because it enables access to large systems of practical relevance, the functional derivatives of typical semilocal functionals, i.e., their corresponding xc potentials, miss important features of the exact xc potential, in particular discontinuities [3,4] and step structures [5][6][7][8][9] that are relevant, e.g., in charge-transfer situations [10][11][12] and ionization processes [5,[13][14][15][16]. Many attempts have been made to incorporate some of the missing features into semilocal DFT [17][18][19][20][21][22][23][24][25][26][27]. In recent years, it was the Becke-Johnson (BJ) model potential [28] in particular that sparked interest in this respect [22,[29][30][31][32][33][34][35].…”