Matrix H in Subspace 27 with Dimension 2

H is the Hamiltonian. Here we set µ=0. For grand-canonical calculations add -6µ to the main diagonal.
-2*(h - U) 0
0 -2*(h - U)