I try diagonalize the Hamiltonian for a 1D wire with proximity-induced superconductivity. In the case without a superconductor is all fine. However, with a superconductor I don't get the correct result for the energy spectrum of the Hamiltonian (arxiv:1302.5433)
H=\eta(k)τz+Bσ_x+αkσ_yτ_z+Δτ_x...