Intrastate coupling constants were obtained as numerical gradients; interstate coupling constants were derived from the change in the wavefunction overlaps. Spin-orbit couplings S6 were calculated at the triplet ground-state geometry at the CASSCF level of theory. Finally, vertical excitation energies were calculated at the quartet ground-state geometry at the CASPT2 level of theory. 29