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

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