overall ferromagnetism in this system due to the presence of multiple unpaired electrons per metal center. We point out that the canonical orbital energies remain unchanged between the CASSCF and NEVPT2 calculations, indicating that the additional dynamic correlation is captured only in the overall J coupling and not reflected in the orbital energies. This introduces a source of error when performing scans using perturbation theory methods such as CASPT2 or NEVPT2, in addition to the unaccounted 1p1h effects. Nevertheless, qualitative trends for different coupling strengths along the σ and δ paths can be reliably extracted from both CASSCF and NEVPT2 calculations.