Matrix H in Subspace 87 with Dimension 4

H is the Hamiltonian. Here we set µ=0. For grand-canonical calculations add -6µ to the main diagonal.
2*U (-I/2)*(-I + Sqrt[3])*t Sqrt[2]*t 0
(I/2)*(I + Sqrt[3])*t 2*U 0 Sqrt[2]*t
Sqrt[2]*t 0 2*U 0
0 Sqrt[2]*t 0 3*U