Does Q-Chem have any sort of root-following procedure available during geometry optimizations with EOM methods?
We have a pretty standard case where two excited states cross, but optimizing “state #3” just leads to oscillations as these states alternate identities.
Thanks!