The Hamiltonian
$$H = \left[ \begin{array}{cccc} a & 0 & 0 & -b \\ 0 & 0 & -b & 0\\ 0 & -b & 0 & 0\\ -b & 0 & 0 & -a \end{array} \right] $$
commutes with the qubit exchange operator
$$ P = \left[ \begin{array}{cccc} 1 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0\\ 0 & 1 & 0 & 0\\ 0 & 0 & 0 & 1 \end{array} \right] $$
So I would expect the two to have the same eigenvectors. The eigenvectors of $P$ are easily seen to be $(1,0,0,0)^T; (0,0,0,1)^T; (0,1,1,0)^T ; (0,1,-1,0)^T$. The latter two are also eigenvectors of $H$, but the first two are not. Why? I thought commuting operators shared the same eigenbasis?