Suppose we have two spin 1/2 particles. Then we have a common eigenbasis $\{|m_{s1},m_{s2}\rangle\}$ to the operators $\hat S_{1z}$ and $\hat S_{2z}$. The permutation operator is defined as:
$$ \hat P_{21}|m_{s1},m_{s2}\rangle=|m_{s2},m_{s1}\rangle $$
My job is to find the eigenvalues of $\hat P_{21}$ and its eigenkets, as well as to prove that:
$$ \hat P_{21}=\hat 1+\frac{\hat 1}{2}+2\hat{\vec{ S_1}}\cdot\hat{\vec{ S_2}}, \tag{1} $$
Here's my attempt:
We know that for some vector $(x,y)^T$ there is a matrix that can put it like $(y,x)^T$ given by $ M= \left( {\begin{array}{cc} 0 & 1 \\ 1 & 0 \\ \end{array} } \right) $, so I would expect that, after writing the (1) in $\{|m_{s1},m_{s2}\rangle\}$ basis something of the sort would come up. But after some tweeking I get that:
$$ \hat{\vec{ S_1}}\cdot\hat{\vec{ S_2}}=\frac{\hbar^2}{4}\bigg[(\sigma_x \otimes \hat 1)(\hat 1 \otimes \sigma_x)+(\sigma_y \otimes \hat 1)(\hat 1 \otimes \sigma_y)+(\sigma_z \otimes \hat 1)(\hat 1 \otimes \sigma_z) \bigg] $$
which reduces to:
$$ \hat{\vec{ S_1}}\cdot\hat{\vec{ S_2}}=\frac{\hbar^2}{4}\bigg[(\sigma_x \otimes \sigma_x)+(\sigma_y \otimes \sigma_y)+(\sigma_z \otimes \sigma_z) \bigg] $$
As we know that $ \sigma_x \otimes \sigma_x= \left( {\begin{array}{cc} 0 & 0 & 0 & 1 \\ 0 & 0 & 1 & 0 \\ 0 & 1 & 0 & 0 \\ 1 & 0 & 0 & 0 \\ \end{array} } \right) $ , $ \sigma_y \otimes \sigma_y= \left( {\begin{array}{cc} 0 & 0 & 0 & -1 \\ 0 & 0 & 1 & 0 \\ 0 & 1 & 0 & 0 \\ -1 & 0 & 0 & 0 \\ \end{array} } \right) $ and $ \sigma_x \otimes \sigma_x= \left( {\begin{array}{cc} 1 & 0 & 0 & 0 \\ 0 & -1 & 0 & 0 \\ 0 & 0 & -1 & 0 \\ 0 & 0 & 0 & 1 \\ \end{array} } \right) $, we get that $ \sum_i \sigma_1 \otimes \sigma_1= \left( {\begin{array}{cc} 1 & 0 & 0 & 0 \\ 0 & -1 & 2 & 0 \\ 0 & 2 & -1 & 0 \\ 0 & 0 & 0 & 1 \\ \end{array} } \right) $. Plugging this all back in (1) we see that what we end up with is something like this:
$$ \hat P_{21}= \left( {\begin{array}{cc} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ \end{array} } \right) + \frac{1}{2}\left( {\begin{array}{cc} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ \end{array} } \right) + \frac{\hbar^2}{2} \left( {\begin{array}{cc} 1 & 0 & 0 & 0 \\ 0 & -1 & 2 & 0 \\ 0 & 2 & -1 & 0 \\ 0 & 0 & 0 & 1 \\ \end{array} } \right) $$
$$ \hat P_{21}= \left( {\begin{array}{cc} 3/2+\hbar^2/2 & 0 & 0 & 0 \\ 0 & 3/2-\hbar^2/2 & 3/2+\hbar^2 & 0 \\ 0 & 3/2+\hbar^2 & 3/2-\hbar^2/2 & 0 \\ 0 & 0 & 0 & 3/2+\hbar^2/2 \\ \end{array} } \right) $$
But this has no resemblance to the matrix M and tells me nothing. How can one prove (1)? Why is my approach wrong?