Normally, the Px, y and 3 Pz states should be degenerate, as should the 1 Px, y and 1 Pz states. However, this exact spin-contamination has produced unphysical states. The same issue was observed in SF-TDDFT. In contrast, MRSF variants (UMRSF and MRSF) eliminated the spin contamination and produced state energies of 2.652 (2.667) eV and 4.668 (4.690) eV for Px, y and 1 Px, y, respectively. While these are not identical to the Pz and 1 Pz energies, they are very close, within approximately 0.2 eV. Compared to the values from MRSF, all state energies in UMRSF are slightly reduced by about 0.01 eV, which is again attributed to the variational flexibility of UKS. It was discussed that the bulk of the slight splitting between the multiplet components is caused by the self-interaction error of the density functional.