This system is more challenging for MR methods due to the size of active space, and as a reminder, DFT predicts that all Fe centers are identical. We calculate S = 0, 1, ..., 7 spin states for each active space for all methods. )] on each fragment in the given rootspace and it's all possible spin polarizations.