How do you evaluate the canonical ensemble average of a product of spins, e.g.:
$$[S_zS_x]$$
Where:
$$S_x = \frac{\hbar}{2} \begin{pmatrix} 0 & 1\\ 1 & 0\\ \end{pmatrix}$$ $$S_y = \frac{\hbar}{2} \begin{pmatrix} 0 & -i\\ i & 0\\ \end{pmatrix}$$ $$S_z = \frac{\hbar}{2} \begin{pmatrix} 1 & 0\\ 0 & -1\\ \end{pmatrix}$$
The matrices are too simple, but then how is the ensemble average of the resultant matrix defined? EDIT: I knew it is about the density matrix, but I guess I was confuseing this with something else. Thanks anyway.