Thank you for the answers! And how to determine the spin of the transitions from the output? As you specify "unrestricted" it should calculate both singlets and triplets (my system is in singlet (S0) initial state).