H is the Hamiltonian. Here we set µ=0. For grand-canonical calculations add -6µ to the main diagonal.

-2*(h - U)	(-I/2)*(-I + Sqrt[3])*t
(I/2)*(I + Sqrt[3])*t	-2*(h - U)