Yes, this is the exactly the next thing I'd like to better understand.
What is xc hfexch ? Is it DFT with Optimized Effective Potential or, equivalently, EXX ?
Then it is not HF (HF is not DFT), which leads me to the original question of how to run HF with fractional occupations.