When we use the formula to calculate two-qubit entanglement, like these:
$$ C(\rho)=\max \left\{\sqrt{e_{1}}-\sqrt{e_{2}}-\sqrt{e_{3}}-\sqrt{e_{4}}, 0\right\}\tag{18} $$
with the quantities $e_{i}\left(e_{1} \geq e_{2} \geq e_{3} \geq e_{4}\right)$ are the eigenvalues of the operator
$$ R=\rho\left(\sigma^{y} \otimes \sigma^{y}\right) \rho^{*}\left(\sigma^{y} \otimes \sigma^{y}\right),\tag{19} $$
where $\rho^*$ is the complex conjugate of the reduced density matrix $\rho$ given by Eq. (12), and $\sigma^y$ is the Pauli operator.
Why do we use the complex conjugate of the density matrix instead of its complex conjugate transpose?