X*R
is not necessarily symmetric. So what does your use of cyclic permutation invariance of trace prove? You use of R^{1/2} seems to imply you are assuming R is symmetric PSD, which you did not state.
If X*R
were symmetric PSD, you could use the approach in Generalizing "trace_inv" for matrix quadratic forms .