Ok you were right, one of the eigenvectors is (0,-1,1,0), as you were about the ##S^2## operator; what I still do not understand is why my singlet of 0 spin is the first excited state and not the ground state, as it usually is, since the energy eigenvalue of ##(2\lambda - \sqrt{4\lambda^2...