Dear Sir, My intention was that there is no coupling between up and down spins. That's why I didn't apply any hopping between up and down spin.
As you said, " There is no coupling between the two spin lattices. The only hoppings you put are in the lines: sys[lat_u.neighbors()] = -1 sys[lat_d.neighbors()] = -1 " then the transmission coefficients for up spin channel (smatrix(1,0) (when there is only two leads for up spins) will be same as smatrix(2,0) (when we add both up and down spin channels). With regards, Sudin