Suppose I have a bipartite system (with Hilbert space $H = H_a \times H_b$) and the following state:
$$\sigma = \sum_{n} \frac{e^{-\beta E_n}}{Z} \rho_n$$
where $Z = \sum_n e^{- \beta E_n}$ and $\rho_n$ are arbitrary density matrices in $H$.
I'm interested in understanding if there are entanglement measures that can separate the 'thermodynamical' entropy (i.e., not knowing in which state the system is in) and the 'quantum' entropy (i.e., measuring the correlations and entanglement between sub-systems).
So, to $\sigma$ I'd associate a 'thermodynamical' entropy $$S = - \beta \frac{\partial \log Z}{\partial \beta} + \log Z,$$ but I don't know if there's an easy way to calculate the Von Neumann entropy, for instance.
If each $\rho_n$ was $| E_n \rangle \langle E_n|$ then the Von Neumann entropy and the 'thermodynamical' entropy would coincide, but I think this is the only case they would.
So, my question is the following: is there some entanglement measure that would divide into something like classical entropy + quantum entropy?