Restricted open-shell DFT

Is there a way to run restricted open-shell DFT in Psi4? I know you’re not really supposed to run RODFT, but it would be useful to be able to initialize open-shell CAS calculations with a KS reference.

Sure, its pretty easily done if you look at libscf_solver/rhf.cc and libscf_sovler/rohf.cc. I don’t think this will be a supported feature however as in the general case this operation does not make much sense due to the needed divergence of near-nuclear spin densities.

The interesting thing is that ROKS is a partly valid option, e.g. in psi4/driver/procrouting/proc.py:

    elif core.get_option('SCF', 'REFERENCE') in ['ROHF', 'ROKS', 'UHF', 'UKS']:
        core.set_local_option('SCF', 'GUESS', 'GWH')