Do your electrostatics / population analysis from some expensive post-HartreeFock calculation look suspiciously similar to the HF result? You probably forgot to add ‘Density=Current’ to your input stack.
I don’t understand why this isn’t the default, as when you are doing (hybrid)DFT calculations, the SCF density is being varied, thus the pop / electrostatics output change with functional.
