I sometimes need the reduced density matrix function for tracing out non-normalised states, like a difference of two valid density matrices. In such a setting, re-normalising the trace to 1 does not make sense. I've added an option to turn off this behaviour by passing 'normalize = False' to the method.