Koopmans' theorem is a useful approach to calculate the global reactivity parameters from the HOMO-LUMO energies. My question is, does it apply to the open-shell systems where we get two sets of singly occupied (alpha and beta) HOMO-LUMO energies?
I am working on rare earth systems. Rare earth (RE) elements generally show +3 oxidation states which result in open shell configurations for $\ce{Ce^3+}$, $\ce{Eu^3+}$, $\ce{Gd^3+}$ etc. I am performing unrestricted open shell DFT calculations for these systems, as restricted open shell calculations for RE systems require larger computational time and cost. In unrestricted open shell calculations, we get two sets of alpha and beta HOMO and LUMO. My question is, can I apply Koopmans' theorem here to calculate the global reactivity parameters?
Could anyone help with some literature on the reactivity parameters for open shell systems? Thank you!