The accurate treatment of excited electronic states of molecules is a long-standing and active area of research in computational chemistry. In particular, the accurate treatment of excited states is difficult when a single-determinant ground state provides a poor reference for computing the excited states (e.g., double excitations or strongly correlated systems ). A useful form of wave function for overcoming such difficulties is the complete active space configuration interaction (CASCI) trial function: