Matrix H in Subspace 111 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