Quantum mechanical method for the calculation of open-shell excited states. The method accurately describes low-lying multi-configurational electronic states of @[email protected] and triradicals in an efficient and robust single-reference scheme. The target @[email protected] states (e.g., S = 0 or S = 1/2) are described as spin-flipping excitations (S = −1) from a well-behaved @[email protected] reference state (e.g., S = 1 or S = 3/4). By employing theoretical models of increasing complexity for the reference (e.g., SCF, MP2, CCSD), the accuracy in the target states’ description can be systematically improved. The SF methods result in multistate single-step computational schemes, e.g., several low-lying states can be computed in a single calculation that includes both dynamical and non-dynamical correlation effects.