for any two operators A and B and any real τ (we are working with finite-dimensional Hilbert spaces after all). We used the fact that the density matrix commutes with any function of (H − μ N) and that the trace is cyclic.