I'm following the Topological Condensed Matter online course by TU Delft. For systems with spin 1/2, time-reversal symmetry as $$ \mathcal T = i\sigma_y\mathcal K, $$ where $\sigma_y$ is the second Pauli matrix. A Hamiltonian with this kind of symmetry satisfies the equation $$ \mathcal H = \sigma_y \mathcal H^* \sigma_y, $$
Based on the authors, the following is an example of such Hamiltonian
where all the entries are real and $i$ is imaginary. I tried to compute the above equation (in Mathematica) by redefining $\sigma_y$ to be $\sigma_y\otimes\text{I}$ and calculating $(\sigma_y\otimes\text{I}) \mathcal H^* (\sigma_y\otimes\text{I})$. I obtained
which clearly is a different Hamiltonian from the original. I uploaded the Mathematica sheet to GitHub if someone wants to check it: https://github.com/ManuelGmBH/Time-Reversal-Hamiltonians.
The calculation should be straightforward, but I don't know why I don't obtain the desired result. Any help or comment is appreciated.

