I am working on a computational paper involving transition metal dopants. One of the reviewers asked that I explain the results in terms of the DOS. However, the results not only seem strange, but contradictory to the basic rules of crystal field theory (CFT). My Fe dopant sits in an octahedrally coordinated doping site. So according to CFT the degenerate d orbitals should split into the lower energy t2g and higher energy eg levels. However, I'm not sure how to interpret the lm decomposed projected density of states. In addition, the data seems to contradict CFT as it appears that the eg level is lower than the t2g. Any advice would be appreciated.
The results were calculated using DFT as implemented in VASP using PAW pseudopotentials. The decomposition was accomplished using LORBIT=11 I'm sure the calculation is converged.